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ABSTRACT 


A  comprehensive  theoretical  model  for  ignition  of  composite  solid  propel¬ 
lants  was  proposed  as  a  framework  for  formulating  ignition  models  for  differ¬ 
ent  propellants  under  a  wide  range  of  operating  conditions. 

Detailed  chemical  kinetics  information  for  the  AP/PBAA  propellants  was 
incorporated  into  the  proposed  theoretical  model  for  ignition  of  composite 
solid  propellants,  enabling  a  more  complete  description  of  source  terms  in 
the  gas-phase  energy  and  species  equations.  Solid-phase  energy  equations 
were  coded  and  implemented  on  a  computer.  Numerical  solutions  of  the  solid- 
phase  subprogram  were  successfully  checked  against  analytic  solutions  for 
some  limiting  cases.  Using  actual  measured  heat  flux  near  the  propellant 
surface  as  an  input  to  the  solid-phase  equations,  predictions  for  the  ignition 
delay  time  were  made. 

Effects  of  pressurization  rate,  crack-gap  width,  and  igniter- flame 
temperature  on  the  ignition  of  AP-based  composite  solid  propellants  located 
at  the  tip  of  an  inert  crack  were  studied.  The  ignition  process  was  observed 
by  using  a  high-speed  camera  (-.40,000  pictures/s)  and  a  fast-response  photo¬ 
diode  system.  Heat  flux  to  the  propellant  surface  was  measured  with  a  thin- 
film  heat- flux  gage. 

A  bright  luminous  zone  behind  the  reflected  compression  wave  was  observed 
at  high  pressurization  rates.,.  It  is  believed  that  this  luminous  zone  is 
caused  by  the  combustion  of  unreacted  species  or  particles  from  the  igniter 
system.  Experimental  results  indicate  that  the  ignition  delay  time  decreases 
and  the  heat  flux  to  the  propellant  surface  increases  as  the  pressurization 
rate  is  increased.  No  distinguishable  effect  of  the  crack-gap  width  on  the 
ignition  process  was  evident.  Limited  results  obtained  by  using  an  aluminized 
propellant  as  the  igniter  show  that  the  ignition  delay  time  is  somewhat  lower 
for  higher  flame  temperature  igniter  gases.  Theoretical  predictions  employ¬ 
ing  the  solid-phase  equations  are  in  reasonable  agreement  with  experimental 
data. 

The  decrease  in  ignition  delay  with  increasing  pressurization  is  caused 
by  enhanced  heat  feedback  to  the  propellant  surface  at  higher  pressurization 
rates. --This  augmentation  in  heat  feedback  to  the  propellant  at  higher 
pres  surnTzat  ion  is  a  result  of  a  combination  of  the  following  mechanisms: 
heating  due  to  compression-wave  reflection  at  the  closed  end;  heat  transfer 
due  to  recirculating  hot  gas  near  the  tip;  heat  release  due  to  combustion  of 
unreacted  igniter  species  (ot  particles)  near  the  tip,  behind  the  compression 
wave;  and  increase  in  temperature  of  the  gas  adjacent  to  the  propellant  sur¬ 
face  due  to  continued  pressurization  of  the  crack  cavity. 
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I.  INTRODUCTION 


This  report  summarizes  progress  made  during  the  period  August  1,  1980  to 
July  31,  1981,  under  the  project  entitled  "Transient  Ignition  Mechanisms  of 
Confined  Solid  Propellants  under  Rapid  Pressurization"  (Contract  No.  N00014-79-C- 
0762). 

The  overall  objective  of  this  investigation  is  to  achieve  a  better  under* 
standing  of  ignition  and  flame-spreading  processes  of  a  solid  propellant  under 
rapidly  varying  pressure  conditions.  These  fundamental  studies  are  expected 
to  help  in  the  design  of  solid  rocket  motors  and  in  the  reduction  of  hazavds 
caused  by  convective  burning  and  deflagration-to-d.  tonation  transition  (DDT) 
processes.  The  research  includes  both  experimental  and  theoretical  studies  of 
solid  propellant  ignition.  Specific  objectives  of  this  project  are: 

1.  To  formulate  a  comprehensive  ignition  model  for  composite  solid 
propellants; 

2.  To  observe  the  detailed  ignition  phenomenon  and  to  measure  the 
ignition  delay  and  instantaneous  heat  flux  to  the  propellant  sample 
surface ; 

3.  To  study  the  effect  of  pressurization  rate  (dP/dt),  igniter  gas 
temperature  (T£) ,  and  propellant  type  on  ignition  delay;  and 

4.  To  predict  the  ignition  delay,  using  measured  values  of  heat 
flux  to  the  propellant  sample  surface,  and  to  compare  calculated 
and  measured  values. 

Prior  to  formulating  the  ignition  model,  an  extensive  literature  survey 

was  conducted.  Details  of  this  survey  were  included  in  our  last  annual  report,1 

2 

and  were  also  presented  at  the  AIAA  16th  Joint  Propulsion  Conference.  Theoretical 
and  experimental  work  performed,  and  major  results  obtained  during  this  reporting 
period,  appear  in  Sections  II  to  IV. 
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II.  THEO?ETICAL  WORK 


2.1  Ignition  Model  for  AP/PBAA  Propellant 

A  comprehensive  model  was  formulated  in  order  to  study  ignition  of  the  composite 
AP/PBAA  propellant  under  rapid  pressurization  conditions.  The  detailed  mathematical 
description  of  this  model  was  reported  in  the  last  annual  report.*  A  brief  summary 
of  the  formulation  of  the  model  is  given  below  to  facilitate  discernment  of  further 
developments  and  solution  procedures. 

The  model  describes  the  ignition  process  for  an  AP-based  composite  propellant. 

To  make  the  problem  mathematically  tractable,  a  typical  oxidizer  crystal  is  considered 
to  be  embedded  in  a  fuel  binder  in  the  configuration  shown  in  Fig.  1.  Dimensions  of 
the  fuel  and  oxidizer  are  determined  by  statistical  averaging  of  the  bulk  of  com¬ 
posite  AP/PBAA  propellant.  The  model  considers  two-dimensional  axisymmetric  geometry, 
includes  the  effect  of  sudden  pressurization  in  the  surrounding  region,  and  makes 
no  a  pnLonL  assumption  concerning  the  site  of  ignition.  The  model  allows  ignition 
to  occu;  in  any  region  (gas  phase,  solid  phase,  or  interface)  by  considering  surface, 
subsurface,  and  gas-phase  chemical  reactions.  Governing  heat  equations  for  the  oxi¬ 
dizer  and  the  fuel  binder  in  the  solid  phase  are 


Oxidizer:  o  c  a*?  -,2_ 

P0x,sC0x,8  -  IK  -2-4  +  k  i-I  +  -  ill  + 

Ox,s  3z2  kOx,s  3r2  r  9r]  +  q0x. 


Fuel: 


ax 


&  +  K  a  +  i- 


3  —  *  rk_  -  A  +  Ir  Cf,  .1  ,  t  ,8  ST, 

F,SCF,S  at  [kF,3  3z2  +  kr,s^2+  r  3?1  +  qF,s 


(1) 


(2) 


where  q'"  - 


Ox,s  q0x, radiation  +  q0x, photochemical  +  q0x, pyrolysis 


km 


Ox,  relation  Ox  p0x 


*  T0_GB_e  Sox  *,  S  i 

Ox  z 


C3) 


*0x,  photochemical 


T0xIzQ0x,pcAn0xxdX  (see  Re£-  9) 


(5) 


^  °0x, pyrolysis  "  “  Z0x,py  Q0x,py  ^^Ox.py^V^  ^See  Re£*  95  (6) 


Fig.  1  Cross-sectional  View  of  a  Statistically 
Averaged  Element  Considered  in  the  Model 


. .  .  *:.  r h  J  >•• 


Initial  and  boundary  conditions  for  the  solid  phase  (fuel  and  oxidi2er)  are 
at  t  *  0  :  T(0,r,z)  « 


on  r  =  0 


on  r  *  R, 


§-° 


at  the  solid-gas  interface,  z  »  z  : 

s-g 

‘SL  ♦  ■  kg  Ir  Z  -  *  “s-gG  -  SEb 

s-g  *  s-g  5  S 


+  r.  P  T(c  -  c  1  +  q" 
V«  s  pgJ  Vg 


at  z  =  -  L  : 

-  k  II 

kF,s  3z 


3T 

“  kOx,s  37  _  L+  +  ^x-F  +  a0x-FT0xGe 


“&0xL 


T  ,  -  ~  T 

z  =  -  L  „  .+ 

z  *  -  L 


at  r  =  R, 


CF,s  3r  R1+ 


Q  l  ,f 

0  x,s  3F  R~  \)x-F 


r  =  R, 


r  *  R, 


at  z  -*■  -  oo  ;  T  *  T. 

i 


Governing  equations  for  the  gas  phase  are 


Continuity  Eq: 


.  8(P.V 

at  a  z 


Energy  Eq: 


OT  ,  3T 

Cp  -T-+PCV  — 

p  g  3t  g  p  z  3z 


JP  „  I  JL  ev  ~  ill  x  3  /t,  3Tv  J.  • , 
3t  r  3r  ^kg  arJ  +  3z  (kg  Sz^  +  qc 


mss* 


1 


Species  Eq: 

3Y 


j 


3Y, 


■'g  St  +  pgvz  3z 


i  a  3Y< 

i  A  (rpp8  ir> +  «T 


(18) 


1 1 


i 

i 


where  j  ■  1,  2,  3,  or  4  (for  the  gas  phase  species:  Oxidizer, 

NH^ ,  HCIO^,  and  Fuel) 


Equation  of  State: 


Ps  -  PgCP.T) 


(19) 


The  momentum  equation  is  replaced  by  a  measured  pressure-time  trace  near  the 
propellant  surface. 


Initial 

and  boundary 

conditions  for  the  gas  phase  are 

at 

t  -  0 

* 

V  «  V 

z  z,  i 

(20) 

at 

t  -  0 

• 

T(0,r,z)  *  T  . 

g»i 

(21) 

at 

t  -  0 

« 

• 

Yj  "  Yj,i 

(22) 

on 

r  -  0 

|S.  o 

3r 

(23) 

on 

r  »  0 

: 

Tr  "  0 

(24) 

on 

r  -  R2 

: 

(25) 

on 

r  «  R2 

: 

O 

« 

s1« 

(26) 

on 

z  ■  z 

a-g 

: 

P  v  -  p  r, 

g  z  s  bg 

(27) 

on 

z  *  z 

s-g 

: 

Tj  -  T 

z  +  z 

1  3-g  S-g 

(28) 

on 

z  -  z 

S-g 

• 

v*v) .  •  V»y3  ,  -  °8 1  ^  -  *1 

3-g~  S-g+ 

(29) 

on 

z  -*•  00 

• 

T  -  Tjt) 

(30) 

2 . 2  Reaction  Kinetics  of  AP/P8AA  Propellants 

Five  chemical  reactions  are  considered  for  this  family  of  propellants,  viz., 
exothermic  degradation  of  solid  AP  particle  into  oxidizer  gases;  dissociative 
sublimation  of  AP  into  NHj  and  HCIO^;  reaction  between  premixed  NHj  and  HC104 
gases;  fuel  pyrolysis;  and  diffusion  flame  resulting  from  chemical  reaction 
between  oxidizer  and  fuel  gases  to  form  products.  Detailed  reaction  steps  for 
tho  AP/PBAA  propellants  are  given  in  Table  1.  Thu  sixth  reaction  shown  at  the 
bottom  of  Table  1  represents  the  heterogeneous  reaction  between  the  fuel  vapor 
and  solid  oxidizer.  However,  this  reaction  is  not  considered  in  the  current  model 

because  no  data  is  presently  available. 

Table  1  Reaction  Steps  for  AP/PBAA  Propellant 


Chemical  Reactions 
Considered  in  Mathematical  Model 


AP 


2.  AP 


(s)  Exothermic 
Degradation 


Ox 


(8) 


(s)  Dissociative  3 
Sublimation 


NH,  +  HC2.0. 


(8) 


(8) 


3.  NH  +  HCJIO, 
J(g) 


(g) 


P remixed 
A/ PA  Reaction 


Ox 


(8) 


4.  PBAA 


(s)  Pyrolysis  Fuel(g) 


°X(g)  +  Fuel(g)  Diffusion 
_ Flame 


Products 


(8) 


Reaction 
Appears  in 


B.C. 


B.C. 


Gas  Phase  Eq. 


B.C. 


Gas  Phase  Eq. 


6.  Fuel ,  ■.  +  AP.  .  ;;  -  Products 

(g)  (s)  Heterogeneous 

Reaction 


B.C. 
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A  specific  reaction  rate  k  is  given  by  the  Arrhenius  expression 

-E/RT 

k-ZTe  (31) 

Typical  values  for  the  above  chemical  kinetic  constants,  together  with  the 
sources,  are  presented  in  Table  2. 

2. 3  Numerical  Solution  of  a  Simplified  Case 

A  simplified  case  of  the  generalised  model  described  above  was  solved  initially. 
Instead  of  using  calculated  heat  feedback  from  the  gas-phase  equations,  measured 
heat  flux  to  the  propellant  surface  was  used  as  an  input  to  the  solid-phase  heat 
equation. 

There  are  several  advantages  in  obtaining  a  partial  solution  to  the  generalized 
model.  Since  the  heat  flux  is  measured  at  a  location  adjacent  tp  the  propellant 
surface,  it  provides  realistic  values  of  heat  feedback  to  the  propellant  surface, 
even  for  such  a  complex  physical  situation.  The  simplified  model,  that  considers 
only  solid-phase  equations  and  employs  measured  heat  feedback  to  the  propellant 
surface,  will  yield  reasonably  accurate  values  of  ignition  delays.  Because 
solutions  to  several  limiting  cases  of  the  solid  phase  equations  are  well  docu¬ 
mented,  results  obtained  from  the  simplified  model  will  make  it  simpler  to 
independently  check  the  numerical  solution  procedure  for  the  solid-phase  portion 
of  the  complete  model.  The  detailed  solution  metuod  for  the  solid-phase  governing 
equations,  and  comparisons  with  some  analytic  results  to  limiting  cases,  are 
given  in  the  following  subsections.  It  should  be  noted  that  even  though  the 
major  advantage  of  the  partial  solution  is  its  simplicity,  results  thus  obtained 
depend  upon  experimental  determination  of  the  heat  flux,  which  is  expensive 
and  not  always  feasible.  Therefore,  solution  of  the  complete  model  (both 
solid-  and  gas-phase  equations)  is  essential  for  detailed  prediction  of  the 
ignition  process. 
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2.3.1  Numerical  Scheme 

A  finite  difference  scheme  was  used  for  the  numerical  solution  of  solid-phase 
equations.  A  schematic  of  the  grid  system  in  r  and  £  coordinates  is  shown  in 
Fig.  2.  The  following  features  have  been  incorporated  into  the  numerical  solu¬ 
tion  procedure  for  improved  accuracy,  faster  convergence,  and  reduced  computation 
costs. 

i)  Variable  mesh  size 

ii)  Central  difference  scheme  to  approximate  the  spatir?  derivatives 

iii)  Generalized  Crank-Nicolson  representation  of  the  inite-difference 
equations 

iv)  Quasi linearization  of  the  inhomogeneous  or  source  terms. 

A  successive  overrelaxation  iterative  scheme  was  used  to  solve  the  resulting  set 
of  simultaneous  algebraic  equations. 

Figure  3  shows  a  stencil  diagram  for  a  general  node  (I,  J,  K) .  As  indicated 
by  the  arrows,  calculation  of  variables  at  the  new  time  step,  i.e.,  at  (I,  J, 

K+l),  is  influenced  by  all  nodes  surrounding  the  node  (I,  J,  K+l),  as  well  as 
those  surrounding  the  node  (I,  J,  K) .  Finite  difference  expressions  for  various 
derivatives  are  given  below.  Here,  9  is  the  Crank-Nicolson  parameter. 


3T  3T 

It  "  at 


3^T  _  e^T 

as2  a?2 


W  »  V  1 

~2  7T 
ar  (jr 


-K+l  ,  tK+1  _K+1  tK  .  _K  tK 

e  JILL  1 1  L.J _ 1L LI  .(i.e)  JILLLLIlI _ Jihl 


(A Sj)' 


(A5j)‘ 


K+1  -K+l  K+l  K  _  K  K 

0  JjJJL  "  l,J  I.J-1  .  (1  .  0)  JLiili  2  I.J  I.J-1 


(Arj)' 


(Arj)‘ 


r  K+l  JC+1 

1  3T  1  3T  °  1  „  ‘l.J+1  '  1 I.J-1 

r  IF  F  3F  ’  FJ  *  24*, 


♦  (1  -  8) 


TK  -  tk  * 
I.J+l  *I,J-1 


OXIDIZER 


ADIABATIC  WALLS  £(I) 


Fig.  2  Finite  Difference  Grid  for  Solid  Phase 
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2.3.2  Procedures  Used  In  Checking  Numerical  Solution 

Three  major  checks  were  made  on  computed  solutions.  First,  numerical  results 
were  compared  with  the  analytic  solution  for  the  simple  case  of  a  semi-infinite 
homogeneous  solid  at  uniform  initial  temperature,  with  constant  external  heat 
flux  starting  st  t  ■  0.  Conparison  of  calculated  and  analytic  distribution  of 
temperature  within  the  solid  at  various  times  was  excellent  (see  Fig.  4a).  Surface 
temperature-time  history  values  were  also  very  close  (see  Fig.  4b).  In  the  second 
check,  a  global -energy  balance  was  made.  The  total  energy  input  from  outside 
at  various  time  steps  was  within  the  acceptable  error  limits  of  the  net  increase 
of  internal  energy  obtained  from  computed  temperature  distributions.  The  third 
check  involved  energy  balance  in  elemental  volumes  around  some  arbitrarily  chosen 
nodes.  The  local  energy  balance  also  showed  excellent  agreement. 


2.4  Solution  of  the  Complete  Model 

Governing  equations,  as  outlined  in  section  2.1,  represent  a  set  of  transient., 
second-order,  coupled,  inhomogeneous,  nonlinear  partial  differential  equations. 

An  exact  analytic  solution  of  this  set  of  equations  is  not  possible,  and  approxi¬ 
mate  solution  methods  employing  asymptotic  expansion  or  similarity  do  not  appear 
to  be  very  useful.  For  a  complete  solution,  therefore,  the  numerical  method 
is  used.  Several  key  steps  of  the  solution  procedure  ire  briefly  described  in 
the  following  subsections. 

2.4.1  Coordinate  Transformation 

Since  the  mathematical  domain  in  z -direct ion  is  infinite,  it  was  transformed 
into  a  finite  region  using  the  transformation  given  below.  This  transformation 
also  makes  possible  a  finer  grid  size  near  the  propellant  surface  in  the  real 
spatial  coordinate,  where  large  temperature  gradients  exist,  while  using  a  uniform 
grid  in  the  transformed  coordinate.  This  transformation  alleviates  the  instability 


Axial  Distance,  yra 


Fig.  4a  Analytical  and  Numerical  Temperature  Distribution 
in  an  Inert  Solid  with  a  Constant  Heat  Flux 


Temperature , 
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Fig,  4b  Comparison  of  Analytical  and  Numerical 
Surface  Temperature -Time  History 
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associated  with  nonunifora  grid  spacing,  which  would  have  been  necessary  for  the 
z  coordinate.  The  following  exponential  coordinate  transformation  is  employed 
in  the  gas  phase. 


This  gives 


s  ■  1  -  exp  (-Ax) 
or,  z  *  -  ~  In  (1  -  s) 

z  »  *  :=>  s  »  1 

z  ■  0  — >  s  ■  0 


(36) 

(37) 

(38) 


lz  *  A  (1  ’  S)  I? 


(39) 


a2  .2  >2  a2  .2  %  a 

3z2  3s2  dS 

Note:  The  r  axis  remains  unchanged. 

The  transformed,  gas-phase,  governing  equations  t 
Continuity  Eq. 

3(Pa  v  ) 

.  A  Cl  -  s)  — fc-2-  -  0 


Energy  Eq. 


3T 


3T  3P 


pg  £  *  pgViA(1  ’  '  IF 


P  8 


8  P 

c  ’ 

8 


(40) 


(41) 


■ft?  (kgr  * A  (l  *  s)  h  [V  C1  • s)  It1  *  * 


Species  Eq.  (j  *  1,  2,  3,  4) 

3Y, 


IN 

8 

(42) 


3Y. 


"i  13  v*i 

*  3t  *  Vl*  C»  -  *>  Si2-  ■  r  5?  (r%  3T>  *  “j  (43) 


All  boundary  conditions  remain  unchanged  except  the  mass  balance  of  individual 
species  at  the  interface  eq.  (29): 


°n  5  *  Ss-g  :  PgVj 


pgvzYj 


s-8 


3Y. 

s  ♦  -  va  {1  • s)  tt  -  ay  (44) 

s-g  •  J 
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2.4.2  Modeling  of  Source  Terms 

•  • 

In  order  to  solve  the  governing  equations,  source  terms  (q",  qm,  (d'.*,  and  w'.'O 

SI  J 

must  be  expressed  as  functions  of  other  variables  and  known  parameters.  In 
modeling  the  source  terms,  it  is  implicitly  assumed  that  the  rate-controlling 
reactions  are  single-step  and  irreversible.  It  is  further  assumed  that  the 
reaction  rate  constant  has  Arrhenius  dependence  on  local  temperature. 

If  the  gas-phase  mass  fractions,  Y^  Y2,  Yg,  Y^,  Yg,  represent  the  oxidizer, 
ammonia,  perchloric  acid,  fuel,  and  products,  respectively,  the  source  terms  can 
be  expressed  as  follows. 


“i"  <-VR„T) 

-  Z5  S^g  Y1Y4  '’*0  (-W) 
“7*  -  Z3  Y2Y3T«I>  (-E3/RuT> 

“3- “2  Wj 

“4--Z5^»g  Y1Y4  e!t^  (-W> 
*^g  a  -  (w2*AH3  +w^'AH5) 


(45) 

(46) 

(47) 

(48) 

(49) 


i/here  AHgis  the  heat  of  reaction  per  unit  mass  of  ammonia  for  reaction  3,  and.  AHg 
is  the  heat  of  reaction  per  unit  mass  of  fuel  for  reaction  S.  Other  source  terms 
are  evaluated  in  a  similar  fashion.  It  should  be  noted  that  the  reaction  rate 
expression  for  reaction  3  is  taken  from  Guirao  and  Williams.*** 

A  quasilinearization  technique  was  used  to  linearize  the  inhomogeneous  (or 
source)  terras  of  the  governing  equations.  To  illustrate  the  quasilinearization 
procedure,  F(T,  Y.j  represents  an  inhomogeneous  term  of  the  governing  equations. 
Using  the  Taylor  series  expansion  technique,  F  at  any  time  k  +  9  is  evaluated  in 


13 


the  following  manner: 

Fk+9  -  Fk  ♦  (Tk+  0  .  Tk)  k  +  (Y  k+6  -  Y.k)  |£-  k  (50) 

1  J  0  j 

L  i  0  J»  ^  1  1^ 

where  for  any  variable  U,  U  *  6U  +  (1-9)  U  .  The  partial  derivatives  of  F 
are  obtained  analytically. 

It  should  be  noted  that  for  the  five  species  considered  in  the  model,  only 
three  species  equations  have  to  be  solved  independently.  This  is  the  case  since 
species  2  and  3,  i.e.,  ammonia  and  perchloric  acid,  can  be  shown  to  have  identical 
governing  equation  and  boundary  conditions.  Mass  fraction  of  products  Yg  is 

obtained  from  the  following  algebraic  relation. 

4 

Ys  ■  1  -  2  Y.  (SI) 

j  =  l  J 

The  finite  difference  representation  of  the  gas-phase  equations  is  very 
similar  to  that  of  the  solid-phase  equations.  As  noted  earlier,  the  coordinate 
transformation  in  the  gas  phase  allows  one  to  use  uniform  grid  size  in  the  trans¬ 
formed  coordinates,  thereby  bypassing  some  of  the  instabilities  associated  with 
nonuniform  grid  spacing.  A  general  layout  of  the  computation  procedure  is  given 
in  Fig.  5.  First,  a  temperature  and  specios  mass-fraction  distribution  is 
assumed  in  the  entire  region  for  the  new  time  step.  Since  the  pressure  is  known, 
the  density  is  easily  computed  from  equation  of  state  using  the  assumed  temperatures. 
The  continuity  equation  is  then  integrated  to  obtain  v  at  each  node  point.  The 
gas -phase  energy  equation  is  then  solved  using  an  iterative  procedure.  Because 
the  species  equation  for  Y^  is  independent  of  the  mass  fractions  of  other  species, 
it  is  integrated  first  to  obtain  the  new  value  of  Y^.  Species  equations  for  Yj 
and  Y4  are  then  solved  iteratively  to  obtain  mass  fractions  of  fuel  and  oxidizer. 
Using  the  new  values  of  mass  fractions,  the  energy  equation  is  solved  again  to 
obtain  gas-phase  temperature,  until  both  temperature  and  mass  fractions  converge. 
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From  the  gas-phase  temperature  distribution,  heat  flux  to  the  propellant  surface 
is  computed,  and  the  solid-phase  heat  equations  are  solved.  The  solid-gas  interface 
temperatures  thus  calculated  are  used  as  the  new  boundary  conditions  to  solve  the 
gas-phase  equations  again.  This  procedure  is  repeated  until  the  solid-gas 


interface  temperatures  at  a  given  time  step  converge. 
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III.  EXPERIMENTAL  WORK 


3.1  Test  Apparatus 

The  test  apparatus  was  designed  to  obtain  high -temperature,  high-pressure 
gases,  which  closely  simulate  conditions  in  an  actual  rocket  motor.  Pressuriza¬ 
tion  rates  obtained  in  this  setup  are  close  to  those  encountered  before  the  onset 
of  transition  to  detonation.  A  solid-propellant  igniter  system  was  used  to  obtain 
the  hot-gas  composition  similar  to  that  in  actual  rocket  motors.  Using  this 
apparatus,  it  is  possible  to  obtain  pressures  up  to  SO  MPa  (-7,000  psi) ,  and 
pressurization  rates  in  the  order  of  100  GPa/s  (-106  atm/s). 

Figure  6  shows  a  schematic  of  the  igniter  system,  consisting  of  electric 
primer  (FA  874)  and  its  housing  unit;  booster  propellant  and  booster  nozzle; 
solid  propellant  igniter  charge  and  igniter  chamber;  and  mult iperf orated  exit 
nozzle.  Exploded  view  of  the  igniter  chamber  is  shown  in  Fig.  7.  When  voltage 
is  applied  to  the  electric  primer,  it  produces  hot  gases  which  flow  over  the 
propellant  strip  and  ignite  it.  The  product  gases  flow  through  a  multiperforated 
nozzle  into  the  main  chamber. 

A  schematic  of  the  test  section  is  shown  in  Fig.  8.  A  crack-like  cavity, 
formed  between  a  transite  slab  and  a  plexiglass  window,  is  situated  normal  to 
the  flow  direction  of  gases  in  the  main  chamber.  Thickness  of  the  interchange¬ 
able  metal  spacers  determines  the  crack-gap  width.  In  the  present  study,  the 
gap  width  was  varied  between  1.17  mm  and  3.18  mm.  The  length  of  the  cavity  is 
145  mm.  As  shown  in  Fig.  8,  propellant  test  sample  is  situated  at  the  tip 
of  this  cavity.  A  heat-flux  gage  is  symmetrically  placed  adjacent  to  the 
propellant  sample  at  the  tip.  Two  pressure  ports  are  provided  in  the  crack  cavity 
one  near  the  test  sample,  and  the  other  at  the  crack  entrance.  Exploded  view  of 
the  test  section  is  presented  in  Fig.  9.  Figure  10  shows  the  partially  assembled 
test  section  in  the  chamber. 
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Fig.  6  Schematic  Diagram  of  Igniter  System 
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B  Primer  Housing  Unit 
C  Electric  Primer 
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E  Dielectric  Material 


F  Booster  Chamber 
G  Booster  Nozzle 
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Fig.  7  Exploded  View  of  Igniter  System 
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A  Retainer  Block 
B  Heat  Flux  Gage 

C  Propellant  Sample  and  Heat-Flux 
Gage  Housing  Block 

D  Propellant  Sample  and  its  Holder 


E  Pressure  Transducer  Housing  Unit 
F  Transite  Plate 
G  Metal  Spacers 
H  Plexiglass  Window 


Fig.  9  Exploded  View  of  Test  Section 
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The  ignition  event  is  observed  through  a  set  of  transparent  plexiglass  windows 
The  window  assembly  is  held  in  place  by  the  window  retainer,  which  forms  the 
top  half  of  the  chamber.  With  the  exception  of  the  exit  nozzle  through  which 
gases  are  discharged  into  the  atmosphere,  the  chamber  is  completely  sealed  during 
tests.  For  safe  operation,  the  chamber  is  also  equipped  with  a  port  for 
the  burst  diaphragm;  the  diaphragm  ruptures  at  a  pressure  of  about  70  MPa  (~1 0,000 
psi).  Figure  11  shows  the  assembled  test  chamber,  along  with  the  igniter  system. 

Several  parameters  can  be  varied  to  alter  chamber  pressurization  and  igniter- 
gas  temperature.  The  rate  of  chamber  pressurization  is  varied  by  altering  any 
one  or  more  of  the  following:  a)  mass  of  the  propellant  igniter  charge,  b)  mass 
of  the  booster  propellant,  c)  chamber  volume,  d)  insertion  of  a  bursting  dia¬ 
phragm  between  the  igniter  chamber  and  the  test  chamber,  and  e)  diameters  of 
exit,  booster,  and  multiperforated  igniter  nozzles.  In  order  to  change  the 
igniter  gas  temperature,  igniter  propellants  with  different  flame  temperatures 
were  used. 

3 . 2  Data  Acquisition  System 

A  block  diagram  of  the  data  acquisition  system  used  in  the  present  study  is 
shown  in  Fig.  12,  This  system  is  comprised  of  1)  pressure  measuring  system, 

2)  heat-flux  measuring  system,  3)  photodiode  setup,  4)  high-speed  photography 
system,  and  5)  transient  wave-form  recording  system.  The  ignition  event  was 
remotely  controlled,  using  a  relay  circuit  connected  to  the  footage-controlled 
event  switch  on  the  camera. 

Pressure  measurements  were  made  at  three  locations  in  the  chamber:  one  each 
near  the  propellant  sample,  at  the  entrance  of  the  crack  cavity,  and  in  the  main 
chamber.  Piezoelectric  quartz  transducers  with  a  rise  time  of  1.5  us  and  natural 
frequency  of  300  kHz  were  used  to  measure  the  pressure.  Signals  from  the 
transducers  were  amplified  through  a  charge  amplifier  before  recording.  Heat-flux 
measurements  were  made  with  thin- film  heat -flux  gages.  The  heat -flux  gage  is 
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Fig.  12  Block  Diagram  of  Data  Acquisition  System 
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essentially  a  thermocouple  in  which  0. 5-1.0  ym  thick  films  of  thermoelements  are 

deposited  on  a  ceramic  base  (mullite).  The  response  time  of  these  thin-film 
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thermocouples  is  in  microseconds.  In  the  present  study,  platinum  and  rhodium 
were  used  as  thermoelements.  Output  of  the  thermocouple  is  amplified  by  a  wide¬ 
band  amplifier.  The  procedure  for  obtaining  the  heat  flux  from  temperature 
measurements  is  outlined  in  a  later  subsection. 

The  photodiode  and  high-speed  camera  systems  were  used  simultaneously  to 
monitor  the  ignition  event.  A  16-mm  Hycam  movie  camera  recorded  the  complete 
ignition  event  in  the  chamber.  With  a  quarter-frame  optical  head,  the  camera 
was  capable  of  filming  at  about  40,000  pictures/s.  The  photodiode  system  consists 
of  a  photography  lens  mounted  on  an  optical  bed,  a  fast-response  photodiode, 
and  an  amplifier.  In  order  to  observe  ignition  of  the  propellant  surface,  a  small 
region  around  the  propellant  surface  is  focused  onto  the  photodiode  !>y  the  lens. 
Output  of  the  photodiode  can  be  amplified  before  recording.  Use  of  these  two 
independent  systems  for  the  observation  of  the  ignition  event  reduces  the  possi¬ 
bility  of  error.  Figure  13  shows  the  assembled  test  chamber  and  igniter  system, 
the  Hycam  caraora,  and  the  photodiode  system. 

All  transient  signals  are  recorded  on  a  digital  transient  wave-form  recorder 
(Physical  Data  Model  515),  as  well  as  on  a  high-frequency  FM  tape  recorder.  The 
transient  wave-form  recorder  has  a  sampling  rate  up  to  2  x  106  samples/s,  a 
maximum  amplitude  resolution  of  0.1%,  and  a  4k  word  memory  per  channel.  A  light- 
emitting  diode  (LED)  driver  unit  is  used  to  generate  a  trigger  or  common-time 
pulse,  which  is  simultaneously  recorded  on  the  film  and  on  the  transient  recorders. 
This  common-time  pulse  is  used  for  time  correlation  between  data  on  the  film  and 
on  the  transient  recording  devices.  Digitized  output  of  the  transient  recorder 
can  be  displayed  on  an  oscilloscope  or  an  x-y  plotter,  or  can  be  transferred  to 
the  miniconqputer  (POP  11/23)  for  data  storage  and  processing. 
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3.3  Determination  of  Surface  Heat  Flux 

The  procedure  for  determining  heat  flux  from  the  thin-film-gage  temperature 

data  is  similar  to  that  outlined  by  Vidal  et  al . 20,21  Since  the  thickness  of 

the  thermoelements  is  less  than  1  pm,  temperature  measured  by  the  thermocouple 

is  assumed  to  be  the  surface  temperature  of  the  insulating  ceramic  base.  3ecause 

the  duration  of  the  event  is  in  the  order  of  1  ms,  the  ceramic  base  is  treated 

as  a  semi-infinite  solid  for  heat-flux  computations.  Tho  expression  for  heat 
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flux  for  a  semi-infinite  solid  with  known  surface  temperature-time  history  is 
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where  T(t)  is  the  measured  surface  temperature. 

The  integral  appearing  in  Eq.  (S2)  must  be  evaluated  numerically.  A  simple 

trapezoidal  rule  or  Simpson's  rule  techniques  do  not  provide  accurate  answers, 

since  the  integrand  becomes  infinite  at  the  upper  limit.  A  method  proposed  by 
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Cook  and  Felderman'"  was  employed  to  alleviate  this  problem.  T(t)  is  approximated 
by  a  piecewise  linear  function.  The  final  expression  for  heat  flux  is  given  by 
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where  At  =  t  /n. 

n 

The  steps  involved  in  computing  heat  flux  from  voltage-time  data  is  as  follows. 
Voltage-time  data  is  converted  into  temperature-time  data  using  tabulated  temperature- 
EMF  values  for  platinum- rhodium  thermocouples.  The  temperature-time  curve  is  then 
smoothed,  using  a  cubic  spline  fit  through  the  data.  This  smoothed  temperature 
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time  information  is  used  in  Eq.  (53)  to  compute  the  surface  heat  flux  as  a  function 
of  time.  Properties  of  mullite,  which  are  a  function  of  temperature  as  shown  in 
Fig.  14,  are  averaged  over  the  temperature  range. 


Specil 
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IV.  RESULTS  AND  DISCUSSION 

The  ignition  of  an  AP-based  composite  solid  propellant  with  PBAA-EPON  binder 
was  a  subject  of  this  investigation.  The  weight  percent  of  oxidizer  in  the  propellant 
was  75%,  and  the  average  oxidizer  particle  size  was  76  pm.  The  flame  temperature 
of  this  propellant  is  1920  K.  In  most  of  the  experiments,  the  same  propellant  was 
also  used  as  the  igniter.  In  all  tests,  initial  pressure  in  the  chamber  was  0.1 
MPa  and  initial  temperature  was  about  298  K. 

The  sequence  of  events  that  occur  during  the  test  can  be  described  generally 
as  follows.  Hot,  combustion  product  gases  from  the  igniter  flow  into  the  main 
chamber  and  pressurize  it.  The  pressure  gradient  causes  a  part  of  the  gases  to 
penetrate  the  crack-like  cavity.  Hot  gases  and  the  pressure  front  propagate  along 
the  crack,  reach  the  tip,  and  are  reflected  from  the  closed  end.  As  the  process 
continues,  additional  hot  gases  are  driven  into  the  crack  because  of  the  continually 
rising  chamber  pressure.  As  the  hot  gases  reach  the  propellant  surface,  energy  is 
transferred  from  hot  gases  to  propellant.  Heat  feedback  from  the  gases  can  be 
extremely  high  because  of  compression  wave  reflection  from  the  tip;  flow  recircu¬ 
lation  near  the  closed  end;  heat  release  from  the  unreacted  igniter  species  (or 
particles)  behind  the  reflected  compression  wave;  and  increased  enthalpy  of  the 
hot  gases  because  of  continued  pressurization  of  the  cavity.  This  heating  eventually 
causes  the  propellant  to  ignite. 

In  the  experimental  portion  of  this  study,  ignition  is  defined  as  the  onset 
of  emission  of  luminous  light  from  the  propellant  surface.  A  go/no-go  type  of 
ignition  criterion  cannot  be  applied  to  the  test  configuration  of  this  study  since 
it  was  impossible  to  remove  the  energy  stimulus  following  onset  of  ignition. 

Although  the  criterion  of  onset  of  luminous  light  is  usually  associated  with  onset  of 

surface  ablation,  especially  for  nitramine  propellants  where  considerable  time 

lag  exists  between  surface  ablation  and  ignition,  the  ignition  criterion  employed  in 
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this  study  is  quite  suitable  for  AP-based  composite  propellants,  and  it  has  been  widely 
used  in  the  past.  In  all  but  one  of  approximately  sixty  tests  conducted  in  this 
study,  the  propellant  sample  burned  completely.  Figure  15  is  a  photograph  taken 
following  testing  of  the  burned  plexiglass  window  next  to  the  propellant  sample. 

As  mentioned  in  the  previous  section,  pressure  measurements  were  made  at  two 
locations  in  the  cavity,  one  near  the  propellant  sample  at  the  tip  and  another 
near  the  crack  entrance.  Heat  flux  from  the  gases  to  the  propellant  was  deduced 
from  a  heat-flux  gage  situated  adjacent  to  the  test  sample.  The  ignition  event 
was  simultaneously  recorded  on  film  and  a  photodiode  focused  on  the  propellant 
surface.  High-speed  photography  was  conducted  by  using  either  a  streak  attachment 
or  a  quarter-frame  optical  head.  Maximum  filming  speed  during  the  tests  was  about 
40,000  pictures/s.  It  was  noted  that  light  intensity  obtained  through  streak  record¬ 
ing  near  the  tip  was  a  replica  of  the  photodiode  signal.  Therefore,  in  order  to 
achieve  a  more  complete  understanding  of  the  ignition  event,  the  quarter-frame  head 
was  used  in  most  of  the  tests.  Even  though  the  time  interval  between  each  picture  is 
about  25-30  us,  ignition  delay  can  be  accurately  determined  by  first  locating  the 
ignition  event  on  the  film  and  then  obtaining  its  more  precise  location  from  the 
photodiode  signal. 

A  typical  set  of  time-correlated  pressure  traces  at  the  crack  entrance  and  at 
the  tip  is  presented  in  Fig.  16.  The  curves  shown  pertain  only  to  the  uprising  part 
of  the  P-t  trace  because  that  is  the  region  of  interest.  Figure  16  shows  that  the 
first  discernible  pressure  rise  at  the  tip  occurs  some  time  later  than  the  corresponding 
pressure  rise  at  the  entrance,  due  to  the  finite  time  required  by  the  pressure  front 
to  propagate  from  crack  entrance  to  tip.  Also,  except  for  a  short  period  immediately 
following  arrival  of  the  pressure  front  at  the  tip,  the  average  pressurization  rate 
is  generally  lower  near  the  tip,  perhaps  as  a  result  of  frictional  losses  in  the  channel 
as  gases  travel  to  the  closed  end.  The  pressurization  part  of  this  curve  is  quite 
linear.  Ignition  was  observed  to  occur  during  the  uprising  part  of  P-t  traces. 


Another  tine-correlated  pressure  trace,  for  which  the  initial  pressurization  rate  was 
extremely  high  (dP/dt^)  *1260  GPa/s,  is  shown  in  Fig.  17.  From  the  P-t  trace  at 
the  tip  in  Fig.  17,  a  very  steep  shock-like  compression  wave  is  clearly  evident  at 
the  instant  at  which  compression  waves  arrive  there. 

Typical  tine -correlated  pressure  at  the  tip,  light  intensity  at  the  propellant 
surface,  and  temperature  signals  from  the  heat  flux  gage  are  shown  in  Fig.  18. 

The  temperature  trace  shown  is  the  amplified  signal  through  a  wide  band  amplifier 
at  a  gain  setting  of  50  db.  The  ignition  event  is  also  shown  in  the  figure.  As 
mentioned  earlier,  the  time  for  ignition  is  first  obtained  from  the  film.  A  more 
accurate  location  of  ignition  is  obtained  from  the  photodiode  trace,  where  a  trough 
followed  by  an  uprising  of  the  curve  exists  in  the  neighborhood  of  the  point  obtained 
from  the  film.  In  this  study,  ignition  delay  is  defined  as  the  time  lag  between 
arrival  of  the  pressure  front  at  the  crack  tip  and  subsequent  ignition  (i.e.,  emission 
of  luminous  light  from  the  propellant  surface). 

Figures  19  and  20  show  two  streak  photographs  of  tests  conducted  at  low 
pressurization  rate  (17  GPa/s)  and  high  pressurization  rate  (55  GPa/s) ,  respectively. 
One  startling  difference  between  the  two  photographs  is  evident.  At  a  higher  pressuri 
zation  rate  (Fig.  20) ,  immediately  after  the  hot  gases  reach  the  crack  tip,  a  luminous 
region  appears  there.  As  time  progresses,  this  luminous  region  moves  downstream 
toward  the  crack  entrance.  Depending  upon  test  conditions,  the  luminous  region  may 
or  may  not  travel  all  the  way  to  the  crack  entrance.  This  bright  region  may  be 
caused  by  heat  released  from  combustion  behind  the  reflected  compression  waves  of 
unreacted  species  or  particles  from  the  igniter  system  which  are  carried  to  the  crack 
tip  by  strong  pressure  gradients.  At  higher  pressurization  rates,  the  pressure  gradi¬ 
ent  is  steeper,  and  the  igniter  species  have  less  time  to  burn  or  react  completely 
before  entering  the  combustion  chamber.  Both  of  these  p-.  >cesses  would  contribute  to 
the  appearance  of  the  bright  luminous  region  at  the  tip. 


Propellant  Sample: 
AP/PBAA-EPON 


Fig.  18  Time  Correlated  Pressure,  Light 

Intensity,  and  Thermocouple  Signals 
(Test  No.  DNI  3-42) 
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The  observed  luminous  (reaction)  region  could  also  be  caused  by  reaction  of 
igniter  species  with  oxygen  present  in  the  chamber.  Because  of  rapid  pressurization 
of  the  crack  cavity,  the  air  present  in  the  cavity  could  be  adiabatically  compressed 
causing  a  flash  reaction  between  the  oxygen  in  the  air  and  some  of  the  igniter  species 
To  eliminate  the  possibility  that  this  reaction  may  have  been  caused  only  by  the 
presence  of  oxygen,  tests  were  also  conducted  in  the  chamber  purged  with  nitrogen. 
Figure  21  shows  photographs  for  a  test  in  which  air  was  present  in  the  chamber, 
whereas  the  photograph  in  Fig.  22  is  for  a  test  in  which  the  chamber  was  purged  with 
nitrogen.  It  is  evident  in  both  cases  that  the  luminous  region  is  present,  even 
though  the  luminosity  is  somewhat  lower  when  oxygen  is  not  present.  Therefore, 
oxygen  present  in  the  chamber  enhances  luminosity  but  is  not  the  only  cause  of  it. 

Tests,  in  which  both  the  igniter  grain  and  a  sufficient  quantity  of  powdered 
propellant  were  placed  in  the  igniter  chamber,  were  conducted  to  further  investigate 
the  effect  of  igniter  species  or  particles  in  causing  the  luminous  zone.  Figure  23 
presents  a  set  of  pictures  for  this  case.  The  chamber  was  purged  with  nitrogen. 

A  comparison  of  Figs.  22  and  23  shows  that  the  luminosity  of  this  zone  is  higher 
for  the  case  shown  in  Fig.  23.  If  more  particles  that  can  be  entrained  in  the 
flow  are  present,  the  probability  that  the  particles  may  be  carried  to  the  tip 
with  the  gas  flow  and  burn  there  increases.  This  is  especially  true  at  higher 
pressurization  rates  because  the  particles  may  not  have  sufficient  time  to  react 
completely  in  the  igniter  chamber  or  the  main  chamber. 

Figure  24  shows  smoothed  temperature-time  traces  measured  by  using  the  thin 
film  heat-flux  gage  for  three  different  pressurization  rates.  It  can  be  seen  from 
the  figure  that  the  rate  of  increase  of  surface  temperatue  increases  as  the  pressuri¬ 
zation  rate  is  increased.  Deduced  heat  flux  as  a  function  of  time  for  the  conditions 
corresponding  to  Fig.  24  is  given  in  Fig.  25.  The  shapes  of  the  heat  flux-tine 
traces  are  dictated  by  their  respective  temperature-time  plots.  It  should  be 
noted  that  the  heat  flux  to  the  propellant  surface  is  extremely  high  under  these 
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operating  conditions  (in  the  order  of  10  MW/m  or  250  cal/cm"-s) ,  and  that  the 
heat  flux  to  the  propellant  surface  increases  with  the  increased  rate  of  pressuri¬ 
zation.  For  the  time  period  specified  in  Figs.  24  and  25,  measured  temperature 
continually  increases  and  heat  flux  to  the  propellant  surface  is  always  positive 
because  the  hot  igniter  gases  continue  to  flow  into  the  crack  cavity  during 
this  short  period  of  time. 

The  increase  in  heat  flux  to  the  propellant  surface  at  an  increasing  rate 
of  pressurization  can  be  attributed  to  a  combination  of  the  following.  At  higher 
pressurization,  heat  transfer  to  the  propellant,  from  both  the  reflected 
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compression  wave  and  the  hot  recirculating  gases  near  the  tip,  will  be  higher. 
Photographs  of  the  ignition  event  presented  in  Figs.  20-23  show  that  at  higher 
pressurization  some  of  the  unreacted  species  or  particles  from  the  igniter  grain 
or  booster  propellant,  which  did  not  have  sufficient  time  to  burn  completely, 
may  become  entrained  and  carried  by  the  hot  gases  to  the  tip  and  burn  there. 

Because  of  continued  pressurization  of  the  crack  cavity,  the  enthalpy  of  gases 
adjacent  to  the  propellant  surface  will  increase,  which  in  turn  increases  the 
rate  of  heat  transfer  to  the  propellant  surface. 

The  effect  of  pressurization  rate  at  the  tip  on  measured  ignition  delay 
time  is  shown  in  Fig.  26.  It  can  be  seen  that  the  ignition  delay  decreases  as 
the  pressurization  rate  increases  because  the  rate  of  heat  transfer  to  the 
propellant  surface  increases  as  the  pressurization  rate  is  increased.  It  is 
interesting  to  note  that  the  ignition  delay  is  in  the  submillisecond  range  for 
pressurization  rates  of  the  order  of  10  GPa/s  or  higher.  Scatter  in  the  experi¬ 
mental  data  appears  to  be  more  pronounced  at  lower  pressurization  rates;  this 
may  be  attributed  to  the  lower  rate  of  heat  transfer  to  the  propellant  under 
these  conditions.  Figure  26  also  shows  that  crack-gap  width  has  no  distinguishable 
effect  on  the  ignition  process.  Data  for  the  smallest  gap  width  (6  =  1.17  mm) 
appear  to  be  concentrated  at  lower  pressurization  rates;  because  no  diaphragm  was 
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inserted  between  the  igniter  and  combustion  chambers  for  tests  with  the  smallest  gap 
width,  lower  pressurization  rates  were  obtained. 

Experiments  were  also  conducted  to  evaluate  the  effect  of  igniter  gas 
temperature  on  the  ignition  process.  An  aluminized  AP-based  composite  propellant 
was  used  as  the  igniter  for  these  studies.  The  flame  temperature  of  this  propellant 
is  3600  K.  Because  of  the  very  high  luminosity  of  the  product  gases  from  this  pro¬ 
pellant,  it  was  difficult  to  accurately  measure  the  ignition  delay.  Limited  results 
obtained  from  these  tests  indicate  that  the  ignition-delay  time  is  somewhat  lower 
for  high- flame  temperature  igniter  gases.  This  is  to  be  expected,  since  the  rate  of 
heat  transfer  to  the  propellant  will  increase  with  higher  gas  temperatures. 

Table  3  lists  the  values  and  sources  of  physical  properties  and  input  parameters 
used  in  the  calculations.  Thermal  properties  of  AP  were  taken  from  Ref.  26.  Thermal 
properties  of  PBAA  are  approximate,  and  are  identical  to  those  used  by  Varney  and 
Strahle.16  Density  of  the  fuel  binder  in  the  table  is  weighted  average  density  of 
PBAA  and  EPON.  Conductivity  of  PBAA  was  deduced  from  this  average  density  and  the 
value  of  thermal  diffusivity  used  in  Ref.  16 for  fuel  binders.  The  values  of  in-depth 
radiation  absorption  coefficient  B  are  considered  to  be  equal  for  both  oxidizer  and 
fuel,  and  are  identical  to  those  used  in  Refs.  29-31.  Separate  values  of  $  for  AP 
and  PBAA  were  not  found  in  propellant  literature.  In  the  present  study,  radiation 
is  less  than  2.5%  of  the  total  heat  flux  to  the  propellant;  in-depth  radiation 
absorption,  therefore,  does  not  play  a  significant  role  in  the  overall  ignition  process. 

The  ignition  criterion  for  the  theoretical  analysis  was  based  upon  attainment 

of  a  critical  temperature  at  the  oxidizer  fuel-binder  interface  on  the  propellant 

surface.  This  temperature  was  denoted  as  623  K  (350°C) ,  which  is  the  accepted 

14  34 

decomposition  temperature  of  AP.  ’  Because  of  the  heterogeneous  nature  of  the 
propellant,  the  AP  and  the  fuel  binder  have  different  temperatures  at  any  axial 
location;  therefore,  the  ignition  criterion  was  based  upon  the  temperature  at  the 
AP-fuel  interface.  Photochemical  reactions  were  not  considered  in  the  computations. 
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Table  3  Physical  Properties  and  Input  Parameters 


Property 

Units 

AP 

PBAA 

c 

J/kg-K 

1305. 426, 27 

1255. 216 

V 

kJ/gmol 

13428 

14216 

k 

W/m-K 

0.4626 

0.1216 

L 

pm 

50 

— 

^py 

kj/kg 

-50210’1S 

81S916 

R1 

pm 

30.62 

-- 

R2 

pm 

— 

39.76 

zpy 

s'1 

1.0  x  109  (Ref.  28) 

2.0  x  107  (Ref.  16) 

8 

m'1 

1.0  x  104  (Ref.  29-31) 

1.0  x  104  (Refs.  29- 

P 

kg/m3 

19S032 

94733 
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Figure  27  shows  the  comparison  between  measured  and  predicted  ignition  delays. 
Theoretical  values  appear  to  predict  accurately  the  functional  dependence  of 
ignition  delay  on  the  pressurization  rate.  The  variation  in  the  predicted  ignition 
delay  when  t  .<sre  is  a  +25  K  change  in  the  ignition  temperature  is  also  shown  in 
the  figure.  At  low  pressurization  rate,  the  variation  is  large  because  of  lower 
heat  flux  to  the  propellant.  It  should  be  noted  that  only  representative  values  of 
predicted  ignition  delays  are  shown  in  Fig.  27;  there  is  also  some  scatter  in  the 
predicted  results  because  actual  measured  heat  flux  is  an  input  to  the  theoretical 
model . 

Figure  28  presents  a  plot  of  the  calculated  radial  variation  of  propellant 
surface  temperature  at  three  different  times.  The  temperature  profiles  at  these 
specified  times  are  quite  similar.  Temperature  in  the  fuel  binder  is  higher  than 
that  in  the  oxidizer  crystal  because  the  thermal  diffiisivity  of  fuel  is  lower. 

Even  though  the  pyrolysis  reaction  for  the  fuel  is  endothermic,  it  has  little 
effect  on  the  temperature  profiles  since  appreciable  pyrolysis  reaction  does  not 
occur  in  the  temperature  range  shown  in  Fig.  28.  It  is  also  evident  that  except 
near  the  oxidizer-fuel  interface,  the  radial  temperature  gradients  are  very  small 
as  a  result  of  adiabatic  radial  boundaries. 

Figure  29  shows  the  predicted  axial  variation  of  the  oxidizer  fuel  interface 
temperature  at  various  times  for  conditions  corresponding  to  Fig.  28.  The  thennal 
wave  penetration  depth  is  extremely  small  (<  20  um)  because  the  heat  flux  to  the 
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propellant  surface  is  very  large  (in  the  order  of  10  MW/m”)  and  the  time  duration 
is  very  short  (submilliseconds).  If  a  significant  portion  of  the  incident  heat 
flux  is  radiative,  the  thermal  wave  will  penetrate  deeper  than  that  shown  in  the 
figure.  In  the  present  study,  as  noted  earlier,  a  very  small  portion  (<2.5%)  of 
the  incident  heat  flux  is  radiative. 

Figure  30  is  a  plot  of  the  calculated  propellant  surface  temperature  at  the 
oxidizer  fuel  interface  as  a  function  of  time  for  pressurization  rates  of  15.35, 


Fig.  28  Calculated  Radial  Variation  of  Propellant 
Surface  Temperature 
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Fig.  30  Calculated  Variations  of  Propellant  Surface 
Temperature  at  the  Oxidizer  Fuel  Interface 
with  Tims 


59 


49.72,  and  112  GPa/s.  The  input  heat  flux  for  these  cases  is  presented  in  Fig.  25. 
Since  the  heat  flux  to  the  propellant  surface  increases  with  higher  pressurization, 
the  average  slopes  of  the  curves  in  Fig.  30  also  increase  with  higher  pressurization. 
The  time  variation  of  surface  temperature  fjr  i-.  given  pressurization  rate  depends 
strongly  on  the  corresponding  heat-flux  time  data. 
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V.  SUMMARY  OF  PROGRESS  ANC  CONCLUSIONS 

Progress  made  during  the  past  year  and  the  results  obtained  can  be  summarized 
as  follows: 

1.  Detailed  chemical  kinetics  information  for  the  AP/PBAA  propellants 
was  incorporated  into  the  proposed  theoretical  model  for  ignition  of  composite 
solid  propellants.  This  enables  a  more  complete  description  of  source  terms  in  the 
gas-phase  energy  and  species  equations. 

2.  Solid-phase  energy  equations  were  coded  and  implemented  on  a  computer. 

The  numerical  scheme  employs  a  stable,  implicit  scheme  for  solution  of  the  governing 
equations.  Numerical  solutions  of  the  solid-phase  subprogram  were  successfully 
checked  against  analytic  solutions  for  some  limiting  cases.  Using  actual,  measured 
heat  flux  near  the  propellant  surface  as  an  input  to  the  solid-phase  equations, 
predictions  for  the  ignition  delay  time  can  be  made. 

3.  Numerical  solution  of  the  complete  model  is  in  progress.  An  overall  solution 
scheme  was  established,  and  a  portion  of  the  gas-phase  finite  difference  equations 
was  coded. 

4.  The  test  chamber  was  modified  to  house  the  heat  flux  gage,  pressure  trans¬ 
ducer,  etc.,  near  the  propellant  sample.  A  new  igniter  system,  which  uses  electric 
primers,  was  designed  and  fabricated.  With  the  present  setup,  pressurization  rates 
in  the  o^der  of  100  GPa/s  (~106  atm/s)  were  obtained. 

5.  The  data  acquisition  system  was  updated  by  acquiring  and  installing  a 

2  MHz  transient  waveform  recorder  and  digitizer;  by  interfacing  the  digitizer  with 
a  minicomputer;  and  by  upgrading  the  movie  camera  system  to  reach  40,000  pictures  per 
second.  A  photodiode  system  was  installed  to  detect  onset  of  ignition. 

6.  Test  firings  were  conducted  to  investigate  effects  of  pressurization  rate 
dP/dt,  crack-gap  width  6,  and  igniter  flame  temperature  on  the  ignition  process 
The  ignition  delay  time  was  simultaneously  recorded,  using  the  Hycam  camera  and  the 
fast -response  photodiode  system. 
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7.  Heat  flux  to  the  propellant  surface  was  deduced  from  the  temperature-time 
profile  measured  by  a  fast-response  (5-10  ys)  thin  film  thermocouple,  using  a 
numerical  integration  procedure. 

8.  Experimental  results  indicate  that  the  ignition  delay  time  decreases  and 
the  heat  flux  to  the  propellant  surface  increases  as  the  pressurization  rate  is 
increased.  No  distinguishable  effect  of  the  crack-gap  width  on  the  ignition  process 
was  evident.  Limited  results  obtained  by  using  an  aluminized  propellant  as  the 
igniter  show  that  the  ignition  delay  time  is  somewhat  lower  for  high-flame  temperature 
igniter  gases. 

9.  Theoretical  predictions  made  using  the  solid-phase  equations  show  reasonable 
comparison  with  experimental  data.  Theoretical  results  also  indicate  that  the 
ignition  delay  time  decreases  as  the  pressurization  rate  is  increased. 

10.  The  decrease  in  ignition  delay  with  increasing  pressurization  is  caused  by 
enhanced  heat  feedback  to  the  propellant  surface  at  higher  pressurization  rates. 

This  augmentation  in  heat  feedback  to  the  propellant  at  higher  pressurization  is  a 
result  of  a  combination  of  the  following  mechanisms:  heating  due  to  conpresuion  wave 
reflection  at  the  closed  end;  heat  transfer  due  to  recirculating  hot  gas  near  the 
tip;  heat  release  by  burning  of  unreacted  igniter  species  (or  particles)  near  the 
tip,  behind  the  compression  wave;  and  increase  in  enthalpy  of  the  gas  adjacent  to 
the  propellant  surface  due  to  continued  pressurization  of  the  crack  cavity. 
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VI.  NOMENCLATURE 
Description 

axial  transformation  constant  for  the  gaseous  region 
specific  heat 

binary  mass  diffusion  coefficient 

activation  energy 

black  body  radiation  intensity 

external  radiation  heat  flux 

local  radiation  flux  at  z 

thermal  conductivity 

rate  constants  for  reaction  i 

thickness  of  tne  oxidizer  pellet 

constant  defined  in  Eq.  (31) 

pressure 

heat  flux  (.energy  per  unit  time  per  unit  area) 

heat  generation  rate  (energy  per  unit  time  per  unit  volume) 

rate  of  heat  generation  due  to  photochemical  process  per  unit 
wavelength/ local  radiation  flux  at  z 

heat  of  pyrolysis  per  unit  mass 

radial  distance  from  the  center  of  the  statistically  averaged 
element 

burning  rate 

universal  gas  constant 

radii  of  outer  surfaces  of  oxidizer  particle  and  fuel  binder, 
respectively 

transformed  nondimensional  axial  coordinate  in  gas  phase 
time 

temperature 

gas  phase  velocity  in  z  direction 
average  molecular  weight 

mass  fraction  of  species  j  it.  the  gas  phase 
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Symbol 


Subscripts 


1,2, 3, 4,5 


,  various  reactions  or  species 
fuel 

fuel -gas  interface 
gas  phase 
initial  value 
Species  j 
oxidizer 

oxidizer-fuel  interface 
oxidizer-gas  interface 
photochemical 
pyrolysis 
solid  phase 
solid-gas  interface 


in  z  direction  (perpendicular  to  the  initial  position 
the  interface) 

far  away  from  the  interface  in  z  direction 
wavelength 
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SUMMARY 


The  effect  of  interaction  between  combustion  processes  and  structural 
deformation  of  solid  propellant  was  considered.  The  combustion  analysis  was 
performed  on  the  basis  of  deformed  crack  geometry,  which  was  determined  from 
the  structural  analysis.  On  the  ocher  hand,  input  data  for  the  structural 
analysis,  such  as  pressure  distribution  along  the  crack  boundary  and  ablation 
velocity  of  the  crack,  were  determined  from  the  combustion  analysis.  The  inter¬ 
action  analysis  was  conducted  by  combining  two  computer  codes,  a  combustion 
analysis  code  and  a  general  purpose  finite  element  structural  analysis  code. 


INTRODUCTION 


In  recent  years,  much  attention  has  been  focused  on  the  investigation  of 
the  coupling  effect  between  combustion  phenomenon  and  mechanical  behavior  of  solid 
propellant.  The  solution  of  problems  of  this  type  can  further  better  understand¬ 
ing  of  the  transient  combustion  processes  inside  solid  propellant  cracks,  which 
may  significantly  affect  the  performance  of  a  rocket  motor.  The  combustion  pheno¬ 
menon  inside  the  crack  of  solid  propellant  is  strongly  influenced  by  the  crack 
geometry  as  the  material  is  being  deformed  and  burned  away.  Generally,  there  are 
two  major  reasons  for  alteration  of  the  crack  geometry:  1)  mass  loss  due  to 
gasification  of  propellant  surface  along  the  crack  during  the  comb  scion  process, 
and  2)  mechanical  deformacion  of  the  propellant  due  to  pressure. 

On  one  hand,  both  the  burning  rate  and  mechanical  deformation  are  governed 
by  pressure  acting  on  the  crack  surface.  On  the  ocher  hand,  a  change  in  crack 
sice  will  cause  che  pressure  distribution  to  vary.  The  pressure  distribution 
will  strongly  influence  che  deformation  and  stress  concentration  at  che  crack 
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Cip,  which  In  earn  will  affect  the  manner  of  the  crack,  propagation.  It  is, 
therefore,  apparent  that  the  pressure  distribution  and  the  change  in  crack  geo¬ 
metry  are  strongly  Interdependent. 

In  the  past,  combustion  and  structural  analyses  of  solid  propellant  were 
conducted  independently,  with  the  result  that  Interaction  effects  were  completely 
ignored.  As  noted  above,  such  Interaction  effects  can  be  quite  Important,  especi¬ 
ally  when  the  deformation  is  large  as  compared  to  the  original  crack-gap  width. 

The  deformation  response  of  the  material  is  categorized  as  linearly  viscoelastic. 

It  is,  therefore,  the  Intent  of  this  paper  to  present  a  method  of  analysis  for  the 
combustion-structural  interaction  in  a  linear  viscoelastic  medium.  To  this  end, 
three  major  tasks  are  Involved:  1)  combustion  analysis  to  model  the  transient 
combustion  process,  2)  viscoelastic  analysis  in  conjunction  with  moving  boundary, 
and  3)  linkage  of  the  two  analyses. 

For  the  combustion  analysis,  investigations  of  certain  aspects  of  combustion 
processes  have  been  made.  Taylor  [1]  conducted  experimental  tests  to  study  the 
convective  burning  of  porous  propellants  with  closed-  and  open-end  boundary  condi¬ 
tions.  Belyaev  et  al.  [2]  showed  that  the  burning  of  propellant  Inside  a  narrow 
pore  may  lead  to  an  excess  pressure  buildup.  In  a  later  study,  Belyaev  at  al.  [3] 
made  a  series  of  experimental  tefits  to  determine  the  dependence  of  flame-spreading 
rate  on  crack  geometry,  propellant  properties,  boundary  conditions,  and  combustion 
chamber  pressures.  Cherepanov  [4]  stated  that  as  a  result  of  the  impeded  gas  flow 
in  a  sufficiently  narrow  and  long  cavity,  the  pressure  reaches  such  high  values 
that  the  system  becomes  unstable.  From  his  work,  Godai  [5]  indicated  Chat  there  is 
a  threshold  diameter  or  critical  width  of  a  uniform  cavity  below  which  flame  will 
not  propagate  into  the  crack.  Krasnov  et  al.  [b]  investigated  the  rate  of  pene¬ 
tration  of  combustion  into  the  pores  of  an  explosive  charge.  Jacobs  et  al.  [7,8] 
studied  the  pressure  distribution  in  burning  cracks  that  simulate  the  debonding 
of  solid  propellant  from  the  motor  casing. 

Although  results  of  previous  experiments  were  of  interest,  no  sound  theore¬ 
tical  model  was  developed.  In  this  study,  a  theoretical  model  was  established 
for  predicting  rate  of  flame  propagation,  pressure  distribution,  and  pressuriza¬ 
tion  rata  inside  the  crack.  Two  secs  of  coupled  partial  differential  equations 
were  obtained:  one  from  mass,  momentum,  and  energy  conservation  of  the  gas  phase 
of  the  propellant  product  in  the  void  region  adjacent  to  the  crack  surface;  the 
other  from  consideration  of  solid-phase  heat  conduction.  Due  to  the  mathematical 
complexity  of  governing  equations  and  boundary  conditions  involved,  the  finite  ■ 
difference  method  was  used  to  obtain  the  solution  for  the  combustion  analysis. 

In  the  numerical  solution,  the  boundary  conditions,  which  vary  with  time,  are 
specified  in  terms  of  the  changing  crack  geometry,  which  in  turn  is  found  from  the 
structural  analysis.  In  addition,  the  pressure  distribution  along  the  crack  sur¬ 
face,  varying  as  a  function  of  time,  was  obtained  from  the  analysis  and  was  used 
as  input  for  the  structural  analysis. 

For  structural  analysis,  different  approaches  have  been  taken  previously 
in  solving  (analytically  or  numerically)  several  moving  boundary  problems  in  linear 
viscoelasticity.  Lee  et  al.  [9]  obtained  a  solution  for  the  pressurization  of 
an  annihilating  viscoelastic  cylinder  contained  by  an  elastic  casing  in  which 
the  material  was  assumed  to  be  a  Kelvin  model  in  shear  and  incompressible  in  bulk. 
Arenz  et  al.  [ 10 ]  performed  a  similar  analysis  for  a  sphere.  Corneliusser.  et 
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al.  (11,12]  presen cad  solutions  for  a  spinning,  annihilating,  viscoelastic 
cylinder  with  free  outer  boundary.  Since  the  constraining  case  was  not  included 
in  their  analyses,  the  stress  distribution  was  independent  of  the  material  pro¬ 
perties.  With  the  assumption  of  a  standard  linear  solid  model,  Shinozuka  [13] 
presented  the  analytical  solution  for  a  case-bonded  pressurized  viscoelastic 
cylinder.  More  generalized  solutions  were  obtained  by  Rogers  at  al.  [ 14]  for  a 
class  of  linear  viscoelastic  problems  by  using  the  numerical  integration  scheme. 
Schapery  [IS]  also  developed  a  general  method  for  solving  moving  boundary  problems. 
In  his  approach,  the  moving  boundary  condition  was  replaced  by  a  fictitious  non- 
toovlng  boundary  subjected  to  a  time-dependent  pressure.  Later,  Christensen  et 
al.  [16]  obtained  a  series  solution  for  the  stresses  of  the  same  problem.  As 
noted  above,  most  of  the  analytical  solutions  were  available  for  viscoelastic  prob¬ 
lems  of  simple  geometry.  For  complex  geometry,  the  finite  element  method  has 
proven  to  be  most  useful. 

Application  of  the  finite  element  method  for  solving  viscoelastic  problems 
is  not  new;  reports  of  such  work  can  be  found,  for  instance,  in  references  [17-20]. 
However,  most  of  the  previous  work  did  not  consider  the  effect  of  moving  boundary, 
an  important  feature  for  the  structural  analysis  of  solid  propellant.  Sankaran 
and  Jana  [21]  presented  a  technique  for  the  solving  of  axisymmetrlc  viscoelastic 
solids  with  moving  boundary.  In  their  approach,  the  finite  element  mesh  corres¬ 
ponding  to  the  new  boundary  was  re-generated,  while  the  stress-strain  histories 
and  material  properties  were  assumed  to  be  carried  over  from  chose  of  the  previous 
time  increment.  This  assumption  is  valid  only  if  the  time  increment  is  very  small. 
An  algorithm  for  automatically  cracking  ablating  boundaries  was  given  by  Weeks  and 
Cost  [22].  All  previous  work  dealing  with  moving  boundary  viscoelastic  problems 
lacks  both  the  appropriate  treatment  of  material  properties,  and  stress-strain 
histories  for  the  newly  generated  mesh.  It  is  Che  purpose  of  this  paper  to  pre¬ 
sent  such  a  treatment. 

Three  major  features  must  be  included  in  the  structural  analysis  for  a  solid 
propellant:  1)  proper  modeling  of  viscoelastic  behavior,  2)  tracking  of 
ablating  boundary  in  order  to  generate  new  finite  element  meshes,  and  3)  treat¬ 
ment  of  the  material  responses  (i.e.,  stress-strain  histories  and  material  pro¬ 
perties)  for  the  new  mesh.  All  of  these  features  have  been  incorporated  ince  a 
nonlinear  finite  element  program  called  NFAP  [23].  Combustion  and  structural 
programs  were  combined  in  order  to  make  possible  an  interaction  analysis.  Numeri¬ 
cal  results  are  presented  to  demonstrate  the  effect  of  interaction  beeveen 
combustion  and  structural  responses  of  the  material. 


COMBUSTION  ANALYSIS 


The  theoretical  model  was  developed  to  simulate  the  combustion  phenomenon 
inside  a  propellant  crack,  which  is  located  in  a  transverse  direction  to  the 
main  flow  of  the  rock  chamber.  During  the  course  of  derivation,  the  following 
assumptions  are  made; 

1)  All  chemical  reactions  occur  near  the  propellant  crack  surface,  and  the 
combustion  zone  is  so  thin  that  it  is  considered  a  plane. 

2)  Rate  processes  at  the  propellant  surface  are  quasi-steady  in  the  sense  that 
characteristic  times  associated  with  the  gaseous  flame  and  preheated  pro¬ 
pellant  are  short  ir  comparison  to  that  of  pressure  transient  variation. 
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3)  Gases  flowing  in  Chs  propallant  crack  obey  the  Clausius  or  Nobla-Abel 
aquation  of  staca. 

4)  Bulk  flow  in  tha  pore  is  one-dimensional  [24]. 

To  describe  gas-phase  behavior  inside  a  solid  propellant  crack,  mass, 
momentum,  and  energy  equations  in  unsteady,  quasi-one-dlmensional  forms  have 
bean  developed,  based  upon  the  balance  of  fluxes  in  a  control  volume  within  the 
propellant  crack. 

The  mass  conservation  equation  is 
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The  momentum  conservation  equation  is 
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The  energy  conservation  equation  written  in  terms  of  the  total  stored  energy 
(internal  and  kinetic)  per  unit  mass,  E,  is 
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The  conservation  equations  are  further  simplified  by  an  order  of  magnitude 
analysis  in  which  the  following  terms  are  negligible:  1)  forces  between  mole¬ 
cules  due  to  viscous  normal  stress  in  axial  direction;  2)  viscous  dissipation 
and  rate  of  work  done  by  the  force  caused  by  viscous  normal  stresses  in  the  energy 
equation;  and  3)  axial  heat  conduction  batween  gas  molecules  in  the  energy 
equat ion. 


The  propellant  surface  temperature  at  a  fixed  location  along  the  crack  before 
the  attainment  of  ignition  is  calculated  from  the  solid-phase  heat  conduction 
equation  written  in  unsteady  one-dimans ional  form: 


(4) 
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where  the  length  variable  y  Is  measured  perpendicular  to  the  local  propellant 
crack  surfaca.  Initial  and  boundary  conditions  ara 
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Tha  heat  conduction  aquation  Is  solved  by  using  an  Integral  method  [23] 
which  employs  a  third-order  polynomial,  or  by  direct  numerical  solution  of  Eqs. 
(4-7)  with  variable  mesh  size  in  the  subsurface. 


For  the  gas  phase,  the  Noble-Abel  equation  Is  used  for  the  equation  of 
state: 


P(i  -  b)  -  RT  (8) 

The  gas-phase  equations,  l.e.  Eqs.  (1),  (2)  and  (3),  ere  non-linear.  Inhomo¬ 
geneous,  partial  differential  equations.  Along  with  the  partial  differential 
equation  for  the  solid  phase  (Eq.  (4)),  they  are  solved  simultaneously,  using  the 
finite  difference  method.  The  derivation  described  above  was  implemented  into  a 
computer  program,  crack  combustion  code  (CCC)  by  Kuo  at  al.  [26]. 

STRUCTURAL  ANALYSIS 


To  conduct  the  structural  analysis  of  the  solid  propellant,  three  main 
features  must  be  included  in  the  numerical  formulations:  1)  modeling  of 
viscoelastic  material  behavior,  2)  simulation  cf  ablating  boundary,  and 
3)  treatment  of  material  responses  by  an  interpolation  scheme.  Each  feature 
is  outlined  below. 


Viscoelastic  Material  Model 

The  material  behavior  of  the  solid  propellant  is  assumed  to  be  visco- 
elastlc  in  shear  and  elastic  in  bulk.  Only  the  isothermal  condition  is  considered 
The  stress-strain  relations  with  zero  initial  conditions  ars  written  in  two  parts. 


1)  Shear  behavior: 
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where  is  Che  relaxation  modulus  in  shear.  For  most  viscoelastic 
materials,  it  is  usually  considered 


2)  Bulk  behavior: 
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As  discussed  in  (27],  the  incremental  stress-strain  relations  in  matrix  form 
are  written  as 
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and  the  term  C^,  has  a  recursive  relationship,  i.e., 
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The  advantage  of  Eq.  (19)  is  chac  all  of  the  strain  history  can  be  obtained  by 
referring  only  to  information  in  the  previous  time  step,  thus  reducing  computer 
storage  and  numerical  calculations. 

From  the  virtual  work  principle  and  the  relationship  of  Eq.  (12),  the  finite 
element  equilibrium  equations  for  a  typical  time  interval  [t,  t  +  At]  can  be 
derived  as 
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where 
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Simulation  of  Ablating  Boundary  and  Mesh  Generation 

Burning  of  the  propellant  causes  a  significant  change  in  geometry,  thus 
presenting  complications  in  finite  element  structural  analysis.  The  effect  of 
ablating  boundary  is  accounted  for  by  redefining  the  finite  element  mesh  at  speci¬ 
fied  time  intervals.  This  involves  two  stages  of  calculations:  1)  tracking 
of  the  ablating  boundary,  and  2)  generation  of  new  finite  element  mesh.  With 
some  modifications,  the  procedures  adopted  herein  are  similar  to  those  presented 
in  [22]. 

Consider  a  structural  geometry  with  ablating  boundary.  The  spatial  posi-  . 
tions  of  the  ablating  boundary  are  determined  by  the  ablation  velocities  which 
are  found  from  the  combustion  analysis  at  discrete  times.  It  is  assumed  that 
the  ablation  occurs  always  in  the  directon  normal  to  the  boundary.  For  struc¬ 
tural  analysis,  the  entire  surface  is  divided  into  an  ablating  part  and  a  non¬ 
ablating  part;  each  part  is  formed  by  discrete  line  segments  joining  at  the  nodes 
of  the  finite  element  mesh.  The  new  position  of  each  line  segment  is  located  from 
the  given  ablating  velocity.  Consequently,  the  new  boundary  nodes  are  determined 
by  calculating  the  intersections  of  two  subsequent  new  line  segments.  Likewise, 
the  nodes  at  the  intersections  of  new  ablating  and  non-ablating  boundaries  are 
then  determined. 

During  the  locating  process,  however,  some  of  the  boundary  nodes  may  not  lie 
on  the  new  boundary  and  thus  must  be  eliminated.  If  the  distance  from  the  tip 
of  the  normal  vector  at  a  new  nodal  position  to  any  node  on  the  original  boundary 
is  less  than  the  value  of  the  normal  itself,  the  node  is  removed. 


75 


In  general,  the  total  number  of  nodes  on  the  boundary  at  discrete  times 
will  be  different  because  some  of  the  nodes  have  been  removed.  However,  in 
the  analysis  it  is  more  convenient  to  generate  a  finite  element  mesh  similar  to 
Che  original  one  so  that  interpolation  of  material  response  can  be  made.  One 
way  to  accomplish  this  is  by  keeping  the  number  of  boundary  nodes  constant. 
Consequently,  the  boundary  nodes  are  redistributed  between  two  discontinuity 
points  which  are  specified  in  the  input  data  in  such  a  way  that  the  lengths 
of  the  new  line  segments  have  the  same  ratio  as  those  of  the  original  lines. 

Once  the  new  boundary  nodes  are  defined,  an  automatic  mesh  generation 
scheme  is  used  to  create  the  interior  nodes  for  further  analysis.  Because  of 
its  flexibility  in  obtaining  a  desirable  mesh,  a  Laplacian-isoparametrlc  grid 
generation  scheme  [28]  is  utilized.  However,  this  method  is  limited  to  a 
geometry  bounded  by  four  sides.  A  finite  element  mesh  is  shown  in  Fig.  1.  The 
coordinates  of  the  i-th  interior  node  can  be  expressed  in  terms  of  those  of 
neighboring  nodes  by 


*i 

zi 


1 

4(2-w) 

1 

4(2-w) 


(2(yu  +  y12  *  y13  *  yi4> 

-  “''IS  +  y16  +  y17  +  y18)l 

(24) 

[2<2tl  +  2 12  +  2 13  +  V 

-  V(215  +  216  *  217  +  218>I 

(25) 

where  w  is  the  weighting  factor  for  adjusting  the  distribution  of  Interior  nodes, 
and  0  £  w  <_  1. 

Setting  up  the  equations  for  each  interior  node  yields  two  systems  of  simul¬ 
taneous  equations.  It  is  observed  that  the  resulting  systems  of  equations  are 
banded  and  symmetric.  The  Gaussian  elimination  scheme  is  employed  to  solve  for 
Che  coordinates  of  the  interior  nodes. 

Interpolation  of  Material  Responses 

As  seen  from  Eq.  (16),  the  stress  increment  Ao  for  the  time  interval 
[t,  t+At ]  varies  with  material  properties  and  with  the  strain  history  at  both 
current  and  previous  time  steps.  When  the  region  of  an  element  changes  over  a 
period  of  time  due  to  ablation,  the  material  response  history  of  the  new  elements 
is  lost  and  must  be  determined  by  an  interpolation  procedure  from  Che  old  ele¬ 
ments  at  previous  time  steps.  Accordingly,  the  Interpolation  procedure  is  carried 
out  on  the  element  level.  For  calculations,  the  material  responses  are  separated 
into  cm  groups:  the  first  includes  such  variables  evaluated  at  the  Gaussian 
integration  points,  i.e.,  Aa  ,  Ac  and  Cj  ;  the  second  Includes  the  nodal 
displacements  which  are  evaluated  at  noda?  points.  In  the  present  calculations, 
two  limitations  are  Imposed:  1)  eight -node  quadrilateral  elements  are  used 
throughout  the  analysis;  and  2)  the  four  sides  of  each  element  remain  straight 
before  and  after  ablation. 

1)  Interpolation  of  Gaussian  variables  -  It  is  noted  that  the  quadratic  dis¬ 
placement  approximation  of  an  eight-node  element  yields  a  linear  strain  varia¬ 
tion.  With  this  fact  in  mind,  the  quantities  of  Gaussian  variables  at  -odal 
points  are  first  evaluated  for  every  old  element.  As  shown  in  Fig.  2a,  b,  this 
can  be  done  by  using  the  linear  isoparametric  shape  functions,  namely, 
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where  -  1/4(1  +  r  r)(l  +  SjS)  i  »  1,  2,  3,  4 

For  each  new  element,  the  local  coordinates  (r  ,  s)  of  the  k-th  Gaussian 
point  are  known.  The  global  coordinates,  (y^,  z^).  of  that  point  are.  therefore, 
computed  by  using  the  following  equations: 
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After  (y  ,  z  )  are  found,  the  old  element  to  which  the  point  belongs  must  be 
identified.  A  search  process  based  upon  the  values  of  r’  and  s'  is  developed  for 
this  purpose.  The  search  9tart3  from  the  old  element  which  corresponds  to  the 
neighboring  elements.  Equations  for  such  calculations  are  given  by 
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Fig.  2c  shows  how  to  identify  the  element  to  which  the  points,  (r* ,  s’),  belong. 
Once  the  location  of  the  point  is  verified,  an  Interpolation  procedure  is  per¬ 
formed,  using  the  relationship 
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2)  Interpolation  of  nodal  displacements  -  A  similar  procedure  to  that  explained 
above  f-  also  used  to  determine  the  position  of  the  node  in  question  with  refer¬ 
ence  to  the  old  element.  However,  the  interpolation  procedure  in  Eq.  (26)  Is  no 
longer  necessary  since  Che  nodal  displacements  are  known.  The  nodal  displacements 
of  the  new  mesh  are  computed  from 
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where  h  are  the  standard  quadratic  isoparametric  shape  functions. 


(30) 


All  formulations  discussed  in  this  section  have  been  implemented  into  a 
general  purpose  nonlinear  finite  element  program  called  NFAP  for  conducting 
viscoelastic  analysis  of  solid  propellant  with  ablating  boundary.  Some  numer¬ 
ical  examples  are  presented  in  a  later  section. 


COUPLING  EFFECT 


For  structural  analysis,  the  boundary  condition  along  the  crack  geometry 
is  defined  by  pressure  distribution  which  varies  with  time,  and  ablation  velocity; 
both  are  determined  from  combustion  analysis.  In  the  combustion  analysis,  the 
regression  rate  of  the  propellant  is  dependent  on  the  deformed  crack  geometry. 
Therefore,  the  two  processes  are  strongly  interdependent.  Such  a  coupling  effect 
is  obtained  by  combining  the  analysis  of  two  computer  programs:  a  crack  combus¬ 
tion  code  (CCC)  and  a  structural  analysis  code  (NFAP).  Both  codes  were  developed 
independently  to  facilitate  program  verifications.  Linkage  of  the  two  codes  was 
made  subsequently. 

The  coupling  effect  considered  in  the  present  analysis  is  limited  to  the 
major  parameters,  namely  pressure  loading,  ablation  velocity,  and  crack  deforma¬ 
tion.  Pressure  and  ablation  velocity  are  calculated  by  the  CCC  at  each  nodal  point 
located  on  a  or*- dimensional  grid  along  the  length  of  the  crack.  The  analysis 
of  crack  combustion  incorporates  the  crack  geometry  variation  caused  by  both 
mechanical  deformation  and  mass  loss  through  gasification  of  the  propellant 
surface.  Once  the  gas-phase  equations  are  solved  and  the  pressures  and  ablation 
velocities  along  the  crack  are  calculated  for  a  given  time  t,  the  data  are  trans¬ 
ferred  to  the  NFAP  as  the  input  information.  NFAP  then  simulates  the  updated  crack 
geometry  from  the  ablation  velocities  and  generates  a  new  finite  element  mesh. 

With  the  new  mesh  and  pressure  data,  NFAP  updates  the  stiffness  matrix  and  inter¬ 
polates  material  responses  for  conducting  a  quasi-static  analysis  at  time  t.  After 
obtaining  the  deformation,  the  change  in  the  crack  width  at  each  finite  different 
node  is  calculated  and  added  to  the  existing  crack  width.  Since  the  crack  width  i3 
the  input  of  the  combustion  analysis,  one  cycle  of  calculations  is  thus  completed. 
The  same  procedure  is  followed  for  every  specified  time  increment. 


EXAMPLES 


For  program  verification  and  demonstration  of  its  analysis  capability,  three 
sample  problems  were  run  either  by  NFAP  alone  or  in  the  combined  NFAP/CCC  program. 
The  results  of  the  analysis  are  discussed  in  the  following.  The  numerical  results 
obtained  from  CCC  alone  are  contained  in  reference  [26], 

1.  A  Reinforced  Thick-walled  Cylinder 

Figure  3  shows  a  cylinder  of  viscoelastic  material  bonded  by  a  steel  casing 
and  subjected  to  a  step-function  internal  pressure.  The  example  was  selected 
because  it  is  competed  of  two  different  materials  and  the  analytical  results 
are  readily  available  for  comparison.  Only  five  eight-node  axisymmetric  elements 
were  used  to  model  the  cylinder.  The  material  properties  of  the  elastic  casing 
are 

E  -  2.068  x  106  MPa 


v  -  0.3015 
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The  material  properties  of  the  viscoelastic  core  are  defined  by 


K  -  689.5  MPa  C  -  51.71  *  exp  (-0.lt)  MPa 

In  Fig.  4,  the  variations  of  circumferential  stresses  with  time  are  plotted 
for  comparison  with  the  analytical  solution  obtained  in  [9].  It  is  observed  that 
both  solutions  agree  very  closely.  This  problem  was  analyzed  previously  by 
Zienkiewicz  at  al.  [18],  using  strain  rate  formulation  of  the  finite  element 
method.  However,  the  formulation  presented  in  the  present  paper  is  more  easily 
incorporated  into  the  NFAP  program. 

2.  A  Star-shaped  Solid  Rocket  Motor 

As  an  application  of  the  present  approach  in  dealing  with  the  moving  boundary, 
a  star-shaped  solid  rockec  motor  was  analyzed  by  assuming  both  a  constant  and 
ablating  inner  boundary.  The  configuration  and  finite  element  me3h  are  shown 
in  Fig.  5,  and  the  material  properties  of  outer  casing  and  inner  propellant  are 
identical  to  those  of  the  first  example.  Taking  advantage  of  the  symmetry 
condition,  only  a  30°-sector  was  modeled  by  finite  element  mesh.  The  contours 
of  maximum  compressive  stress  analyzed  by  constant  inner  boundary  at  various 
times  are  shown  in  Fig.  5.  Comparing  the  present  results  with  those  of  [18],  it 
is  evident  that  the  general  pattern  is  quite  similar  but  that  some  small  differences 
do  exist.  Since  the  geometry  of  the  rocket  motor  in  [18]  was  not  clearly  defined, 
the  difference  in  dimension  used  in  these  two  analyses  could  be  the  cause  of  such 
deviations. 

The  actual  case  of  a  solid  rocket  motor  can  be  modeled  more  closely  by  con¬ 
sidering  the  inner  boundary  being  ablated.  Figure  6  shows  the  contours  of  maximum 
compressive  stress  predicted  by  NFAP,  using  the  option  of  moving  boundary.  The 
results  obtained  are  quite  different  from  chose  of  [18].  However,  observing  the 
differences  between  Figs.  5  and  6,  we  can  conclude  that  the  results  obtained  by 
NFAP  are  quite  reasonable.  The  solutiou  reveals  that  the  high  stress  region 
obtained  for  ablating  boundary  propagates  faster  than  chat  with  non-ablating 
boundary . 

3.  A  Propellant  Crack  Specimen 

As  a  final  example,  a  propellant  crack  sample  was  analyzed,  using  the 
combined  NFAP/CCC  program  to  demonstrate  the  coupling  effect.  The  initial  geo¬ 
metry  and  finite  element  mesh  generated  by  NFAP  is  given  in  Fig.  7.  The  crack 
is  0.15  m  long  and  the  Initial  gap-width  is  0.89  am.  The  web  chlckness  is 
8  an  along  the  crack  and  20  mm  at  the  tip.  Because  of  symmetry,  only  half  of  the 
sample  was  modeled  by  80  plane  scrain  elements.  The  shear  relaxation  modulus  of 
the  propellant  was  assumed  to  be 

(t)  -  1.461  +  7.43  *  exp  (-.095c)  MPa;  and  K  -  4,826  MPa. 

Calculated  pressure  distributions  at  various  times,  from  the  CCC  alone,  are 
given  in  Fig.  7.  The  burning  phenomenon  of  the  propellant  can  be  briefly  des¬ 
cribed  as  follows.  The  pressure  in  Che  chamber  increases  with  time,  causing 
the  hot  gases  to  penetrate  further  into  the  crack.  As  time  passes,  the  pressure 
wave  travels  along  the  crack  and  is  reflected  from  the  closed  end.  Ac  about 
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200  us,  Che  pressure  front  has  already  reached  the  tip  and  is  reflected,  causing 
pressure  at  the  tip  to  increase. 

Figure  8  shows  Che  results  obtained  from  the  combined  NFAP/CCC  program. 

During  the  initial  period,  the  general  trend  of  the  pressure  distribution  is  simi¬ 
lar  to  that  from  convective  burning  analysis  alone.  However,  as  time  progresses, 
noticeable  differences  between  the  two  cases  begin  to  appear.  Up  to  200  u»*  the 
pressures  obtained  from  the  combined  analysis  are  lower,  except  near  the  crack 
entrance  region.  At  t  ■  300  Us,  two  pressure  peaks  appear.  At  t  -  325  ps,  three 
pressure  peaks  appear.  These  pressure  peaks  are  caused  by  the  partial  closure  of 
the  gap.  The  deformation  pattern  of  the  propellant  is  quite  Irregular  because 
of  the  uneven  distribution  of  the  pressure  along  the  crack  surface.  The  elements 
at  the  crack  entrance  are  compressed  by  the  high  chamber  pressure,  which  results 
in  the  propellant  being  pushed  into  the  crack.  Since  chamber  pressure  Increases 
more  quickly  than  pressure  Inside  the  crack,  the  propellant  is  pushed  toward  the 
lower  pressure  region  inside  the  crack.  The  mechanical  deformation  of  the 
propellant  causes  narrowing  of  the  crack  width,  and  consequently  results  in  a  local 
crack  closure.  This  local  gap  closure  manifests  itself  in  a  pressure  peak.  The 
localized  pressure  peaks  or  gap  closures  move  along  the  crack.  At  t  *  325  us,  this 
localized  pressure  phenomenon  becomes  evident  at  x/L  «  0.167,  0.433,  and  0.633. 


CONCLUSION 


The  computer  program  for  evaluating  the  coupling  effect  between  convective 
burning  and  structural  deformation  was  developed  by  combining  the  Crack  Combustion 
Code  and  a  Nonlinear  Finite-Element  Analysis  Program.  In  structural  analysis, 
the  linear  viscoelastic  material  model,  together  with  the  capabilities  of  simu¬ 
lating  ablating  boundary  and  interpolating  material  responses,  was  considered. 
Also,  the  coupling  effect  estimated  by  the  combined  analysis  shows  some  signifi¬ 
cant  interaction  between  the  combustion  and  mechanical  deformation.  This  pheno¬ 
menon  will  be  verified  further  by  future  experiments. 
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SYMBOLS 

1.  Combustion  Analysis 

cross-sectional  area  of  crack 
body  force 

A 

b  *  co-volume 

0^  -  specific  heat  at  constant  pressure 
C  -  total  stored  energy 

h  »  local  convective  heat-transfer  coefficient 

c 

local  convective  heat-transfer  coefficient  over  propellant  surface 
local  convective  heat-transfer  coefficient  over  nonpropellant  port  wall 
enthalpy  of  combustion  gas  at  adiabatic  flame  temperature 
burning  perimeter 
wetted  perimeter  of  port 
static  pressure 

specific  gas  constant  for  combustion  gases 

burning  rate  of  solid  propellant,  Including  erosive  burning  contribution 
temperature  (without  subscript,  static  gas  temperature) 
adiabatic  flame  temperature  of  solid  propellant 
initial  propellant  temperature 
propellant  surface  temperature 
nonpropellant  wall  surface  temperature 
t  -  time 
u  *  gas  velocity 

V  velocity  of  propellant  gas  at  burning  surface 
axial  distance  from  propellant  crack  opening 
perpe  .dicular  distance  from  propellant  surface  into  solid 
a  ■  thermal  diffusivlty 
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Y  ■  ratio  of  specific  heats 
X  *  thermal  conductivity 
U  “  gas  viscosity 

p  ■  density  (without  subscript,  gas  density) 

t  ■  shear  stress  on  port  wall 
w 

T  ■  normal  viscous  stress 

XX 

©w  •  angle  measure,  in  a  counterclockwise  direction,  at  lower  side  of 
propellant,  degree 

Subscripts 

i  *  initial  value 

pr  a  solid  propellant  (condensed  phase) 
c  ■  rocket  chamber 
2.  Structural  Analysis 

■  stress  deviators 
e^  ■  strain  deviators 

■  stress  tensor 

>  shear  relaxation  modulus 
K  ■  bulk  modulus 
(  )C  ■  a  quantity  at  time  t 
{Ao}  ■  incremental  stress 
{Ae}  ■  incremental  strain 

{oq}  ■  equivalent  initial  stress  vector  due  to  viscoelastic  behavior 
[K ]  *  stiffness  matrix 

[  ■  transpose  of  matrix 

M  -  number  of  terms  of  series  in  relaxation  modulus 


material  constants  in  relaxation  modulus 


-  increment  of  nodal  displacement  vector 

■  viscoelastic  material  matrix 

■  strain-nodal  displacement  transformation  matrix 

■  values  of  Gaussian  variables  at  k-th  Integration  point  referred 
to  old  element 

*  values  of  Gaussian  variables  at  i-th  nodal  point  referred  to 
old  element 

■  local  coordinates  of  k-th  Integration  point  referred  to  old 
element 

*  global  coordinates  of  i-th  nodal  point  referred  to  new  element 

*  i-th  nodal  displacement  referred  to  new  element 

*  1-th  nodal  displacement  referred  to  old  element 

*  local  coordinates  of  point  in  equation  referred  to  old  element 
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Figure  3.-  Finite  element  mesh  of  a  reinforced  thick-walled  cylinder 
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Finite  element  mesh  of  a  propellant  crac*  and 
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Figure  8.-  Calculated  pressure  distributions  for  various 
times  from  the  combined  crack  combustion  and 
non-linear  finite-element  analysis  program. 
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APPENDIX  B  (Ref.  No.  4) 

EFFECT  OF  PROPELLANT  DEFORMATION  ON  IGNITION 
AND  COMBUSTION  PROCESSES  IN  SOLID  PROPELLANT  CRACKS* 

M.  Kumar ^  and  K.  K.  Kno^ 

Department  of  Mechanical  Engineering 
The  Pennsylvania  State  University 
University  Park,  Pennsylvania  16802 

ABSTRACT 

A  comprehensive  theoretical  model  was  formulated  to  study  the  develop¬ 
ment  of  convective  burning  in  a  solid  propellant  crack  which  continually 
deforms  due  to  burning  and  pressure  loading.  In  the  theoretical  model,  the 
effect  of  interrelated  structural  deformation  and  combustion  processes  was 
taken  into  account  by  considering  a)  transient,  one-dimensional  mass, 
momentum,  and  energy  conservation  equations  in  the  gas  phase,  b)  a  trans¬ 
ient,  one-dimensional  heat  conduction  equation  in  the  solid  phase,  and 
c)  quasi-static  deformation  of  the  two-dimensional,  linear  viscoelastic 
propellant  crack  caused  by  pressure  loading.  This  set  of  coupled,  nonlinear, 
partial  differential  equations  was  solved  numerically.  Several  regions  of 
partial  crack  closures  were  observed  experimentally  in  narrow  cracks 
(“450  ym) .  Predicted  results  indicate  that  the  partial  closures  may  generate 
substantial  local  pressure  peaks  along  the  crack,  Implying  a  strong  coupling 
between  chamber  pressurization,  crack  combustion,  and  propellant  deformation, 
especially  when  the  cracks  are  narrow  and  the  chamber  pressurization  rates 
are  high.  Predicted  results  for  ignition  front  propagation  and  pressure 
distribution  are  in  good  agreement  with  the  experimental  data.  Both  theore¬ 
tical  and  experimental  results  indicate  that  the  maximum  pressure  in  the 
crack  cavity  is  generally  higher  than  that  in  the  chamber.  It  was  found 
that,  under  the  conditions  studied, the  initial  flame-spreading  process  is 
not  substantially  affected  by  propellant  deformation. 

INTRODUCTION 

High-energy  propellants,  which  are  used  to  obtain  increased  specific 
impulse  in  rocket  motors, generally  contain  high  solids  loading  of  energetic 
materials,  e.g.,  cvclotetramethylenetetranitramine  (HMX) .  As  the  density 
of  solids  in  the  propellant  is  increased,  the  probability  that  cracks  and 
flaws  will  develop  in  the  propellant  grain  is  also  increased.  Defects  in 
propellant  grains  can  originate  during  manufacture,  storage,  or  handling, 
during  ignition  and  combustion,  or  during  loading  tf  different  segments  of 
a  segmented  rocket. 

The  extent  to  which  cracks  can  reduce  the  reliability  of  solid  rocket 
motors  has  been  a  major  concern  in  the  development  of  high-energy  propellant 
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grains.  It  is  known  that  cracks  in  solid  propellant  grain  can  allow  hot, 
high-pressura  gases  to  penetrate  the  cavity,  thereby  providing  additional 
surface  area  for  combustion.  The  phenomenon  of  such  a  rapid  flame  propaga¬ 
tion  into  a  propellant  cavity  and  the  subsequent  rapid  regression  is  called 
convective  burning.  The  convective  burning  rate  often  exceeds  the  normal 
(conductive)  burning  rate.  Mechanical  deformation  and  crack  propagation 
may  result  when  combustion  processes  inside  the  crack  produce  much  higher 
pressure  than  the  designed  maximum  pressure.  If  the  local  pressure  rise  due 
to  the  gasification  is  sufficiently  rapid,  it  may  produce  strong  compression 
waves,  or  even  shock  waves,  which  can  initiate  detonation.  The  generally 
understood  mechanism  of  transition  from  deflagration  to  detonation  (DDT)  is 
given  in  Ref.  1.  DDT  is  more  likely  to  occur  in  gas-permeable  propellants, 
such  as  granular  propellant  beds  or  propellant  grains  with  flaws  or  cracks, 
which  have  large  burning  surface  areas.  Even  if  the  convective  burning  rate 
does  lead  to  catastrophic  failure,  the  combustion  inside  the  propellant 
crack  can  cause  the  performance  of  a  rocket  motor  to  deviate  significantly 
from  expectations. 

Convective  burning  Inside  a  crack  is  complicated  because  it  involves 
several  interdependent  processes  such  as  a)  pressure  wave  phenomenon, 
b)  convective  heating  of  the  propellant  crack  surface  and  subsequent  flame 
propagation  along  the  crack,  c)  pressurization  of  the  crack  cavity  due  to 
burning,  which  may  result  in  flow  reversal  when  the  pressure  inside  the 
crack  exceeds  that  of  the  chamber,  d)  change  in  the  crack  geometry  due  to 
burning,  as  well  as  mechanical  (^formation  caused  by  pressure  loading,  and 
e)  propagation  of  the  crack  due  to  regression  and  mechanical  fracture. 
Parameters  which  may  significantly  affect  the  combustion  process  are:  cham¬ 
ber  pressure  and  pressurization  rate,  geometry  of  the  crack,  erosive  burning 
effects  due  to  high  gas  velocities,  physicochemical  properties  of  the  pro¬ 
pellant,  composition  of  the  igniter  gas,  and  initial  and  boundary  conditions. 

This  paper  deals  with  the  formulation  of  a  comprehensive  theoretical 
model,  which  Includes  the  effect  of  propellant  deformation  to  predict  flame 
spreading  and  combustion  processes  during  the  development  of  convective 
burning  in  solid  propellant  cracks.  Predicted  results  are  also  compared 
with  some  of  the  experimental  investigations2 >3  conducted  by  the  authors  for 
model  validation. 

As  mentioned  earlier,  convective  flame  propagation  and  subsequent  crack 
propagation/branching  may  lead  to  anomalous  burning  or  even  DDT.  Since  the 
propellant  is  deformable,  burning  inside  a  solid  propellant  crack  is  basic¬ 
ally  a  coupled  solid  mechanics  and  combustion  phenomenon.  The  combustion 
phenomenon  inside  the  crack  is  strongly  influenced  by  the  crack  geometry. 

The  geometry,  or  gap-width,  can  be  altered  for  the  following  two  reasons: 

1.  Mass  loss  due  to  gasification  or  ignition  of  the  propellant 
surface  along  the  crack  during  the  combustion  process,  leading 
to  a  variation  in  the  geometry. 

2.  Mechanical  deformation  of  the  propellant  crack  due  to  a)  pres¬ 
sure  loading  on  the  Interior  surface  of  the  crack,  b)  pressure 
force  of  the  combustion  chamber  acting  on  the  propellant  grain 
near  the  crack  entrance,  and  c)  stresses  in  the  propellant. 
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caused  by  temperature  gradients  in  the  grain,  which  may  also 
cnuse  changes  in  Che  crack  geometry. 

Both  burning  rate  and  the  mechanical  deformation  are  governed  by  pres¬ 
sure  acting  on  the  crack  surfaces.  On  the  other  hand,  a  change  in  gap  width 
or  geometry  will  cause  pressure  distribution  to  vary.  Pressure  distribution 
will  strongly  influence  deformations  along  the  crack  and  stress  concentra- 
ion  at  the  crack  tip.  It  is  obvious,  therefore,  that  pressure  and  crack 
geometry  are  strongly  interdependent  in  a  burning  solid  propellant  crack. 

Convective  burning  inside  a  solid  propellant  crack  has  been  e  subject 
of  Interest  in  recent  years.  An  extensive  literature  review  In  this  area 
was  recently  conducted  by  Bradley  and  Boggs. ^  Belyaev  et  al.  have  compiled 
most  of  the  recent  Russian  work  oh  convective  burning  and  DDT.  Both  theore¬ 
tical  and  experimental  studies  in  the  area  of  crack  combustion  have  been 
conducted  at  The  Pennsylvania  State  University.®”*®  Convective  burning 
studies  can  be  subdivided  into  two  broad  categories:  a)  onset  of  convective 
burning,  and  b)  development  of  convective  burning.  This  study  is  restricted 
to  the  development  of  convective  burning  in  isolated  propellant  cracks.  The 
following  paragraphs  list  some  of  the  important  research  conducted  in  this 
area. 


Experimental  and  theoretical  studies  on  the  development  of  convective 
burning  are  limited.  Belyaev  et  al.^>^  made  preliminary  experimental  ^ 
Investigations  of  the  development  of  combustion  in  single  pores.  Kuo  et  al. 
developed  a  detailed  theoretical  modal  for  predicting  the  development  of  con¬ 
vective  burning  in  isolated  propellant  cracks.  However,  this  model6  does 
not  include  the  mechanical  deformation  of  the  crack  during  the  development 
of  convective  burning.  Kiml3  studied  the  possibility  of  shock  to  detonation 
transition  in  propellant  cracks  using  the  model  of  Kuc  et  al.®  and  an 
extremely  simplified,  one-dimensional,  elastic  model  for  propellant  deforma¬ 
tion.  Pilcher1^  has  also  reported  some  computations  conducted  on  a  deform¬ 
able  crack,  but  has  given  no  details. 

High  pressures  present  in  the  crack  during  tha  development  of  convec¬ 
tive  burning,  as  well  as  associated  mechanical  deformation  and  atresa  concen¬ 
trations,  may  lead  to  crack  propagation.  Most  of  the  earlier  work  in  tha 
area  of  crack  propagation  was  dona  In  tha  USSR.*'"*'  Kirsanova  and  Leipun- 
skii*5  examined  the  mechanical  stability  of  propellant  cracks,  using  steady- 
state  approximations.  Only  continuity  and  momentum  aquations  vara  considered 
in  tha  gas  phase.  Tha  erosive  burning  affects  wars  ignored,  lha  material 
was  assumed  to  be  elasto-brlttle,  and  stress  intensity  was  used  to  determine 
tha  Initiation  of  crack  propagation.  Tha  stablo  crack  length  obtained  from 
their  simplified  modal  was  a  function  of  tha  crack  length,  chamber  pressure, 
burning  rata,  gas  enthalpy,  atraaa  intensity  factor,  and  propellant  density. 
Tha  eutostabillsation  condition  for  tha  unstable  burning  crack  waa  alto 
obtainad . 

A  greatly  Improved  analytical  modal  for  tha  aama  problem  waa  formula tad 
by  Cherepanov ,  Id  Ha  stated  that  If  tha  cavity  waa  sufficiently  long  and 
narrow,  tha  praasurea  la  tha  crack  would  roach  a  high  valua  and  render  tha 
ayatoa  unstable.  Even  though  hia  gas-phase  equations  contained  nonstoady 
terms,  several  important  terms  warm  ignored.  Tha  propallant  material  was 


assured  to  be  elnsto-plastic ,  and  the  deformation  process  to  be  quasistatic. 
Howevers  Cherepanov  did  not  solve  the  governing  set  of  equations  and,  there¬ 
fore,  did  not  present  any  results.  He  obtained  a  sufficient  condition  for 
the  stability  of  a  steady-state  problem. 

Belyaev  et  al.*^  proposed  modifications  to  the  Kirsanova  and  Leipun- 
skii*'  modal  and  pointed  out  some  of  its  limitations.  It  was  noted  that  an 
elasto-brlttle  material  model  cannot  describe  the  visco-elastic  character, 
which  is  typical  of  solid  propellants.  They  made  maximum  pressure  measure¬ 
ments  in  the  crack,  and  accordingly  changed  the  empirical  correlation  for 
maximum  pressure  In  the  cavity.  The  maximum  pressure  in  the  crack  was  a 
function  of  chamber  pressure  and  the  ratio  of  crack  length  to  half-width. 

They  proposed  that,  depending  upon  the  chamber  pressure  and  crack  geometry, 
the  stability  of  the  crack  can  be  subdivided  into  four  regions:  a)  a  region 
of  absolute  mechanical  stability;  b)  a  region  of  practical  stability;  c)  a 
region  of  hounded  growth;  and  d)  a  region  of  self-sustained  growth. 

18-20 

Jacobs  et  al.  studied  the  critical  nature  of  cracks  and  debonds, 

assuming  a  quasi-steady  model.  The  propellant  was  considered  to  be  linear 
viscoelastic.  They  also  measured  the  pressure  distribution  along  the  crack 
surface  under  the  condition  of  instantaneous  ignition  over  the  entire  sur¬ 
face.  The  ga?  dynamics  equations  of  Jacobs  et  al.  are  inadequate  because 
they  include  numerous  simplifications  and  assumptions.  Other  related  studies 
on  crack  combustion  have  been  reported . 21-23  jakata  and  Wiedermann22,23 
have  been  investigating  the  initiation  of  detonation  in  a  propellant  crack 
by  introducing  the  effect  of  stress  wave  interaction  with  the  crack  boundary, 
and  assuming  a  foam  layer  to  be  present  on  the  propellant  surface  during 
combustion.  The  propellant  is  considered  to  be  linear.  Isotropic,  elastic 
material  which  does  not  yield  or  fail. 

2A-3Q 

Numerouj  studies,  based  purely  on  solid  mechanics  considerations, 
have  been  conducted  on  crack  propagation.  All  are  based  on  quasi-static 
crack  propagation  in  viscoelastic  material  and  use  the  concept  of  local 
energy  dissipation  at  the  crack  tip.  Work  on  dynamic  crack  propagation  is 
currently  underway,  with  Swanson^O  reporting  some  progress  in  this  area. 
However,  there  is  at  present  no  theoretical  model  capable  of  realistically 
predicting  stress  concentrations  at  the  tip  of  a  burning  crack,  which  is  a 
necessary  input  for  all  crack  propagation  theories. 

It  is  apparent  that  a  theoretical  model  is  needed  to  predict  the  devel¬ 
opment  of  convective  burning,  while  also  taking  into  account  the  effect  of 
mechanical  deformation.  Such  a  model  would  also  be  extremely  useful  in 
making  realistic  predictions  of  crack  propagation;  at  present,  these  predic¬ 
tions  are  based  purely  on  solid  mechanics  considerations.  Even  though  this 
investigation  will  not  attempt  to  study  crack  propagation  caused  by  mechani¬ 
cal  fracture  at  the  tip,  the  analysis  and  program  developed  here  can  be 
extended  easily  for  such  a  study. 

The  specific  objectives  of  this  study  are: 

1.  To  develop  a  theoretical  model  to  study  ignition,  flame  spread¬ 
ing  and  convective  burning  in  a  solid  propellant  crack  which 
continually  deforms  due  to  burning  and  pressure  loading. 
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2.  To  study  th*  effaces  of  such  parameters  as  chamber  pressuriza¬ 
tion  rate,  dP/dt,  and  crack  gap  width  on  the  flame  spreading 
and  combustion  process  in  propellant  cracks. 

3.  To  teat  the  validity  of  the  theoretical  model  by  comparing  the 
predicted  results  with  the  experimental  data. 

Actual  propellant  defects  may  consist  of  several  Irregularly  branched 
cracks.  Howevar,  because  of  the  complexity  of  physical  processes  Involved 
in  the  study  of  combustion  in  propellant  cracks,  and  for  mathematical  tract- 
ability,  this  Investigation  Is  limited  to  a  single  isolated  crack. 

ANALYSIS 

DESCRIPTION  OF  PHYSICAL  MODEL 

A  schematic  diagram  of  the  physical  model  chosen  to  simulate  a  propel¬ 
lant  crack  is  shown  in  Fig.  1.  The  location  of  the  crack  is  normal  to  the 
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main  flow  direction  in  the  combustion  chamber  or  in  a  rocket  motor.  The 
current  analysis  is  applicable  to  a  crack  of  variable  geometry.  However,  in 
order  to  closely  simulate  the  assumptions  of  the  physical  model,  and  for  ease 
of  manufacturing  a  reproducible  propellant  crack,  specimens  used  in  the 
experimental  study  were  of  uniform  rectangular  cross-section.  A  schematic 
diagram  of  the  propellant  crack  geometry  used  in  the  experimental  investiga¬ 
tion  is  shown  iu  Fig.  2. 

As  the  hot  combustion  gases  flow  through  the  main  chamber,  a  portion  of 
the  gases  is  driven  into  the  crack  because  of  the  pressure  gradient.  And  as 
the  hot  gases  flow  over  the  propellant  crack,  energy  is  transferred  from  the 
gas  to  the  propellant.  The  rate  of  heat  transfer  to  the  propellant  depr  .o* 
upon  the  local  temperature,  velocity,  and  density  of  the  gas,  and  the  temp¬ 
erature  of  the  propellant  surface.  The  heat  transfer  may  cause  the 
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propeller'.  Co  ignite.  As  additional  gaaas  are  driven  into  Che  crack,  Che 
density  and  velocity  Increase  further,  resulting  in  a  faster  flame  propaga¬ 
tion  along  the  crack.  As  the  ignition  front  moves  downstream,  more  hot  gases 
ar'!  generated  by  combustion  of  the  propellant  surface,  causing  local  pressure 
to  increase.  At  the  same  time,  the  deformation  of  the  propellant  will  also 
alter  the  pressure.  The  net  result  of  these  two  opposing  and  Interdependent 
phenomena  determines  the  local  pressure. 

BASIC  ASSUMPTIONS 


The  following  basic  assumptions  are  made  in  the  derivation  of  the  theo¬ 
retical  model: 

1.  The  gas-phase  reaction  zone  is  considered  quasi-steady.  In 
other  words,  the  relaxation  time  associated  with  the  flame 
is  much  shorter  than  that  associated  with  transient  pressure 
variation;  therefore,  the  flame  adjusts  Itself  immediately 
to  chamber  conditions. 

2.  The  deformation  of  the  propellant  is  quasi-steady,  that  is, 
the  mechanical  deformation  of  the  propellant  can  be  obtained 
by  using  a  static  analysis. 

3.  All  chemical  reactions  occur  near  the  propellant  surface  in  a 
very  thin  planar  zone,  and  the  distance  of  the  combustion  zone 
from  the  propellant  surface  is  small  when  compared  to  the  crack 
gap  width. 

4.  The  bulk  flow  of  the  gases  in  the  pore  is  considered  to  be  one- 
dimensional. 

5.  The  gases  present  in  the  propellant  crack  obey  the  Noble-Abel 
gas  law. 


6.  The  propellant  surface  temperature  is  uniform,  and  the  thermal 
properties  of  the  propellant  are  constant. 

7.  The  propellant  can  be  represented  by  a  linear  viscoelastic 
material  model;  the  material  is  assumed  to  be  isothermal, 
homogeneous  and  isotropic.  The  propellant  grain  is  two- 
dimensional,  with  no  mechanical  fracture  at  the  crack  tip. 


The  validity  of  assumptions  1  and  2  can  be  shown  by  comparing  the  char¬ 
acteristic  times  associated  with  pressure  variation,  gaseous  flame  tone,  and 
compression  wave  propagation  through  the  propellant  grain.  In  the  crack 
combustion  experiments,  the  typical  time  associated  with  the  pressure  excur¬ 
sion,  T  ,  is  in  the  order  of  1  ms.  The  other  characteristic  times  are  evalu¬ 
ated  aspfollows: 


T„  <»/v  .2 

-8.  -  _ 8*  „ 

T  a  /r. 1 
s  pr  b 


X  p  c 

pr  pr  p 


where  is  the  characteristic  time  associated  with  the  solid  phase,  and  is 
given  by 

t  -  a  /r  2  •  0  (1  ms)  (2) 

s  pro 

Since  (Xc  ) / (X  c  )  «  0(1),  and  (p/p  )  »  0(.01),  we  get  T  *  0(.01  ms), 
pr  pr  p  pr  g 


The  time  associated  with  the  stress  waves  in  the  solid  is 


_ web  thickness _  m  0(.Q1  m)  _  0(.Q1  m) 

stress  wave  propagation  speed  o(/E7p)  0(2000  m/s) 


«  0(5  ys)  (3) 

Since  T  and  xy  are  both  <<T  ,  the  quasl-steedy  assumption  of  the  flame 
cone  and  mechanical  deformation  «?•  valid.  The  value  of  tv  calculated  from 
the  shear  modulus  will  be  an  order  of  magnitude  higher;  however,  it  will 
still  be  smaller  than  t  .  Assumption  3  is  reasonable  for  AP-based  composite 
solid  propellants  in  the  pressure  range  of  interest.  Even  though  the  heat 
release  sone  is  distributed,  the  diffusion  flame  stand-off  distance  is 
usually  in  the  order  of  50  ym,  and  is  small  in  comparison  to  the  gap  width 
(gap  widths  are  in  the  order  of  500  ym  in  the  present  experimental  investi¬ 
gation)  . 

Aa sumption  4  is  well  Justified  for  the  crack  dimensions  used  in  this 
investigation.  The  crack  length  to  hydraulic  diameter  ratio  was  approxi¬ 
mately  100.  The  dense  gas  relation  (assumption  5)  can  adequately  describe 
the  deperture  from  the  ideal  gas  law  at  high  pressures.  Assumption  6, 
employed  for  mathematical  simplicity,  allows  the  uce  of  a  one-dimensional 
transient  heat  conduction  equation  for  the  solid.  Since  the  thermal  wave 
penetration  depth  is  of  the  same  order  of  magnitude  as  the  oxidizer  particle 
size  of  the  propellant,  the  actual  heat  transfer  process  is  three-dimensional. 
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However,  researchers  in  the  field  have  not  yet  investigated  the  three- 
dimensional  treatment  of  the  heat  conduction  process  in  composite  solid  pro¬ 
pellants  because  of  the  complexities  associated  with  the  three-dimensional 
structure  of  oxidizer  crystals  and  fuel  bind  .r.  and  the  fact  that  numerical 
solutions  are  cumbersome  and  time  consuming. 


Most  of  the  solid  propellants  can  be  well  characterized  by  a  visco¬ 
elastic  material  model. 24-40  For  simplicity,  the  viscoelastic  material 
model  was  assumed  to  be  linear.  Even  though  the  assumption  of  isothermal 
condition  for  material  response  is  not  strictly  accurate,  it  is  a  good 
approximation,  3ince  the  thermal  wave  penetration  depth  is  less  than  100  ym. 
Fracture  of  the  propellant  at  the  crack  tip  was  considered  beyond  the  scope 
of  the  present  investigation,  since  it  is  a  vast  subject  area  in  itself. 

The  viscoelastic  nature  of  the  material,  i.e.,  the  time  dependence  of  the 
relaxation  modulus,  may  not  be  very  significant  during  the  period  of  experi¬ 
mental  test  firings  of  this  study  because  the  stress  relaxation  time  is  much 
greater  than  the  transient  pressure  variation  time.  Indeed,  the  analysis 
and  model  developed  here  are  applicable  to  a  wide  variety  of  operating  con¬ 
ditions  and  are  not  simply  limited  to  the  experimental  test  conditions  of 
this  study. 

In  this  analysis,  ignition  is  defined  as  the  attainment  of  a  critical 
temperature  at  the  propellant  surface.  It  is  assumed  that  no  solid  or  gas- 
phase  reaction  takes  place  before  onset  of  ignition.  The  tacit  assumption 
is  made  here  that  there  is  little  time  between  the  attainment  of  a  critical 
surface  temperature  and  the  reaction  between  the  fuel  and  the  oxidizer 
species  to  cause  ignition. 

CONSERVATION  EQUATIONS 


To  formulate  the  theoretical  model,  we  first  write  the  unsteady,  vari¬ 
able  area,  one-dimensional  gas-phase  conservation  equations  for  reacting 
compressible  fluid  flow.  The  conservation  equations  are: 

Mass  Conservation: 


D(pA  )  r)(puA  ) 
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Momentum  Conservation: 


-  P  T 


T  cos  0  +  p  A  B  - 
w  W  w  n  x 
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Energy  Conservat ton: 


(|T>A  E)  +  t—  (PA  uE)  »  t—  —  r~  (A  Pu) 

-’t  p  ax  p  ax  dx  p 
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An  order  of  magnitude  analysis  was  made;  the  terms  crossed  with  arrows 
in  Eqs.  (5)  and  (6)  wore  found  to  be  much  smaller  than  the  others,  and  hence 
were  dropped . 

Equation  of  state  for  the  gas-phase: 

P(|  -  b)  -  RT  (7) 

Equation  for  stress  equilibrium  in  the  solid-phase  (neglecting  body 
force) : 


Constitutive  law  for  linear  viscoelastic  material: 


Vc) 


VM>  dV  Eij(T>  dT 
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Gi(t-T)  "dT  dT  +  V  ekk  (9) 


The  transient  heat  conduction  equation  for  the  solid-phase,  before 
ignition: 

3T  32T 

— El  -  a  - El 

3t  V  3y2 

The  initial  and  boundary  conditions  for  Eq.  (10)  are: 
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V  (12) 


{T(t)-T  (C) )  (13) 

ps 

After  rearranging  the  gas-phase  conservation  equations  and  making  some 
simplifications,  velocity-variation,  temperature-variation,  and  pressure- 
variation  equations  were  obtained.  These  equations  were  found  to  be  totally 
hyperbolic  in  nature. 41 

The  gas-phase  and  solid-phase  equations  were  solved  separately;  this 
is  justified  because  of  the  quasi-steady  assumption.  Details  of  numerical 
scheme  for  solving  the  gasdynamics  equations  and  the  solid  mechanics  equa¬ 
tions  are  given  in  a  later  section. 

INITIAL  AND  BOUNDARY  CONDITIONS 
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Initial  and  boundary  conditions  must  be  specified  in  order  to  complete 
the  theoretical  formulation.  The  three  initial  conditions  necessary  for  the 
solution  were  specified  as 


u(0,x)  -  ui 

(H) 

T(0,x)  -  T£ 

(15) 

P(0,x)  -  Pi 

(16) 

The  number  of  physical  boundary  conditions  that  can  be  specified 
depends  upon  the  flow  conditions  at  the  opening  of  the  crack.  When  the 
gases  f lev  into  the  crack  at  a  subsonic  speed,  the  boundary  conditions  are 

P(t,0)  -  Pc(t) 

(17) 

T(t,0)  -  Tc(t) 

(18) 

u(t,xL>  -  0 

(19) 

When  the  gases  flow  out  of  the  crack,  and 
one  boundary  condition  can  be  specified  at  the 
conditions  are  given  as 

the  outflow  is  subsonic,  only 
crack  opening.  The  boundary 

P(t,0)  -  Pc(t) 

(20) 

u(t,xL)  -  0 

(21) 

When  the  gas  flowing  cut  of  the  crack  is  supersonic,  nc  boundary  condition 
can  be  specified  at  the  crack  entrance.  In  such  a  case,  the  only  specified 
boundary  condition  is 

u(t,xL)  -  0 

(22) 
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For  turbulent  flow,  the  final  expression  for  h  is 

f  ,1  „  -0.6  .  0.8  „  0.2  _  -0.67  ,  .  .-0.1 

h  *  0.0J46  Fr  c  (Pu/R)  V,  T  ,  (xd, ) 

c  p  k  ar  h 


(24) 


Before  local  ignition,  h^  *  h  ■  li  :  after  local  ignition,  h  *  0.  The 
Prandtl  numer  in  Eq.  (247pis  cSYculaled  from  Svehla's  equation?^'** 


Pr 


1.77y  -  0.45 


(25) 


The  correlation  for  the  friction  coefficient  for  turbulent  flow  used  in 
this  study  is 


0.4491  (dh/x) 
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where  e^/d  is  the  relative  equivalent  sand  roughness.  This  expression  is 
actually  a  modified  form  of  the  well-knovm  Colebrook  formula^, 48  for  turbu¬ 
lent  flow  in  pipes  with  roughness.  Entrance  effects  are  taken  into  account 
by  modifying  the  friction  coefficient  by  a  power  function  of  the  distance-to- 
diameter  ratio  (the  original  Colebrook  expression  was  obtained  for  x/d  2.  20). 
Additionally,  Eq.  (26)  implicitly  accounts  for  the  temperature  dependence  of 
the  gas  density  and  viscosity,  each  evaluated  at  the  average  film  temperature, 
T  ,.  After  the  propellant  surface  begins  to  burn  locally,  the  value  of  the 
friction  coefficient  is  set  at  zero,  due  to  the  attenuation  of  wall  shear 
stress  caused  by  surface  blowing. 

The  burning-rate  expression  is  obtained  from  the  strand  burning  rate 
data  and  is  given  by 


rfe  -  aPn  (27) 

The  effect  of  erosive  burning  is  taken  into  account  bv  using  either  the  ero¬ 
sive-burning  augmentation  factor  obtained  by  Razdan  and  Kuo,  9,o0 

n  n 

rb  -  aPn  [1  +  K(P)  p(u)  U1  (28) 

or  the  Lenolr-Robillard  burning-rate  formula, ^ 

r,  -  aPn  +  K  h  exp  (-8  r,  p  /u  c)  (29) 
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NUMERICAL  SCHEME  FOR  GASP VNAM ICS  EQUATIONS 

An  implicit  finite  difference  scheme  is  used  to  solve  the  nonlinear 
coupled  partial  differential  equations.  A  central  difference  method  was 
used  for  the  derivatives.  The  predictor-corrector  method  was  used  to  handle 
the  non-linear  nature  of  the  governing  equations.  A  quasilinearization 
method  was  used  to  linearize  the  Inhomogeneous  terms  of  the  governing 
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equations.  The  sec  or  finite  difference  equacions  thus  obtained  is  solved 
simultaneously,  using  a  block  f'idiugcnal  matrix  inversion  method. 

Six  boundary  conditions  are  required  when  the  central  difference  method 
is  used  for  spacewise  derivatives.  This  is  the  case  because  the  difference 
equations  for  the  first-order  spatial  derivative  obtained  by  using  this 
method  correspond  essentially  to  the  iecond-crder  partial  differential  equa¬ 
tion.  The  extraneous  boundary  conditions  needed  for  the  solution  of  the 
finite  difference  equacions  are  derived  from  the  compatibility  relationships 
at  the  boundaries.  These  relationships  were  obtained  by  transforming  the 
governing  hyperbolic  equations  into  their  characteristic  form.  The  equations 
determining  the  extraneous  boundary  conditions  are  all  first-order,  ordinary 
differential  equacions.  They  are  simultaneously  integrated,  using  a  fourth- 
order  Runge-Kutta  integration  technique. 

The  extraneous  boundary  conditions  used  are  as  follows.  When  the  gas 
flows  into  the  crack,  in  addition  to  the  physical  boundary  conditions  (17) 
to  (19),  the  three  extraneous  boundary  conditions  required  are  velocity  at 
the  crack  entrance,  and  pressure  and  temperature  at  the  crack  tip.  When  the 
gas  flows  out  of  the  crack  subsonically ,  the  four  extraneous  boundary  condi¬ 
tions  necessary  are  pressure  and  velocity  at  the  crack  entrance,  and  pressure 
and  temperature  at  the  crack  tip.  Wien  the  gas  flows  out  of  the  crack  super¬ 
sonically,  the  five  extraneous  boundary  conditions  needed  are  pressure, 
temperature  and  velocity  at  the  crack  entrance,  and  pressure  and  temperature 
at  the  crack  tip. 

NUMERICAL  SCHEME  FOR  SOLID-PHASE  HEAT  EQUATION 

The  governing  partial  differential  equation  for  die  solid  phase,  Eq.  (10), 
was  solved  by  using  an  implicit  finite  difference  technique  with  a  variable 
mesh  system.  The  variable  grid  spacing  provides  finer  grid  spacing  near  the 
surface.  The  central-difference  method  was  used  to  approximate  both  time  and 
spatial  derivatives.  The  resulting  set  of  simultaneous  algebraic  equations 
was  solved  by  using  a  standard  tridiagonal  matrix  inversion  method. 


NUMERICAL  SCtiEME  FOR  STRUCTURAL  ANALYSIS 


Mechanical  deformation  of  the  propellant  was  obtained  by  using  a  general 
purpose  structural  analysis  program  which  is  called  the  Nonlinear  Finite- 
Element  Analysis  Program  (NFAP)  .  This  is  an  extended  version  of  the  Non¬ 
linear  Structural  Analysis  Program  (NONSAP) ^  ’  .  Since  NONSAP  is  well- 

documented  and  easily  accessible,  details  of  the  numerical  procedure  used  in 
NONSAP  will  not  be  repeated  here.  NFAP  was  modified  and  consolidated  to  suit 
the  needs  of  the  present  study,  and  for  computaitonal  storage  efficiency. 
Three  major  features  added  to  the  numerical  formulations  are:  1)  modeling  of 
the  viscoelastic  material  behavior,  2)  simulation  of  the  ablating  boundary, 


and  3)  treatment  of  the  material  response 
of  this  formulation  can  be  found  elsewhere. 


by  an  interpolation  scheme.  Details 
8,  56 


To  test  the  numerical  procedure  developed  in  this  section,  the  modified 
version  of  NFAP  was  used  to  compute  several  viscoelastic  problems  for  which 
exact  solutions  arc  known.  The  agreement  between  the  NFAP  solution  and  the 
exact  solution  is  excellent.  Details  of  these  comparisons  are  given  in  Ref. 
56. 
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The  extraneous  boundary  conditions  needed  to  solve  the  equations  are  obtained 
by  using  the  characteristic  equations. 

A  schematic  diagram  of  the  boundary  conditions  for  two-dimensional 
structural  analysis  is  shown  in  Fig.  3.  Because  of  the  symmetry  of  the  crack 
specimens  (see  Fig.  2),  only  half  of  the  crack  is  considered.  All  surfaces 
of  the  propellant  are  either  pressure-loaded  or  have  fixed  boundaries.  The 
propellant  surface  exposed  to  the  chamber  experiences  the  chamber  pressure, 

P  (t),  whereas  the  surface  along  the  crack  experiences  the  local  pressure 
along  the  crack,  P(t,x).  Due  to  the  symmetry  along  the  center  line  of  the 
crack,  the  propellant  boundary  at  the  tip  is  allowed  to  have  deformations 
only  along  the  axial  direction.  Along  the  surfaces  of  the  propellant  in  con¬ 
tact  with  the  brass  mold,  the  boundaries  are  assumed  to  be  fixed. 


Only  Axial  Diiplacmast 


EMPIRICAL  CORRELATIONS 


To  close  the  system  of  governing  equations,  empirical  correlations  are 
needed  for  the  heat  transfer  coefficient,  drag  coefficient,  and  burning-rate 
law.  The  units  of  the  variables  in  the  following  correlation  are  in  cgs 
system,  and  the  units  of  pressure  are  in  gf/cm2. 


The  correlation  for  the  local  convective  heat  transfer  coefficient,  h  , 
is  deduced  from  the  Dlttus-Boelter  correlation  for  flow  in  a  pipe.  The  c 
effect  of  the  entrance  region  is  taken  into  account  by  using  a  power  function 
of  length  to  hydraulic  diameter  ratio.^2**3  The  gas  properties  are  evaluated 
at  an  average  film  temperature,  T  to  take  into  account  the  variation  of 
the  physical  properties  of  the  gas< across  the  boundary  layer.  The  dependence 
of  viscosity  on  temperature  was  included,  using  Svehla's  expression^  for 


viscosity  of  air  at  high  temperature.  Bartz^*  has  sh< 
between  Svehla's  equation  and  the  published  NBS  data 
for  viscosity  is 

P  _  If  n.0.65 


XP 


good  comparisons 
Svehla's  expression 


“  V 


where 


Vk- 


8.699  x  10 


-7  M  0.5 
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NUMERICAL  PROCEDURE 


Of  che  several  interdependent  gasdynamic  and  solid  mechanics  processes 
that  may  influence  combustion  behavior  in  a  propellant  crack,  the  coupling 
between  the  pressure  along  the  crack  and  the  structural  deformation  was  con¬ 
sidered  to  be  the  most  important.  Other  coupled  processes,  such  as  thermal 
stresses  due  to  the  temperature  gradient  in  the  solid,  variation  in  burning 
rate  due  to  compression  of  the  propellant,  etc.,  were  considered  to  be  of 
secondary  importance  end  ware  ignored  in  the  present  analysis. 

The  general  layout  of  the  numerical  procedure  is  shown  in  Fig.  4.  Pres¬ 
sure  is  celculated  by  the  gasdynamic  portion  of  che  computer  program  at  each 
nodal  point  on  a  one-dimensional  grid  along  the  length  of  the  crack.  The 
convective  burning  analysis  of  the  crack  combustion  incorporates  the  crack 
geometry  variation  caused  by  both  mechanical  deformation  and  mass  loss 
through  gasification  of  the  propellant  surface.  Once  the  gas-r  se  equations 
are  solved,  and  pressures  and  burning  rates  along  the  crack  ar  calculated 
for  a  particular  time  t,  program  control  is  transferred  to  the  structural 
analysis  (NFAP)  portion  of  the  combined  program.  NFAP  reformulates  the  geom¬ 
etry  because  of  the  material  loss,  and  updates  the  stiffness  matrix  for  the 
new  time  step.  Surface  elements  of  the  finite-element  mesh  are  loaded  with 
the  pressure  obtained  through  the  gas-phase  equations.  The  propellant  de¬ 
formation  is  calculated  in  NFAP  using  a  static  analysis  at  time  t.  The 
general  stress-strain  equations  are  solved  using  a  plane-strain  analysis; 
this  is  congruous  to  the  experimental  test  configuration.  The  transient 
nature  of  the  pressure  loading  is  considered  by  using  a  static  analysis  at 
incremental  time  steps. 

Figure  5  is  a  diagram  of  the  finite  element  grids  used  for  the  structure 
analysis  and  the  finite-difference  nodes  used  for  the  gas-phase  solutions. 

The  configuration  is  compatible  with  the  geometry  considered  in  the  theoret¬ 
ical  model  as  well  as  with  that  used  lr.  the  experiments.  Two-dimensional, 
eight-node,  isoparametric,  quadrilateral  elements  were  used  to  model  a 
sample  crack.  Because  of  symmetry,  only  half  of  the  crack  needs  to  be  con¬ 
sidered.  Of  the  eight  nodes  in  an  element,  four  are  located  at  the  corners, 
and  one  at  the  midpoint  of  each  side.  Except  near  the  tip  region,  the  finite- 
difference  nodes  used  for  the  gas-phase  solution  are  in  a  one-to-one  corres¬ 
pondence  with  the  finite-element  nodes  which  form  the  propellant  surface  in 
the  two-dimensional  structural  analysis  mesh.  This  is  helpful  in  reducing 
the  computational  effort  in  transferring  the  values  from  the  finite-difference 
nodes  to  the  finite-element  grids,  and  vice  versa. 

RESULTS  AND  DISCUSSION 


EXPERIMENTAL  RESULTS 


Figure  6  shows  a  schematic  diagram  of  the  test  configuration  used  for 
detailed  observation  of  the  flame  front.  The  crack  is  formed  between  a  pro¬ 
pellant  slab  and  an  inert,  transparent  plexiglass  window;  che  crack  is 
perpendicular  to  the  direction  of  flow  in  the  main  chamber.  This  type  of 
configuration  provides  direct  (front  view)  observation  of  the  ignition  front 
propagation  and  burning  propellant  surface,  even  for  very  narrow  cracks. 

Moat  of  the  experimental  results  were  obtained  with  cracks  formed  by  cutting 
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a  slot  of  desired  width  in  a  propellant  slab  (see  Refs.  2  and  3).  Details 
of  the  experimental  setup  and  experimental  results  are  given  in  Refs.  2  and  3. 
The  entrance  of  the  crack  gap  is  rounded  to  facilitate  the  gas  penetration. 

Figure  7  presents  a  sequence  of  pictures  obtained  from  a  test  performed 
on  a  183  mm  long,  narrow  crack  (6  455um),  The  framing  rate  in  this  case 

was  36,700  pictures  per  second.  The  crack  configuration  is  the  same  as  that 
shown  in  Fig.  6.  It  was  observed  that  when  the  width  of  the  crack  gap  is 
very  small  (<  500  um) ,  partial  closure  of  the  crack  occurs  due  to  propellant 
deformation.  This  partial  closure  first  appears  near  the  entrance.  As  time 
progresses,  the  closure  region  moves  downstream  and  a  second  partial  closure 
region  develops  near  the  entrance.  This  process  continues;  at  times,  three 
or  four  such  partial  closures  are  observed  simultaneously  (for  initial  gap 
widths  of  the  order  of  450  um).  Usually,  at  the  moment  when  the  third  closure 
develops  at  the  entrance,  the  second  and  first  are  located,  respectively,  at 
approximately  20%  and  50%  of  the  crack  length  from  the  entrance.  Later,  as 
the  combustion  process  becomes  more  pronounced,  the  partial  closure  regions 
disappear  as  a  result  of  both  propellant  regression  and  higher  pressure  in 
the  cavity. 

The  entire  process  is  believed  to  be  the  result  of  deformation  caused  by 
high  pressure  acting  on  the  propellant  surface  exposed  to  the  chamber,  and  by 
the  complex  interaction  between  propellant  deformation  and  pressure  distribu¬ 
tion  in  the  crack  cavity.  Propellant  surfaces  exposed  to  the  chamber  are 
compressed  by  the  high  chamber  pressure,  which  results  in  the  propellant  being 
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pushed  tnco  the  crack.  Since  pressure  Increases  faster  In  the  chamber  than 
In  the  crack,  during  this  initial  pressure  transient,  the  propellant  is 
pushed  toward  the  region  of  lower  pressure  inside  the  crack.  The  mechanical 
deformation  of  the  propellant  causes  narrowing  of  the  crack  gap,  and  conse¬ 
quently  results  in  local  creek  closure. 

THEORETICAL  PREDICTIONS  AND  COMPARISONS 

The  important  physical  properties  of  the  two  propellants  studied  are 
given  in  Table  1.  Numerical  values  of  the  Input  variables  used  in  the  pre¬ 
dictive  program  are  listed  in  Table  II.  The  recorded  pressure-time  trace 
near  the  crack  entrance  was  used  as  an  input.  The  temperature  of  hot  gases 
in  the  chamber  was  assumed  to  be  the  adiabatic  flame  temperature  of  product 
gases  from  the  Igniter.  An  estimate  of  heat  loss  indicated  that  it  is  less 
than  3%  of  the  chemical  heat  release;  therefore,  this  assumption  is  justified. 
The  flame  temperature  for  propellant  A  was  obtained  from  Price*'  of  NWC;  for 
the  ocher  propellant,  it  was  computed  from  the  CEC  72  thermochemistry 
program.'® 


TABLE  I  Propellant  Properties 


Propellant  Type 

A 

B 

Composition 

AP-hased 

AP/PBAA-EPON 

Weight  percent  of  oxidizer 

- 

75 

Average  particle  size,  d^p,  pm 

- 

76 

Pre-exponetlal  factor  in  Saint  Robert's 
burning  rate  law,  a,  mm/s/(atm) n 

1.62 

0.9591 

Pressure  exponent  in  Saint  Robert's 
burning  rate  law,  n 

.4108 

.41 

Flame  temperature,  ,  K 

3000 

1920 

Propellant  density,  kg/m^ 

1710 

1600 

The  calculated  pressure  distributions  at  various  times  for  a  typical 
test  case  for  Propellant  A  are  shown  in  Fig.  8.  Results  are  presented  only 
for  the  initial  rapid  pressure  rise  portion,  since  that  Is  the  region  of 
interest  for  this  study.  The  average  pressurization  rate  at  the  crack  en¬ 
trance  was  1.75  x  10J  atm/s.  During  the  initial  period,  a  traveling  ptessure 
wave  is  present  in  the  crack  cavity.  This  wave  moves  toward  the  crack  tip 
and  is  reflected  from  the  tip,  resulting  in  higher  pressure  near  the  tip  at 
t  »  0.6  ms.  As  time  progresses,  the  overall  pressure  in  the  crack  cavity 
continues  to  rise.  Thi?  is  the  result  of  a  combination  of  two  simultaneous 
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TABLE  tl  Numerical  Values  of  Input  Variables 


Variable 

Units 

Propellant  A 

Propellant  B 

X 

OB 

0.0 

0.0 

P 

*L 

on 

195.0 

195.0 

Sx 

S/kg 

0.0 

0.0 

b 

»3/kg 

1.0  x  10*3 

*3 

1.0  *  10 

kg/kaole 

76. io 

20.381 

Y 

1.21 

1.26 

Ppr 

kg/*3 

1710 

1600 

X 

kcal/s-a-K 

0.8  x  10~4 

0.3  x  10 

P* 

a 

■> 

o“/a 

0.18  x  10"6 

0.11  x  10 

P* 

cs/dh 

.025 

.025 

Tf 

K, 

3000 

1920 

T. 

K 

850 

850 

ign 

Xi 

K 

295 

295 

a 

(«*/s)/(ata)n 

1.62 

0.9591 

n 

- 

0.4108 

0.41 

K 

m3-R/kcal 

0.70  x  10~3 

• 

e 

3 

- 

53 

• 

T 

K 

1920 

1920 

c 

Gl  ' 

ata 

14.4  +  136  *‘-095c 

Sane  as  A 

X 

aca 

-n  — n 

47619 

47619 

K* 

(MPa)  p(a/s)  u 

- 

2  x  10  4 

"p* 

• 

- 

0.705 

n  * 
u 

*» 

- 

1.252 

“Defines  in  Sq.  (2d) 
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processes:  a)  a  continuous  increase  in  chamber  pressure  during  this  period, 
and  b)  mass  addition  Jue  to  burning  of  the  crack  walls.  After  about  0.8  ns, 

Che  pressure  in  the  crack  becomes  more  uniform  because  the  race  of  increase 
of  chamber  pressure  has  become  quite  small,  and  the  ignition  front  has  prop¬ 
agated  to  the  tip.  U  should  be  noted  that  as  time  progresses,  the  maximum 
pressure  in  the  creek  cavity  becomes  higher  chan  chat  in  the  chamber. 

Predicted  velocity  distributions  at  various  times  are  plotted  in  Fig.  9. 
At  t  •  0.1  ms,  velocity  in  the  latter  half  of  tha  crack  is  almost  zaro. 

Because  the  initial  favorable  pressure  gradient  causes  the  gases  to  scceler- 
ata,  the  gas  velocity  at  any  axial  location  increases  continuously  until 
t  *  0.5  ms.  Die  adverse  pressure  gradient  generated  by  the  compression  wave 
reflection  at  the  crack  tip  causes  the  gates  to  decelerate.  Since  velocities 
continue  to  decrease  after  0.5  ms,  only  one  velocity  distribution  between  0.5 
and  1.0  ms  is  shown  in  the  figure  to  avoid  confusion.  At  1.0  ms,  velocity 
in  the  entire  crack  is  negative,  i.e.,  gases  flow  out  of  the  crack  since 
pressure  throughout  the  cavity  is  higher  at  that  time  than  the  pressure  in 
the  chamber. 

Calculated  temperature  distributions  at  various  times  are  shown  in  Fig. 
10.  During  the  initial  period  (t  »  0.1  ms),  a  steep  temperature  gradient  is 
present.  Dils  is  due  mainly  to  the  fact  that  the  initial  Increase  in  the 
temperature  in  the  cavity  is  caused  by  hot  gas  penetration.  Since  the  gas 
velocity  in  the  second  half  of  the  crack  is  nearly  zero  at  0.1  ms,  the  tem¬ 
perature  in  that  region  is  also  equal  to  the  initial  temperature  of  the  gas. 

As  the  propellant  starts  to  burn,  the  temperature  in  the  crack  cavity  begins 
to  increase,  and  is  dominated  by  the  flame  tea*'-  rature  of  the  propellant. 

Since  the  flame  temperature  of  propellant  A  1$  .<igher  than  that  of  the  igniter 
gases,  and  also  there  is  a  compression  effect  in  the  cavity,  maximum  temper¬ 
ature  in  the  crack  is  higher  than  chat  in  the  chamber. 

Variation  in  the  port  height  (  =  half  of  the  gap  width)  at  various  times 
is  shown  in  Fig.  11.  The  initial  geometry  of  the  crack  is  shown  by  the 
dashed  line.  Because  of  the  rounding  of  the  crack  entrance,  the  gap  width 
at  the  entrance  is  large.  During  the  initial  period,  the  high  chamber  pres¬ 
sure  acting  on  the  crack  surface  exposed  to  the  chamber  causes  the  propellant 
near  the  crack  entrance  to  be  pushed  inside,  resulting  in  a  decreased  port 
height  near  the  crack  entrance.  At  t  ■  0.6  ms,  the  %  are  two  minimum  gap 
widths;  the  second  is  due  to  the  effect  of  compression  wave  reflection  in  the 
gas  phase  (see  Fig.  8).  As  time  progresses,  lovever,  the  pressure  in  the 
cavity  exceeds  that  in  the  chamber,  and  the  propellent  continues  to  burn;  the 
net  result  is  an  increase  in  the  port  height  at  10  ms. 

Comparisons  of  predicted  and  measured  pressure-time  traces  at  x  ■  48, 

138,  and  188  mm  are  shown  in  Figs.  12,  13,  and  14,  respectively.  The  pre¬ 
dicted  values  are  in  good  agreement  with  the  experimental  data.  The  steep¬ 
ening  of  the  pressurization  as  it  moves  along  the  crack  is  evident  from  both 
the  predicted  traces  and  experimental  data.  The  effect  of  the  increase  in 
pressurization  rate  is  most  pronounced  near  the  tip  of  the  crack.  This  is 
due  to  the  combined  effect  of  coalescence  of  the  compression  waves  and  steep¬ 
ening  of  the  pressure  front  caused  by  gasification  of  the  propellant  behind 
the  front. 


Crack  Uul(  ‘«ip  Width 


Non-dlaensloual  Axial  Distance,  x/L 

Figure  10  Calculated  Tesperacure  Distributions 
at  Various  T lisas 


Figure  11  Calculated  Variations  In  Fort  Height 
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Figure  12  Comparison  between  Predicted  and  Measured 
Pressure-Time  Traces  a:  x»i8  am 
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Comparison  of  predicted  and  measured  ignition-front  locations  for  var¬ 
ious  pressurization  rates  is  shown  in  Fig.  15.  Tt  can  be  seer,  that  the 
model  predicts  the  location  of  the  igniton  front  quite  well  except  very 
close  to  the  tip.  Near  the  tip  region,  the  model  overpredicts  the  deceler¬ 
ation  of  the  ignition  front.  This  is  believed  to  be  caused  by  the  one¬ 
dimensional  assumption  in  the  model.  The  one-dimensional  assumption  is 
inadequate  to  describe  the  flow  near  the  tip.  Since  the  calculated  veloci¬ 
ties  near  the  tip  are  small  (sea  Fig.  9),  both  the  gas  temperature  (see  Fig. 
10)  and  heat  transfer  coefficient  are  low;  heat  transfer  to  the  propellant 
is,  therefore,  quite  small.  In  an  actual  case,  the  presence  of  a  secondary 
flow  near  the  closed  end  will  produce  higher  rate  of  heat  transfer  to  the 
propellant,  resulting  in  less  pronounced  flame  front  deceleration.  It  is 
believed  that  Che  discrepancy  between  die  predicted  and  measured  ignition- 
front  locations  near  the  tip  region  of  the  crack  can  be  alleviated  if  the 
effect  of  enhanced  heat  transfer  near  the  closed  end  can  be  incorporated  into 
the  model.  It  should  be  noted  that,  under  cite  conditions  studied,  the  initial 
flame-spreading  process  is  not  substantially  affected  by  propellant  deforma¬ 
tion. 


Figures  16  and  17  show  calculated  pressure  distributions  at  various 
times  for  cases  in  which  propellant  deformation  becomes  very  important.  For 
results  shown  in  Fig.  16  (Case  1),  the  chamber  pressurization  rate  was 
3.8  x  10^  atm/s  and  the  gap-width  was  0.51  mm.  Only  those  curves  of  interest 
are  shown,  and  for  ease  of  explanation,  they  are  not  superimposed.  Between 
t  «  0.25  and  0.275  nr- ,  the  pressure  in  the  crack  shows  a  wavy  distribution 
which  corresponds  to  the  deformation  pattern.  During  this  period  the  crack 
gap  at  the  entrance  is  extremely  small  (partially  closed).  At  t  *  0.3  ms, 
high  pressures  in  the  chamber  and  at  the  crack  entrance  cause  the  partial 
closure  to  move  to  x/L  =  0.1,  the  gasification  in  che  minimum  flow  area 
region  results  in  a  local  pressure  peak.  The  localized  high  pressure  at 
x/L  *0.1  causes  the  propellant  to  deform;  the  region  of  partial  gap  closure 
movee^  to  x/L  0.2,  resulting  In  a  substantial  increase  of  local  pressure. 

The  subsequent  closure  results  in  a  zero  gap  width  and  causes  the  gas  dynamic 
solution  to  blow  up.  The  predicted  movement  of  the  partial  gap  closure  is 
qualitatively  similar  to  that  observed  experimentally  (see  Fig.  7). 

Case  2,  a  ra  in  Fig.  17,  was  computed  for  a  pressurization  rate  of 
5  x  10^  atm/s  and  an  initial  gap  width  of  .89  mm.  Before  0.475  ms,  the 
partial  gap  closure  occurs  near  the  crack  entrance.  Similar  to  case  1,  both 
the  high  pressure  near  the  entrance  in  the  cavipy  and  that  acting  on  the 
side  walls  of  the  propellant  exposed  to  the  chamber  cause  the  partial  gap 
closure  to  move  downstream.  This  area  reduction  and  the  continued  gasifi¬ 
cation  process  result  in  a  pressure  peak  near  x/L  ■  0.17  (at  0.475  ms).  At 
0.5  ms,  two  regions  of  partial  gap  closures  are  obtained,  resulting  in  pres¬ 
sure  peaks  at  x/L  =  0.1  and  0.27,  respectively.  As  the  pressure  peak  moves 
downstream,  another  gap  closure  develops  near  the  entrance  at  0.525  ms.  This 
type  of  process  continues  in  a  repeated  manner  as  shown  by  the  distributions 
at  t  *  0.525  ms,  0,55  ms  and  0.575  ms.  Once  again,  these  results  qualita¬ 
tively  explain  the  behavior  of  the  partial  gap  closures  observed  experimen¬ 
tally. 

Results  obtained  here  indicate  that  closure  of  the  crack  gap,  which  may 
initially  occur  at  the  crack  entrance  because  of  small  gap  widths  and  high 
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chamber  pressures,  may  propagate  along  the  crack  and  result  in  local  pres¬ 
sure  peaks.  The  gap  closure  is  not  observed  for  cracks  with  large  gap 
widths  or  at  low  chamber  pressures  because  propellant  deformation  under  these 
conditions  is  small  and  does  not  affect  the  ignition  process  substantially. 
Gap  closures  and  resulting  pressure  peaks  can  strongly  influence  the  convec¬ 
tive  burning  process  in  the  crack  and  may  contribute  to  the  def lagrat ion-to- 
detonation  transition  process. 

SUMMARY  AND  CONCLUSIONS 

This  research  program  was  undertaken  to  investigate  the  development  of 
convective  burning  in  solid  propellant  cracks  which  continually  deform  due 
to  burning  and  pressure  loading.  Both  theoretical  and  experimental  methods 
were  employed  to  study  the  flame  spreading  and  combustion  processes  in  pro¬ 
pellant  cracks.  In  the  theoretical  model,  the  effect  of  inter-related 
structural  deformation  and  combustion  phenomena  was  taken  into  account  by 
considering  a)  transient  one-dimensional  mass,  momentum,  and  energy  conserva¬ 
tion  equations  in  the  gas  phase,  b)  a  transient,  one-dimens tonal  heat  conduc¬ 
tion  equation  in  the  solid  phase,  and  c)  quasi-static  deformation  of  the  two- 
dimensional,  linear-viscoelastic  propellant  crack  due  to  pressure  loading. 
This  set  of  coupled,  nonlinear,  partial  differential  equations  was  solved 
numerically.  The  gas-dynamic  equations  were  solved  using  a  finite  difference 
analysis  and  the  structure  mechanics  equation  was  solved  using  a  finite 
element  analysis. 

Several  important  observations  and  conclusions  from  this  study  are  sum¬ 
marized  below. 

1.  Results  indicate  that  coupling  between  chamber  pressurization,  crack 
combustion,  and  propellant  deformation  is  quite  important,  especially  in 
the  case  of  very  narrow  cracks,  and  for  high  chamber  pressurization  rates. 

2.  Several  regions  of  partial  gap  closures  were  observed  experimentally 
in  narrow  cracks.  Calculated  results  show  that  these  partial  closures  may 
generate  substantial  local  pressure  peaks  along  the  crack.  The  gap  clo¬ 
sures  generally  do  not  occur  in  the  case  of  large  gap  widths  or  low  chamber 
pressures. 

3.  The  initial  flame-spreading  process,  or  the  time  for  the  ignition 
front  to  reach  the  crack  tip,  is  not  substantially  affected  by  propellant 
deformation.  However,  structural  deformation  may  significantly  affect  the 
pressure  distributions  along  the  crack, 

A.  Both  theoretical  and  experimental  results  show  that,  in  general,  the 
maximum  pressure  in  the  crack  cavity  is  higher  than  that  in  the  chamber. 

5.  Except  near  the  crack  tip  region,  predicted  results  for  the  ignition 
front  propagation  are  in  good  agreement  with  experimental  data.  Predicted 
and  measured  pressure  distributions  are  also  In  good  agreement. 
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Symbol 


a 


B 


x 

b 


c 

c 


f 


E 

E 


a 


1 


i 


P 

P 

r 


NOMENCLATURE 

Uescr  i£t JLon 

■> 

Cross-sectional  area  of  the  crack,  ra“ 

Spatial  change  of  cross-sectional  area  of  crack  with  respect 
to  the  axial  distance,  m 

Pre-exponential  factor  in  the  noneroslve  horning  rate  law, 
apn,  (mm/s)/(atm)n 

Body  force,  N/kg 
3 

Co-volurae,  m  /kg 

2 

Friction  coefficient,  2gTw/u 

Specific  heat  at  constant  pressure,  kJ/kg-K 

Hydraulic  diameter  of  the  crack,  m 

Total  stored  energy  (internal  and  kinetic),  kJ/kg 

Activation  energy  for  surface  reactions 

Relaxation  modulus,  kPa  or  atm 

2 

Local  convective  heat-transfer  coefficient,  kJ/nT-s-K 

Local  convective  heat-transfer  coefficient  over  the  propel¬ 
lant  surface,  kJ/m^-s-K 

Local  convective  heat-transfer  coefficient  over  non¬ 
propellant  port  wall,  kJ/m-s-K 

Bulk  modulus,  kPa  or  atm 

3 

Erosive  burning  constant,  m  -K/kJ 
Length  of  the  crack,  m 
Molecular  weight.,  kg/kmole 

Pressure  exponent  in  the  non-erosive  burning  rate  law 
Pressure,  kPa  or  atm 
Prandtl  number 
Burning  perimeter,  ra 
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?v  Wetted  perimeter  of  the  port,  m 

R  Specific  gas  constant  for  the  combustion  gases,  N-m/kg-K 

Re  Reynolds  number 

r^  Burning  rate  of  the  solid  propellant,  Including  the  erosive 

burning  contribution,  nm»/s 

T  Temperature  (without  subscript,  static  gas  temperature),  K 

Tgf  Average  film  gas  temperature,  (T  +  Tpg)/2,  K 

Tf  Adiabatic  flame  temperature  of  the  solid  propellant,  K 

T  ^  Initial  propellant  temperature,  K 

T  Propellant  surface  temperature,  K 

ps 

Tw8  Nonpropellant  wall  surface  temperature,  K 

t  Time,  s 

u  Gas  velocity,  m/s 

The  coefficient  in  viscosity-temperature  relation 
The  velocity  of  propellant  gas  at  the  burning  surface,  m/s 
x  Axial  coordinate,  m 

xL  Position  at  the  end  of  crack,  m 

Xp  Axial  distance  along  the  crack  at  which  propellant  begins,  m 

y  Perpendicular  distance  from  the  propellant  surface  into  the 

solid,  m 

Greek  Letters 

2 

a  Thermal  diffusivity,  m  /s 

8  Erosive  burning  exponent 

Y  Ratio  of  specific  heats 

5  Gap  width  of  the  crack,  m 

e  Strain  tensor 

eg  Surface  roughness,  w 

X  Thermal  conductivity,  kJ/m-s-K 
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0  Stress  tensor 

V  Gas  viscosity,  kg/m-s 

3 

p  Density  (without  subscript,  gas  density),  kg/m 

Characteristic  tine  of  the  gaseous  flame 

T  Transient  pressure  variation  time  (characteristic  time  of 

p  the  pressure  variation) 

Ts  Characteristic  time  of  the  unburned  solid  phase 

Shear  stress  on  the  port  wall,  kPa 

T  Normal  viscous  stress,  kPa 

6  Weighting  parameter 

6tf  Angle  measured,  in  a  counterclockwise  direction,  at  the  lover 

side  of  the  propellant,  degree 


Subscripts 

c  Rocket  chamber 

cri  Critical  condition  for  suface  ablation 

eff  Effective 

g  Gas 

1  Initial  value 

ign  Ignition  condition 

pr  Propellant 

ps  Propellant  surface 
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IMPROVED  PREDICTION  OF  FLAME  SPREADING  DURING 
CONVECTIVE  BURNING  IN  SOLID  PROPELLANT  CRACKS 


S.  M.  Kovaclc ,*  M.  Kumar,  and  K.  K.  Kuo* 
Department  of  Mechanical  Engineering 
The  Ier.nsyl varila  State  University 
University  Park,  Pennsylvania 


The  development  of  convective  burning  In  solid  propellant 
cracks  was  studied  both  theoretically  and  experimentally. 
Comparison  of  the  experimental  results  of  flame  spreading  with 
the  theoretical  predictions  of  the  previous  crack  combustion 
model  (Kumar  and  Kuo  1?30)  showed  good  comparison  except  near 
the  crack  tip.  The  model  was  modified  to  Include  a)  the 
effect  of  pressure-wave  reflection  near  the  crack  tip  on  local 
Ignition,  b)  the  effect  of  flow  recirculation  In  the  crack-tip 
region  on  local  heat  transfer  to  the  solid  propellant,  and 
c)  turbulent  transport  terms  In  the  f avre-averaged  conservation 
equations  along  the  crack.  Predicted  flame-spreading  rates 
obtained  from  the  modified  modes,  showed  good  agreement  with 
experimental,  data,  indicating  significant  Improvements  (near 
the  crack  tip)  over  predictions  of  the  previous  model.  The 
Inclusion  of  turbulent  transport  terms  along  the  length  of  the 
crack  did  not  significantly  affect  the  flame-spreading  process. 
Results  of  this  study  Indicate  that  the  effects  of  compression- 
wave  reflection  and  reclre  ilut Ion  are  important  In  the  crack 
tip  region. 


INTRODUCTION 

Solid  propellant  grains  sometimes 
develop  cracks  and  flaws.  These  cracks 
or  voids  may  be  caused  by  defects  in 
the  manufacturing  process,  thermal 
stresses  that  occur  during  curing,  the 
affect  of  aging  during  storage,  damage 
received  during  handling  and  transpor¬ 
tation,  and  thermal  and  mechanical 
stress  during  combustion.  The  cracks 
or  flaws  degrade  the  mechanical  pro¬ 
perty  and  structural  Integrity  of  the 
grain.  At  the  same  time,  they  provide 
additional  surface  area  for  combustion. 
Burning  within  cracks  may  result  In 
Increased  pressure  In  the  cavity,  the 
increased  pressure  may  ..auce  crack  pro¬ 
pagation,  development  of  strong  com¬ 
pression  waves  and  possible  generation 
of  shock  waves  which  could  cause  deton¬ 
ation. 

Some  of  the  processes  found  tw  be 
Important  [1-17]  in  ignition  and  com- 
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bustlon  of  propellant  cracks  are:  gas 
penetration  Into  the  crack,  convective 
heating  of  the  propellant  crack  surface 
prior  to  Ignition,  flame  propagation 
along  the  crack.  Increased  pressure  in 
the  crack  cavity  due  to  penetration  of 
product  gases  and  gasification  associ¬ 
ated  with  propellant  burning,  flew  re¬ 
versal  when  pressure  Inside  the  crack 
exceeds  that  of  the  main  chamber,  change 
In  crack  geometry  due  to  burning  as  well 
as  mechanical  deformation,  and  propaga¬ 
tion  of  the  crack  due  to  grain  fracture. 
The  governing  parameters  which  have  been 
found  to  Influence  convective  burning 
f I  —  i7 3  are:  chamber  pressure  and  pres¬ 
surization  rate,  g-oni'-try  of  the  crack, 
the  effects  of  erosl.‘<  burning  caused 
by  hlgh-vel wcity  gases  present  lr.  the 
crack,  physlcoc.html cal  properties  of 
the  propellant,  composition  and  flame 
temperature  of  igniter  gases  In  the  main 
chamber,  and  Initial  and  boundarv  condi¬ 
tions  to  which  the  crack  Is  subjected. 

One  of  the  Important  parameters 
which  exerts  a  stroriir  Influence  on  time 
variation  of  the  pressure  lirtribution 
In  the  crack  cavltv  1 the  instantaneous 
turning  surface  area.  Since  the  flame- 
Sj  reading  rate  determine.;  the  Irnlt-'d 


PASS  IS  EBST  QUALITY  PRACTI 
UJUUSttHD  10  UWJ  - 


fraction  of  the  total  surface  area, 
accurate  determination  of  the  flame- 
propagation  rate  is  important  for 
better  predictive  capability  of  a  mo¬ 
del.  Even  though  the  present  analysts 
is  limited  to  a  non-propagating  crack, 
Improved  prediction  of  the  flame- 
spreading  rate  would  be  helpful  in  dev¬ 
eloping  models  that  Include  effects  of 
crack  propagation  and  branching.  Fur¬ 
thermore,  a  detailed  analysis  of  the 
flame  spreading  is  desirable  In  order 
to  understand  some  of  the  mechanisms 
associated  with  the  multifacet  problem 
of  deflagration  to  detonation  transi¬ 
tion  (DDT). 

In  the  past,  several  researchers 
have  investigated  convective  burning  In 
propellant  cracks.  Prentice  [1,2], 
Margolin  and  tiargulis  (31,  Godal  [A], 
Bobolev  et  ai.  [5],  Payne  [6],  and 
Krasnov  et  al.  [7]  have  studied  the  on¬ 
set  of  convective  burning.  Preliminary 
Investigations  of  the  development  of 
convective  burning  were  conductod  by 
Belyaev  et  al.  [3,9].  Both  theoretical 
and  experimental  investigations  of  the 
development  of  convective  burning  in 
solid  propellant  cracks  have  been  con¬ 
ducted  at  The  Pennsylvania  State  Uni¬ 
versity  [10-15].  Detailed  reviews  of 
previous  research  in  this  area  are 
giver,  in  Refs.  11,  It,  arid  17.  The 
crack  combustion  model  developed  by  Kuo 
et  al .  [10,11]  provided  the  basis  for 
tne  ;r«ser.t  theoretical  study. 


culatlng  flow  near  the  tip  region  wl  !.l 
cause  the  heat-trirsfer  rate  to  be  l igh- 
er  than  that  predicted  by  the  one- 
dlnenslonal  model.  A  recent  study  [lcj 
found  that  the  effect  of  compression 
wave  reflection  and  recirculation  near 
the  crack  tip  is  Important  for  ignition 
of  the  crack  tip.  In  some  of  the  exper¬ 
iments  [151,  it  was  found  that  the  crack 
tip  Ignited  before  the  convective  igni¬ 
tion  front  reached  the  tip,  leading  to 
the  conclusion  that  heat  transfer  pro¬ 
cess  near  the  tip  is  different  from 
that  along  the  crack.  It  is  clear  that 
these  complex  processes  near  the  crack 
tip  must  be  taken  Into  account  to  im¬ 
prove  prediction  of  flame  spreading. 

This  study  was  directed  toward  Im¬ 
proving  the  predictive  capability  of  the 
convective-burning  theoretical  model 
[10,11]  by  including  a)  the  effect  of 
pressure-wave  reflection  at  the  crack 
tip  on  the  ignition  process,  b)  an 
approximation  to  account  for  the  effect 
of  recirculating  flow  near  the  crack 
tip  on  local  heat  transfer  to  propellant 
surface,  and  c)  turbulent  transport 
terms  in  conservation  equations  along 
the  crack.  The  objectives  of  this  study 
are : 

1.  To  improve  the  predictive  capability 
of  the  convectlve-burnlr.g  model  so 
that  the  rate  of  flame  propagation 
in  a  solid  propellant  crack  can  be 
accurately  predicted. 


This  mouel  [lu,ll]  considers  tran¬ 
sient  one-dlmenslonal  mass,  momentum, 
and  energy  conservation  equations  in 
the  gas  phase.  A  one-dlmenslonal  tran¬ 
sient  heat  conduction  is  considered  in 
the  3olid  phase.  As  the  gas  flows  over 
the  crack  surface.  It  heats  the  propel¬ 
lant.  Flame  spreading  is  postulated  to 
result  from  successive  ignitions  of 
adjacent  surfaces  by  convective  heat 
transfer.  A  comparison  of  predicted 
results  from  the  Kuo  model  [10,11]  and 
experimentally  measured  [13 ]  lgnltlon- 
fror.t  locations  Indicated  that  the 
model  tends  to  ur.derpredlct  the  flame¬ 
spreading  rate  near  the  tip  region  [11, 
l1*].  It  appears  that  one  of  the  main 
reasons  for  this  discrepancy  Is  the 
one-dimenslor.al  nature  of  the  model. 

The  assumption  of  one-dimensional  flow 
13  adequate  along  the  crack;  however, 
it  cannot  describe  the  flow  recircula¬ 
tion  near  the  crack  tip.  Predicted 
velocities  near  the  crack  tip  are 
small,  resulting  In  low  gas  temperature 
and  heat  transfer  coef f lcient ;  conse¬ 
quently,  the  predicted  rate  of  heat 
transfer  to  the  surface  is  small. 

In  reality,  the  presence  of  redr- 


2.  To  compare  the  predictions  of  the 

modified  crack  combustion  model  with 
experimental  data  or  flame  spread¬ 
ing  and  pressurization  In  prouellant 
cracks . 


ANALYSIS 

The  present  model  was  developed  to 
simulate  a  transverse  propellant  crack 
inside  a  rocket  chamber,  as  shown  in 
Fig.  1.  Hot,  high  pressure  gases  from 
the  main  chamber  penetrate  the  crack  and 
transfer  heat  to  unhurried  prooellant 
surfaces.  The  rate  of  convective  heat 
transfer  from  the  hot  gases  to  the  pro¬ 
pellant  is  determined  by  the  Instantan¬ 
eous  values  of  gas  velocity,  tempera¬ 
ture,  density,  and  local  propellant  sur¬ 
face  temperature.  As  the  gases  continue 
to  flow  into  the  crack,  gas  pressure, 
temperature,  and  density  along  the  crack 
Increase.  After  being  heated  for  a  per¬ 
iod  of  time,  the  propellant  surface  will 
begin  to  gasify  or  Ignite,  generating 
more  hot  gases.  The  flame  will  spread 
along  the  crack  until  the  -ntlre  crack 
is  Ignited.  It  is  possible  tnat  the 
local  Internal  pressure  may  exceed  the 
main  clumber  pressure,  causing  flow 
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reversal,  with  hot  gases  generated  by 
the  Ignited  propellant  being  forced  out 
of  the  open  end  of  the  crack. 


Fig.  i  Schematic  of  a  Solid  Propellant 
Crack 

The  following  basic  assumptions  are 

employed  in  the  analytic  model: 

(1)  Bulk  flow  of  ga3es  along  the  pro¬ 
pellant  crack  is  quasi -one-dimen¬ 
sional. 

(2)  Gases  flowing  into  the  propellant 
crack  obey  the  Hobel-Abel  law  for 
Imperfect  gases. 

(3)  The  flame  zone  Is  small  compared 
to  the  gap  width  or  the  crack. 

(4)  Rate  processes  Ir.  the  gaseous 
flame  are  considered  qua3i-steadv ; 
that  Is,  the  characteristic  time 
associated  with  the  gas  phase 
reactions  Is  short,  compared  to 
that  of  the  transient  pressure 
varlat Ion. 

(5)  Deformation  of  the  solid  propel¬ 
lant  is  neglected  during  the 
flame-sreading  process.  Any 
changes  in  the  geometry  of  the 
crack  are  due  solely  to  ablation 
of  the  propellant  surface. 

(6)  Thermal  properties  cf  propellant 
are  constant. 


composition  technique  [19]  was  used.  Ne¬ 
glecting  the  third-order  turbulent  cor¬ 
relations,  the  three  Favre-averaged 
conservation  equations  can  be  written 
as  follows: 


Mean  mass  conservation  equation: 

it  (pV  *  h  QV  "  rbcPpb  (1) 

The  mean  momentum  conservation  equa¬ 
tion: 

ft  av  +  h  <p  r,v  +  h 

-  -  h  (F  V  +  lx  cj*  AP  h  <^>1 

+  pk(V  ~\ r  a'Qi  9w 


*  Fapb*  -  (rbpPpb)  V  3ln  ®w  (2) 


The  mean  energy  conservation  equa¬ 
tion: 


fc  (ap  * *  r5  +  \  —  i) 
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Most  of  these  assumptions  are  Jus¬ 
tified  for  crack  combust  Ion  studies 
[11,14].  Even  though  the  propellant  Is 
deformable,  it  has  been  found  that  de¬ 
formation  does  not  significantly  In¬ 
fluence  the  flame-spreading  proce30 
[143. 

Based  upon  these  assumptions,  th*» 
quaal-one-dlmenslonal ,  unsteady  mass, 
momentum,  and  energy  equations  for  the 
gas  phase  were  obtained  in  terms  of  the 
Instantaneous  variables  [18].  Instan¬ 
taneous  gas-phase  equations  were  con¬ 
verted  Into  a  set  of  mean-oenservat icn 
equations  by  uaing  an  averaging  techni¬ 
que  for  turbulent  flows.  Since  the  mas 
is  highly  Incompressible,  the  Fuvre  de¬ 


Thv  equation  of  state  for  the  cas 
phase : 

r  (i  -  m  «  nr  1 4 ) 

Further  approxima' Ions  are  needed  in 
order  to  achieve  turbulence  closure  of 
fhe  model.  In  an  attempt  to  keep  the 
model  3imrle,  and  at  the  same  tine  to 
take  Into  account  most  of  the  Important 
physical  processes,  the  following  rela¬ 
tionships  were  used: 
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These  relationships  coulu  be  r-garued 
aa  the  dtf Initi  unc  fur  11.  and  k,.  h. 
tnialtiop,  3otae  further  simplifications 
were  made;  these  Include 

a  >>  a",  i  »  T"  ('> 

ar.a 

p(Cv T  ♦  -|  )  :»  \  p.i-u*  (B) 

The  three  mean  eouservat  l.u  .-.tu-tt  Ions 
were  rearranged  tu  firm  -»  set  uf  gov¬ 
erning  equations  .fit*  calculating  gaa- 
phaue  temperature,  pressure,  anl  velo¬ 
city. 

The  pro;  ellar.t  surface  tern;  eras  ire 
at  any  location  along  the  crack  before 
attainment  of  ignition  la  calculated 
by  solving  the  transient,  one— dimen¬ 
sional  heat  conduction  equation  in  the 
same  manner  as  .riven  It.  Refs.  10  and 
11.  A  two-temperature  ijnltlvi:  crit¬ 
erion,  simulating  onset  of  ablation 
and  Ignition,  Is  used  to-  let « rain e  sur¬ 
face  gasification  and  bi;i“«.li.,r  candl- 
11  aim.  '  The  temperature  profile  In  the 
soil  i  propellant  -Is  solved  until  igni¬ 
tion  Is  achieved. 

empirical  correlations  are  sup- 
pllei  for  the  heat -t ran.: for  coeffici¬ 
ent  ,  drag  coefficient  propellant  burn- 
in. ‘-file  law,  etc .  d juries  of  the 
correlations  used  In  no  lei  at--.-  1*’“ 

l s el  In  Table  1. 

Table  1.  Empirical  tiorrelatlui.s 
Used  In  the  Mode  1 

Parameter  tlourCo 

lieat  transfer  Modified  I)lltu.i-3oel- 

coefflclent  ter  correlation  f 20 1 

Has  viscosity  live  h  la  express  lor.  [21] 

frandtl  number  Evehla  formula  [21] 


t 10, 11 J .  When  the  gases  enter  the 
crack,  boundary  conditions  are: 

F(t,0)  -  Pc(C)i  T(t ,0)  -  Tfl(t); 

u(t,xL)  •  0  (10) 

When  the  gases  flow  out  of  the  crack 
subsonlcally ,  boundary  conditions  are: 

Ft  t,0)  -  Pc ( t ) ;  u(  t ,xL V  •  0  (11) 

When  gases  exit  the  crack  supersonic¬ 
ally,  the  boundary  condition  la  that 
mean  ga3  velocity  at  the  crack  tip  Is 
zero . 

The  equations  and  boundary  condi¬ 
tions  given  above  describe  gas  flow 
along  the  crack.  As  mentioned  in  the 
Introduction,  heat-transfer  processes 
near  the  tip  region  should  include 
the  effect  of  pressure-wave  reflection 
and  flow  recirculation  in  order  to  Im¬ 
prove  prediction  of  the  f lame -propaga¬ 
tion  rate.  As  reported  by  Kumar  and 
Kuo  [12],  the  effect  of  pressurization 
rate  and  pressure-wave  reflection  Is 
Important  for  heat  transfer  at  the  end 
of  a  crack  and  for  local  Ignition. 

-lince  the  flow  conditions  and  crack 
configuration  of  this  siuuy  are  similar 
to  those  presented  Iri  Ref.  12,  the  tlp- 
Jgnitlcr.  model  [12]  la  compatible  for 
direct  Incorporation  into  the  pre¬ 
sent  theoretical  model.  In  the  tip- 
Igultlon  model,  the  propellant  at  the 
close i  end  is  considered  to  be  therm¬ 
ally  thick.  The  propellant  surface  at 
the  tip  is  heated  by  hot  gases  driven 
into  the  crack  as  a  result  of  rapid 
chamber  pressurization  and  the  presence 
of  strong  pressure  gradient  forces. 

The  tlp-lgnltlon  model  describes  the 
process  of  heating  to  ignition  at  the 
closed  end  of  the  propellant  crack; 
accordingly ,  this  model  is  Incorporated 


Friction 
coe  f  f  1  c  1  er.t 


Modified  Ooletrook 
correlatlon  [22,23] 


Strar.u  burning  dalnt-Rohert  .;xpres- 
rute  Sion  [2-i 

Erosive  turning  Lenoir-Kcbl  Hard  [251 
rate  or 

Ruzdan-Kuo  formula 

_ [26] 

Initial  conditions  f-r  U.c  cus- 
phase  equation  describe  a  uniform  dis¬ 
tribution  in  the  crack: 

u  ( 0 ,  x  )  •  u  t;  T  ( 0 ,  x  )  •  Tj  ; 

f  i  i ,  x )  *  F  ^  ( 9  ) 

Boundary  conditions  depend  up.,n  flow 
c.u.litic-ns  at  the  ..•rack  entrance 


for  calculations  at  the  tin. 

Although  rec lrculatin.-,  flows  in 
gaps  have  been  researched  quite  exten¬ 
sively  [27,23],  most  researchers  have 
only  studied  gaps  with  small  length-to- 
dlumeter  ratios  (In  the  order  of  5). 
Cracks  with  iongth-to-diameter  ratios 
in  the  range  of  30  to  200  were  used  In 
the  present  study.  Kuo  et  al.  [29] 
conducted  flow-visualisation  experi¬ 
ments;  on  gaps  with  aspect  ratios  as 
high  us  29.5,  but  did  not  give  special 
attention  to  recirculation  near  the  tip 
of  the  crack. 

As  an  attempt  to  improve  the  flame 
spro.i.Ui.r,  prediction  in  the  tip  region 
with'  i\  having  to  adve  a  conplex  multi- 
dime:..-  i  dual  problem,  a  simple  method 
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for  approximating  the  rate  of  heat 
transfer  caused  by  flow  recirculation 
was  Incorporated  into  the  theoretical 
formulation,  Aa  a  first  approximation, 
a  condition  of  linear  heat  flux  to  the 
propellant  was  Imposed  between  a  point 
In  the  crack  at  the  edge  of  the  recir¬ 
culation  region  and  the  tip  of  the 
crack.  The  heat  flux,  at  any  point 
before  the  region  of  recirculation, 
was  calculated*  using  the  temperature 
difference  between  the  gas  and  solid 
phase  and  the  heat-transfer  coefficient 
baaed  upon  the  modified  Dlttus-Boelter 
correlation  [20],  The  heat  Tlux  at  the 
crack  tip  was  calculated,  using  the 
tlp-lgnitlon  model  [12].  The  heat  flux 
at  any  point  in  the  recirculation  reg¬ 
ion  was  determined  by  Interpolating 
between  the  flux  at  the  crack  tip  and 
that  at  the  edge  of  the  recirculation 
region.  The  initial  selection  of  the 
size  of  the  recirculation  region  waa 
based  upon  a  good  match  between  experi¬ 
mentally  measured  flame-spreading  rates 
and  theoretically  calculated  resulta 
for  one  test  case.  It  was  found  that 
the  location  of  the  edge  or  the  recir¬ 
culation  zone,  so  long  as  the  length 
of  the  recirculation  zone  was  about 
l/5th  of  the  crack  length,  does  not 
significantly  afreet  results. 

The  numerical  technique  for  solv¬ 
ing  the  gas  and  solid-phase  governing 
equations  along  the  crack  employs  an 
Implicit  scheme,  similar  to  that  given 
in  Refa.  10  and  11.  The  crack  combus¬ 
tion  model  calculates  gas  pressure  and 
temperature  In  the  region  adjacent  to 
the  crack  tip,  and  the  tip  Ignition 
model  calculates  heat  flux  and  surface 
temperature  of  the  propellant  at  the 
crack  tip.  Numerical  solution  or  crack 
combustion  model  equations  and  tip  ig¬ 
nition  model  equations  Is  carried  out 
separately,  but  concurrently. 

RESULTS  ANC  DISCUSSION 

An  experimental  investigation  war 
conducted  In  order  to  validate  results 
obtained  from  the  theoretical  model. 

A  teat  rig  was  designed  tc  simulate, 
as  closely  as  possible,  actual  rocket 
chamber  conditions.  An  igniter  system 
incorporating  a  solid  propellant  charge 
waa  used  as  the  source  of  hot  gas  gene¬ 
ration.  High-speed  photographs  and 
quick  response  pressure  transducers 
were  used  to  observe  name  preparation 
and  to  measure  pressure  along  the 
crack.  Detailed  discus jlons  of  the 
experimental  setup  and  prove  lure  oar. 
be  found  in  Refs.  11  and  12.  Two  types 
of  AP  composite  solid  propellants  were 
atudled.  Table  2  lists  some  of  the 
important  propellant  properties. 


Table  2  Prctellan*  Proeerttes  and 
Input  Parameters 


Parameter  (Unit) 

Prorellant 

frcrel Innt 

A 

3 

Propellant 

AP-baned 

AF/F5AA- 

Composition 

eio:: 

a  (mm/s-atmn) 

1.62 

0 . oh si 

n  (-) 

0.  AlOR 

0  .  'i  1 0o 

Tr  (K) 

’CO" 

)0?0 

Pp  (kg/m?) 

17)0 

IhOO 

(W/m-K) 

P  ? 

Op  (mvs) 

(kg/kmole) 

3  •  3r 

0.l9xl0-'' 

26.1 

?.l 

0. 11x10“* 
20.  35 

Y  (-) 

1.21 

1.26 

b  (mVkg) 

0.001 

o .  00 1 

Vn  «> 

:>0 

<?Sh 

uerf/u  (-) 

100 

100 

Ref.  18 

keff/k  (-> 

.'OO 

inn 

Ref.  12 

Flr>»r“a  a  .and  3  Usplay  two  se^-s 
of  pressure  traces  rrem  test  firings  in 
which  all  teat  conditions  except  th» 
crack-map  width  w*re  approximately  the 
same,  nap  widths  for  Tost  (rlc.  2) 
and  63  (Fir.  3''  v'r*  1.27  and  mm, 

respectively.  i-"ir  rrersiir"*  gauro? 
(qj-'iF)  wer-  uro,j  t;  noas.i!",  ppnsruro 
at  four  axial  locations  along  the  crack. 
Several  differences  can  be  seen  from 
these  pressure  traces  present  In  the 
two  tests.  In  both  tests,  the  rate  of 
pre3surlsatlcn  increases  consecutively 
at  eacji  downstream  location  from  the 
crack  entrance,  but  th«  crack  with 
smaller  gap  width  experiences  a  hither 
peak  pressure.  Furthermore,  the  high¬ 
est  peak  pressure  occurs  in  the  inter¬ 
ior  of  the  narrow  crack,  whereas  peak 
pressure  is  much  mor*  uni  fern  in  the 
wider  crack.  This  can  b»  «xplained  by 
noting  that  in  bi  *  h  cracks  the  burring 
surface  area  is  identical,  l at  * h"  rcrv 
area  In  the  ?'orr  w  ere  ’!■:  *s  rrorert  lon- 
ally  S'nnll-r  than  thot  'f  the  v.dt* 
crack.  As  ‘h'-1  rrvp'llant  alone  t):« 
crack  brains  to  islfy,  local  pr^ssur* 
inside  t. he  crack  r.av  become  -reader 
than  chanter  nr^rmre,  causing  ras  flow 
to  reverse  and  bcr'.n  to  exit  the  crack. 
The  gaseous  macs  generated  tv  the  burn- 
ir.p  prcp*llJir.v  ~p»*civ.  *?r 

than  the  mv*r  which  is  abl“  t c  exit 
fr*!*.  *:  ft  *  *  ft  »  t  \r\  *“ v  A 

r r “ f  2 u t4 jw.'i  ifl  *  j ,.n  r» *» n ** y  ft u i  I  5 r  \  j 
a  maximum.  ror  t‘-c  narrow  orach,  the 
smaller  port  sr-a  will  restrict  “h- 
exiting  flow  ?0'*i  ‘*r.  porultir-  In  a 
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higher  peak  pressure  Ml  chin  th*  crack. 
Higher  preaaure  In  the  interior  of  the 
narrow  crack  alao  causes  the  burning 
rate  to  Increase,  thua  further  increas¬ 
ing  the  preaaure.  Eventually,  the 
preaaure  In  the  Interior  la  higher  than 
the  preaaure  at  the  crack  opening. 

Propellant  B 


0.0  aa 


■  »  ■  ■  ■ 

0.0  0.S  1.0  l.S  2.0  U 1 - 1— J 

15  4S  75  105 

Tine,  m* 


Fig.  2  Heaaured  Pressure-Time  Tracea 
for  Teat  62 


Tiae,  aa 

rig.  3  Heaaured  Praaaure-Tiee  Tracea 
for  Teat  63 


Another  Important  difference  be¬ 
tween  Test  62  and  C 3  1'-  the  effect  of 
crack-gap  width  on  preaaure  waves  tra¬ 
veling  along  the  crack.  Preaaure 
truces  from  Test  62  (shown  In  Fig.  2) 
Indicate  the  presence  or  a  strong  pres¬ 
sure  wave  traveling  along  the  crack, 
resulting  In  the  local  peak  and  valley 
at  05.  The  local  peak  Is  caused  by  the 
compression  wave  reflection  at  the 
crack  tip,  and  the  valley  Is  believed 
to  be  the  result  of  expansion  waves 
generated  by  the  reflection  of  the  com¬ 
pression  wave  at  the  crack  entrance. 
Preaaure  traces  from  Test  63  (shown  In 
Fig.  3)  Indicate  no  strong  pressure 
waves  moving  along  the  crack.  This  can 
be  explained  by  considering  the  mechan¬ 
ical  deformation  of  the  propellant  along 
the  crack.  As  gases  enter  the  crack, 
pressure  waves  traveling  toward  the 
crack  tip  coalesce  and  become  stronger. 
In  the  region  close  to  the  pressure 
wave  front,  the  propellant  walls  of  the 
crack  deform.  Increasing  the  local 
crack  volume  and  lowering  the  local 
pressure  In  the  region  near  the  pressure 
wave  front.  In  a  very  narrow  crack, 
mechanical  deformation  of  the  propellant 
due  to  pressure  loading  can  be  very 
Important.  Propellant  deformation  for 
a  given  pressure  will  be  approximately 
the  same  for  the  same  web  thickness, 
but  the  percentage  change  lti  the  port 
area  will  Increase  tis  the  original 
crack  t*u p  width  decreases .  The  gap- 
width  of  Test  r>3  Is  2.5  times  smaller 
than  that  of  Teat  62;  hence  the  pres¬ 
sure  wave  effect  la  much  lesj  pronounc-d 
in  Test  6j. 

The  effect  of  chamber  pressuriza¬ 
tion  rate  on  the  p  res  sure- 1  l:r.e  history 
of  the  crack  can  be  3een  by  comparing 
two  sets  of  pressure  traces  from  Tests 
6k  and  6\.  as  shown  In  Figs.  k  and  5. 
These  two  tests  were  conducted  with 
the  same  propellant  (Propellant  A)  and 
same  crack  geometry  (L  ■  1  ( > 5  mm,  5  * 

J.89  nun),  but  with  different  chamber 
pressurization  rates.  Test  65,  with  a 
chamber  pressurization  rate  of  about 
A. 5  x  10*  atm/s,  exhibits  a  strong 
pressure  wave  effect  along  the  crack, 
while  Test  6k,  with  a  chamher^pressurl- 
zatlon  rate  of  about  1.3  x  10-  atm/s, 
shows  no  significant  pressure  wave 
effect  In  the  crack. 

Figures  6  to  9  pres-.-r.t  a  compari¬ 
son  of  predicted  and  measured  flame 
front  locations.  The  sash-  l  lines  r*- 
pro-ser.t  predicted  flame  fr  -n*  ! -ar.  1 
obtait.es  J  r  or.  the  prv.-l  ■  r.viel 
11];  tic.  solid  lines  r-r.sivn*  .1  <■ .; 
obtain- 1  from  the  present  study.  The 
plots  si. /w  the  time  required  for  the 
Ignltl  n  front  to  reach  a  particular 
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location  versus  the  nondimens tonal 
axial  distance  alone  the  crack.  '3ut- 
stnntial  improvement  In  the  prediction 
of  the  flame  front  near  the  crack  tip 
region  can  be  aeen  from  theae  figures. 
The  effect  of  Various  parameters  on 
flame  front  propagation  Is  Riven  In 
Refs.  11  and  1?. 


Non-dimensional  Axial  Distance,  x/t 


Trope 11 ant  A 
L  •  19S  mi 


'5  75  105 


Time,  as 

Fig.  3  Measured  Prea aura •Time  Traces 
for  Teat  65 


Fig.  6  Comparison  of  Predicted  and 

Measured  tgnltion  Front  Propa¬ 
gation  for  Test  62 


0  0.:  0.4  0.6  0.8  i.O 


Non-dimensional  Axial  Distance,  x/L 


Fig.  7  Comparison  of  Predicted  and 

Measured  Ignition  Front  Propa¬ 
gation  for  Test  6> 
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A  Measured 

-  Theory  (previous) 

— — '  Theory  (present) 

Propellant  A 
L  •  IPS  mi 
5  *  0.39  a* 

<H»c/«it  ■  4.5x10*  atn/s 


0.2  0.4  O.o  0.3  1.0 

Non-dimensional  Axial  Distance,  x/L 


Fig.  8  Comparison  of  Predicted  and 

Measured  Ignition  Front  Propa¬ 
gation  for  Test  6S 


O  Measured  ! 

-  Theory  (previous)  | 

— —  Theory  (present)  i 


Propellant  B 
b  «  97.5  m 
S  •  C.89  »  , 

dPc/dt  •  1.5x10*  ata/s 


0.2  0.4  0,6  0.5  1.0 

Non-dleensional  Axial  Distance,  x/L 


Fig.  9  Comparison  of  Predicted  and 

Measured  Ignition  Front  Propa¬ 
gation  for  Test  7 A 


’"ho  predicted  pressure  variations 
with  tvspeat  to  time  and  axial  fliatanoe 
f  r  Teat  ?*t  la  shown  in  a  three-dimen¬ 
sional  plot  in  Fir.,  lo.  (Flame-spread¬ 
ing  results  or  ttis  test  are  shown  in 
Pig.  9).  In  Pig-10  v  tfce  traveling 
pressure  wave  front  oan  clearly  be  seen 
moving  back  And  forth  aloitit  the  crack 
as  time  progresses.  The  caloulated 
pressure- tine  trace  la  In  reasonable 
agreement  with  the  measured  data  [18], 
The  physical  refason  for  ths  presence 
or  crests  and  troughs  on  the  plot  Is 
the  some  as  that  given  earlier  in  the 
explanation  of  pressure  peaks  and  val¬ 
leys  present  in  Fir.3.  2  and  5. 

Propel  Ism  B 
L  *  97. S  m 
6  ■  0.89  ■■ 


"  0.4  \  K*“V 

Fig.  10  Predicted  Pressure  Distribu¬ 
tions  Versus  Tima  for  Taat  74 

Predicted  velocity  variations  with 
respect  to  time  and  axial  distance  for 
Te3t  71*  are  also  shown  In  a  three- 
dimensional  plot  In  FI.*.  11.  During 
the  initial  period,  the  favorable 
pressure  gradient  causes  the  hot  gas 
to  penetrate  the  opening  portion  of  the 
crack.  As  the  propellant  begins  to 
burn  and  the  flame  continues  to  spread, 
local  pressure  and  temperature  increase, 
further  facilitating  gas  penetration 
deeper  into  the  crack  and  increasing 
gas  velocity.  As  the  pressure  in  the 
cavity  increases,  the  favorable  pres¬ 
sure  gradient  effi  t  diminishes,  result¬ 
ing  in  a  decrease  .n  «:a a  velocity. 

When  t  la  about  0.it5  ms,  the  adverse 
pressure  gradient  in  the  crack  becomes 
strong  enou-h  to  reverse  flow;  there¬ 
after,  gases  produced  In  the  crack 
cv.'ity  flow  out  of  the  crack.  Velocity 
ulutritut lor.  la  generally  JletatM  b? 
pressure  distribution. 

WI18  PADS  13  ES3T  tyMXOT  ISMS TIGABW 

,  QOdH  iwsxshkd  »  DW  - 
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Propellant  1 
l  >  97. S  m 
i  m  o.n  m 


Fig.  11  Pr«dlct«d  Velocity  Olatri 
butiona  Varau a  Tlaw  for 
TMt  74 


Rosults  of  tho  numerical  compu¬ 
tation*  show  that  although  gas  now 
In  tha  propellant  crack  is  highly 
turbulent  and  high  temperature  gradl- 
•nta  are  present,  as  expected,  the 
effect  of  turbulent  transport  In  the 
axial  direction  was  not  significant 
when  compared  with  other  terrr.3  in 
the  energy  equation  [l8j.  This  is 
due  to  the  fact  that  the  riow  la 
highly  convective,  and  the  energy 
transfer  In  the  axial  direction  Is 
dominated  by  convection  caused  by 
high  gas  velocities.  The  diffusion 
of  energy  in  the  axial  direction  by 
turbulent  transport  Is  much  smaller 
than  that  by  convective  energy  trans¬ 
port. 

CONCLUSIONS 

Several  Important  observations 
and  conclualons  from  this  Investi¬ 
gation  are  summarized  he low. 

(1)  Predicted  Ignition  Trent  pro¬ 
pagation,  using  the  current 
theoretical  model,  was  In  very 
good  agreement  with  experi¬ 
mental  data,  Indicating 
significant  Improvement  over 
predictions  of  the  previous 
model  near  the  crack  tip. 

(2)  The  effects  ot‘  compression  wave 
reflection  and  now  recircula¬ 
tion  are  Important  In  consid¬ 
eration  oT  heat  transfer  between 
gas  and  solid  propellant  near 
the  crack  tip. 


Mi  Beth  experimental  and  theoretical 
results  exhibit  a  traveling  proa- 
sure  wav*'  phenomenon.  This  phe¬ 
nomenon  Is  more  pronounced  with 
tnerrases  In  crack  gap  width  or 
chamber  rressurixation  rate. 

(H)  Since  the  axial  transport  of 
energy  In  the  gas  phase  is 
strongly  dominated  by  convection, 
the  inclusion  of  axial  turbulent 
transport  terms  In  the  governing 
equations  does  not  afreet  the 
flame-spreading  rate  along  the 
crack. 
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nomenclature 


A  Cr^ss-rectlonal  area  of  the 

p  crack,  m' 

n  rro— exponential  factor  in  the 

nonpros Ive  burning  rate  law, 
arn .  (mrt/s  V(ntm)rt 
B  Body  force,  5/kg 

b*  lu:r«  ,  irM'kr 

Cp  ."«'•>  •  if’. !'«at  at  constant 

1  pr*»pr  ir*1 ,  J/kg-K 

Cv  Np^'iric  heat  at  corst.ant 

vol>,p»,  J/kg-K  , 

f  Krlot»>n  e 'of  fl  ■’lent ,  2^/ru1- 

fi  Local  I'onveotive  heat-transfer 

coefficient  over  ^be  propel¬ 
lant  rurface,  W/m*-K 
R  Local  convective  heat- transfer 

cw  cc*>ri'l  -  lent  over  non-propellant 

p-r*  wall,  W/mc-K 

k  Thermal  conductivity,  W/m-K 

k  f  k  *  k. ,  W/m-K 

L  r.enr-h  of  the  crack,  n 

M  Molecular  weight,  ke/kmole 

nw  Treasure  exponent  in  the  non- 

en.'slve  burning  rate  law 
T  T  r->spure ,  Ta 

r  Burning  perimeter,  m 

P“  Wot t.. i.l  porlmoter  of  the  port,  m 

pw  Free  if!  c  gas  con.'-tant-  for  the 

c">mburt !  on  rases ,  N-m/kg-K 
!•  I'urnlnr  nb'  of  the  solid  pr»'- 

^  w  ,e  '  J  “  m  »•  ^  1  r?*  ^ta  9^91 


turning  "■  ntrlbut l^n .  m/s 
T  T  -'r-ature  (without  subscript, 

ctati-'  ms  temperature ) ,  K 
T..  A'H’ba*-lc  flame  temperature 

‘  of  '  I  rolld  propellant,  r 

T  Tr-pe  I  .inr-  surface  temperature. 

PC  y 

?  Il'imruio  11. ant  wall  surface 

WF  t  a’  ure  ,  K 


5SS5SSK* 


PKAfiXlflOU 


» 
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P 


i 


i 

i 

i  • 


fif  Tli.'  velocity  >.P  prcpv:  latst  gae 

at  the  burning  surface,  m/s 

x  Axial  coordinate,  a 

Greek  Letters 

u  Thermal  dlffuslvtty,  ra^-'a 

1  Ratio  of  specific  r.eata 

5  Gup  width  of  the  crack,  tr. 

u  <3aa  viscosity,  Rx/a»» 

M*rf  w  *  Mt»  k</ta-3 

o  Density  (without  subscript, 

gas  density),  kg/m3 

tu  Shear  stress  or.  the  port  wall, 

w  Pa 

e  Angle  measured,  In  a  counter¬ 

clockwise  direct l .n,  at  the 
lower  side  of  the  propellant, 
degree 

Subscripts 


c  Rocket  chamber 

erf  Effective 

£  Gas 

1  Initial  value 

lgn  Ignition  condition 

P  Propellant  (cond-.-n3.ci  phase) 

ps  Propellant  surface 

t  Turbulent 


Superscripts 

Favre  or  mass -averaged  nuinllty 
Reynolds  or  t  litie-n  vt-rure  i  •:<«an- 
t  Ity 

Flutuatlng  quant  It v  !n  Favre 
averaging  method 
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APPENDIX  D  (Ref.  No.  6) 

Flame  Propagation  and  Combustion  Processes 
in  Solid  Propellant  Cracks 
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NoaaencMture 

«  -pre-exponential  factor  in  Saint  Robert's  burning 
rate  relationship  r/** ,  (mm/s)/(atm)" 

An  *  Andreev  number, 

4k  «  hydraulic  diameter  of  crack,  mm 

L  • length  of  crack,  nun 

n  ■‘pressure  exponent  in  Saint  Robert's  burning  rate 
relationship 
P  •pressure,  atm 

P^  •maximum  pressure  in  the  crack  cavity,  mm 
rv  »  burning  rate  of  solid  propellant,  mm/s 
r  •  temperature,  K 

r  •time,  a 

T.  •  adiabatic  flame  temperature  of  solid  propellant,  K 
«  initial  propellant  temperature,  295  K 
vm  -  convective  ignition  front  propagation  velocity,  m/s 
x  -axial  location,  measured  from  entrance  of  crack, 
mm 

a  •  thermal  diffusfvity 

<  -gap  width  of  crack 

a  -density,  kg/m1 


Introduction 

NDBR  certain  operating  conditions,  higk  pressure  and 
high-temperature  gases  horn  the  rocket  chamber 
penetrate  the  defects  or  cracks  that  are  present  in  solid 
propellants.  The  convective  heat  transfer  from  the  hot  gases 
heats  the  propellants  along  tha  cracks,  and  if  a  sufficient 
amount  of  energy  is  transferred  to  the  propellants,  ignition 
also  occurs  inside  the  defects.  Burning  inside  defects  or 
cavities  may  result  in  much  higher  pressures  in  the  cavities; 
indeed,  it  is  believed  that  given  the  proper  stimuli,  them 
defects  or  multiple  cracks  may  initiate  detonation.  Even  if  this 
abnormal  burning  does  not  lend  to  a  catastrophic  failure,  it 
may  cause  the  performance  of  the  rocket  motor  to  deviate 
significantly  from  the  dtsigaed  conditions. 
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This  paper  deals  with  a  detailed  experimental  investigation 
of  the  development  of  convective  burning  in  an  Isolated 
propellant  crack.  Convective  burning  is  defined  ns  enhanced 
burning  of  a  propellant  charge,  caused  by  hot  gaa  penetration 
into  propellant  defects.  Forced  convection  is  the  primary  heat 
transfer  mechanism  during  convective  burning.  Even  though 
the  ignition,  flame  spreading,  and  combustion  processes 
inside  an  Irrlpfl  solid  pr*ipttlinl  crock  represent  an  idaaliaed 
cam  of  tha  actual  bunting  ia  propellant  defects,  a  fun¬ 
damental  understanding  of  the  ignition  and  flame- 
propagation  processes  in  such  a  geometry  is  essential  in  order 
to  study  the  extremely  complex  phenomenon  of  ignition  in 
defects,  lending  to  daflagration-to-detonatioa  transition 
(DDT).  It  should  he  noted,  however,  that  a  single  crock  will 
result  in  n  limited  increase  in  burning  surface  area,  and, 
therefore,  may  not  cause  detonation. 

Prior  to  this  study  several  researchers  have  investigated 
convective  burning  in  solid  propellant  defects.  Recently, 
Bradley  and  Boggs1  have  made  an  extensive  literature  review 
on  convective  burning  in  propellant  defects.  Belyaev  «t  si.1 
have  compiled  most  of  the  recent  work  done  in  this  area  in  the 
USSR.  Them  reviews  indicate  that  tha  convective  burning  in 
solid  propellant  crocks  has  been  investigated  by  a  number  of 
researchers  in  the  past  and  that  both  theoretical  and  ex¬ 
perimental  techniques  have  been  used.  Studies  in  this  area 
have  also  been  conducted  at  The  Pennsylvania  State 
University. **  Experimental  studies  on  convective  burning  in 
solid  propellant  crocks  can  be  subdivided  into  two  broad 
categories1-1:  1)  onset  of  convective  burning,  and  2) 
development  of  convective  burning.  Literature  on  tha  onset  of 
convective  burning  is  extensive;  the  studies  are  usually 
qualitative  in  nature  and  are  characterised  by  tests  of  the 
go/no-go  type.  Results  of  studies  on  the  onset  of  combustion 
are  commonly  presented  in  graphical  form  to  indicate  the 
operating  conditions  (usually  pressure)  under  which  the  (lame 
will  penatrott  a  crock  of  a  given  dimension.  The  following 
paragraphs  list  tome  of  the  important  work  in  (hisarsa. 

Prentice*  conducted  experimental  investigations  on  tha 
onset  of  convective  bunting  ia  a  transparent,  nitrocellulose- 
petrin  propellant.  He  observed  that  it  was  much  easier  to 
flashdown  into  a  crack  with  both  ends  open  than  into  &  crack 
with  one  end  doted.  Catalytic  additives  caused  flashdown  tc 
occur  more  readily.  Flashdown  was  not  observed  up  to  3.9 
MPa  in  i.d-amv-diam  cracks.  Both  motion  pictures  and 
the.-mocouplei  were  used  to  monitor  flame  penetration. 
Later,  Prentice*  extended  Mi  work  to  opaque  composite 
propellants.  Vibration  response  spectroscopy  ( VRS)  was  used 
to  monitor  flashdown.  Double  based  propellants  that 
dispuyed  mau  burning  were  also  investigated  by  Prentice.  He 
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concluded  that  the  critical  pore  dlamttar  for  nathdown  In  a 
given  propellant  U  inversely  proportional  to  the  burning  rata. 
Margolin  and  MargulisM.  noted  that  com  butt  Ion  panttratat 
more  easily  into  a  channel  with  two  open  ends  than  Into  a 
blind  channel,  they  alto  observed  that  if  the  channel  were 
open  at  one  end.  or  if  there  were  one  or  more  outlets  front  the 
channel  to  the  burning  surface,  the  combustion  gates  in* 
stantaneously  penetrated  the  pores  when  the  Andreev  number 
(Ait)  was  beyond  a  critical  value.  However,  the  critical  An* 
dreev  number  is  not  always  constant. 

Oodai11  made  experimental  investigations  of  name 
propagation  in  the  narrow  slit  and  fine  hole  of  solid 
propellant  grain.  Cracks  of  22  mm  in  length  were  formed  with 
S  mm  square  slabs  of  AP*based  composite  propellents.  The 
propellant  was  Ignited  with  a  nkhrome  wire,  and  the  ignition 
event  was  filmed  by  motion-picture  camera  at  a  rate  of  about 
32  frames/s.  Tests  were  conducted  in  a  slab  of  propellant  with 
both  a  single  hole  and  multiple  holes.  Oodai  observed  that 
there  is  a  threshold  gap  width  or  hole  diameter  below  which 
the  flame  does  not  propagate  into  the  defect,  and  that  the 
threshold  crack  hap  is  fundamentally  a  function  of  the 
burning  rate.  Bobolev  ct  a!., |J  Payne, 11  and  Krasnov  ct  al. u 
also  studied  the  onset  of  convective  burning  in  solid 
propellant  cracks. 

Experimental  studies  on  the  development  of  convective 
burning  are  limited.  Belyaev  et  al.l,  u  made  preliminary 
investigations  of  the  development  of  combustion  in  single 
pores.  The  flame  propagation  rate  was  observed  to  increase 
initially,  and  then  to  reach  a  constant  value.  The  pressure  in 
the  cr'jca  was  found  to  be  greater  than  that  of  the  chamber, 
and  tlte  effect  of  erosive  burning  was  found  to  be  important. 
The  maximum  pressure  in  the  crock  was  correlated  with  the 
crack  lengt  i  and  the  ratio  of  the  crack  length  to  crack  width. 

The  focus  of  previous  investigations  of  combustion  in 
propellant  cracks  has  been  on  the  onset  of  convective  burning; 
the  development  of  convective  burning  has  yet  to  be  in¬ 
vestigated  fully.  It  is  generally  believed  that  the  development 
of  convective  burning  in  a  propellant  defect  is  a  necessary 
condition  for  subsequent  transition  to  detonation.  Almost  ail 
of  the  previous  studies  have  been  conducted  at  constant 
pressure  conditions.  The  effect  of  chamber  pressurization  rate 
(dP/dri  on  (lame  propagation,  which  is  of  interest  in  actual 
rocket  motor  and  DDT  studies,  has  not  yet  been  investigated. 

The  aim  of  the  present  research  was  to  narrow  some  of  the 
technological  gaps  that  exist  in  the  study  of  the  development 
of  convective  burning  in  solid  propellant  cracks.  The  specific 
objectives  of  this  investigation  are: 

U  To  experimentally  study  the  effect  of  the  following 
parameters  on  ignition  and  flame  propagation  inside  Am- 
monium-Perchlorate-bascd  composite  solid  propellant 
cracks:  a)  chamber  pressurization  rate.  b)  cractaap  width,  c) 
length  of  the  crack,  and  d)  propellant  compositiSr . 

2)  To  establish  a  data  base  for  the  convective  flame- 
propagation  rate  and  the  maximum  pressure  in  the  crack 
cavity  as  a  function  of  pressurixation  rate,  crack  geometry, 
and  the  physicochemical  pioperties  of  the  propellant. 

The  results  of  this  study  can  alro  be  used  for  verification 
purposes,  in  order  to  compare  till  flame-propagation  rates 
and  pressure  measurements  along  axial  locations  of  the  crack, 
with  the  numerical  predictions  of  a  parallel  theoretical  study 
being  conducted  by  the  authors  at  The  Pennsylvania  State 
University. 

Experimental  Setup  and  Procedures 

Test  Asperates 

The  test  apparatus  and  experimental  conditions  for  this 
study  incorporated  the  following  major  design  con¬ 
siderations:  1)  the  test  apparatus  should  allow  a  wide  range  of 
variation  in  crack  geometry  and  operating  conditions:  and  2) 
test  conditions  should  simulate,  as  closely  as  possible,  those  in 
an  actual  rocket  motor.  A  solid  propellant  igniter  system 
(details  of  which  are  given  in  the  following  section)  was  used 


to  simulate  igr  tar  gas  composition  similar  to  that  in  an  actual 
rocket  motor,  c  ing  this  igniter  system,  it  was  possible  to 
obtain  pressures  up  to  200  atm  in  the  chamber  in  ap¬ 
proximately  2  ms.  The  test  apparatus  is  also  compatible  with 
the  theoretical  model  developed  by  the  authors and. 
therefore,  the  results  obtained  here  could  later  be  used  for 
model  validation. 

/g Hirer  System 

A  solid  propellant  igniter  system  was  the  source  for  hot  gas 
generation.  Percussion  primers  were  used  as  initial  on  for 
ignition  of  a  propellant  charge.  A  remotely  controlled, 
solenoid-activated  firing  pin  was  used  to  trigger  the  primer.  A 
schematic  diagram  of  the  test  chamber  and  igniter  system  is 
shown  in  Fig.  1.  The  igniter  system  consists  of  1)  a  solenoid 
mounted  on  the  retainer  block  of  a  spring-loaded  firing  pin,  2) 
a  percussion  dement  and  its  housing  unit,  31  a  propellant 
igniter  charge  and  igniter  chamber,  and  4)  a  multipcrforated 
nozzle  plate.  Ail  hardware  was  made  of  stainless  steel. 

High-temperature,  high-velocity  gases,  under  rapidly  in¬ 
creasing  pressure  conditions  that  simulate  flow  conditions  in  a 
rocket  chamber,  were  obtained  by  burning  a  50  mm  long, 
rectangular  cross-sectional  ( ~  S.7  x  S.7  mm)  solid  propellant 
chart*.  This  charge  was  placed  inside  the  cavity  of  the  igniter 
chamber,  the  diameter  of  the  cavity  is  S  mm.  When  the  primer 
was  triggered,  it  produced  hot  gases  which  flowed  over  and 
ignited  the  solid  propellant  charge.  The  product  gates  flowed 
through  a  muliiperforated  converging  nozzle  into  the  main 
chamber.  The  pressurization  rate  of  the  chamber  can  be 
varied  by  altering  the  dimensions  of  the  propellant  strip  in  the 
igniter,  by  changing  the  dimensions  of  tlte  multiperforoted 
nozzle,  or  by  changing  the  exit  nozzle  of  the  chamber. 

Te a  CMtmttf  end  Of*  Con/igttiwmMt 

The  test  chamber  shown  in  Fig.  I  was  made  of  stainless 
steel  304.  The  discharge  of  the  igniter  system  passes  through  a 
124  mm  long,  rectangular  cross-sectional  (10x25.4  mm) 
channel.  The  propellant  crack  sample  is  placed  perpendicular 


to  the  flow  direction  in  the  main  chamber.  Cracks  were 
formed  by  cutting  a  slot  of  desired  width  in  a  propellant  slab 
(of  25.4  x  17  inn  cross  section)  glued  into  a  brass  retainer,  it 
‘  as  possible  to  manufacture  cracks  up  to  210  mm  in  length, 
with  gap  widths  as  low  as  450  am.  Figure  2  shows  a  typical 
propellant  crack  sample.  Four  pressure  ports  were  provided 
along  the  center  of  the  crack  in  order  to  measure  pressure  • 
variation  along  the  crack.  An  interchangeable,  convergent 
exit  nozzle  was  used  to  vary  pressure  and  flow  conditions  in 
the  chamber.  The  nozzle  was  made  of  titanium  in  order  to 
alleviate  the  problem  of  metal  erosion  caused  by  highly 
corrosive  gases.  The  exit  diameter  of  the  nozzle  was  varied 
between  2.5  and  3.8  mm.  For  safety,  the  chamber  was  also 
equipped  with  a  port  for  a  burst  diaphragm.  If  for  any  reason 
pressure  inside  the  chamber  exceeded  a  predetermined  value 
(usually  680  atm),  the  burst  diaphragm  would  rupture  to  vent 
pressure  and  avoid  an  explosion. 

Flame  propagation  in  the  crack  was  observed  through  a  set 
of  transparent  plexiglass  windows.  The  window  assembly 
consisted  of  a  sacrificial  window  (238  x  65  x  6.3  mm)  and  a 
viewing  window  (251  x  78  x  25.4  mm).  The  window  assembly 
was  held  in  place  by  a  stainless  steel  window  retainer.  To 
achieve  a  good  seal,  rubber  O-rings  were  used  between  the  two 
halves  of  the  chamber,  and  between  the  inlet  and  exit  nozzles 
and  the  chamber-  During  the  tests,  the  chamber  was  com¬ 
pletely  sealed,  except  for  the  interchangeable  exit  nozzle 
through  which  product  gases  were  discharged  into  the  at¬ 
mosphere.  Since  the  sacrificial  window  is  placed  directly  over 
the  propellant  surface,  it  is  destroyed  after  each  test  firing. 
The  test  sample  is  clearly  visible  through  the  windows  during 
the  test. 

To  obtain  a  more  detailed  observation  of  the  flame-front 
propagation,  an  alternative  crack  sample  shown  in  Fig.  3  was 
used  to  replace  the  one  shown  in  Fig.  2.  The  alternative  crack 
sample  provided  direct  (front  view)  observation  of  the  gas 
penetration  and  ignition-front  propagation  processes.  In  this 
configuration,  the  crack  was  formed  between  a  propellant 
slab  and  the  sacrificial  plexiglass  window,  i.e.,  one  side  of  the 
crack  was  an  inert,  transparent,  plexiglass  window  and  the 
other  side  was  a  propellant  slab  glued  to  a  stainless  steel  base 
plate.  The  gap  width  of  such  a  crack  configuration  is  varied 
.  by  the  amount  by  which  the  propellant  surface  is  recessed 
below  the  side-leg  assembly  (see  Fig.  3).  The  propellant  slabs 
for  these  tests  were  183  mm  long  and  (7.7  mm  wide.  Gap 
width:  of  the  crack  were  varied  between ,0.43  and  1.5  mm .  It 
should  be  noted  Jhat  for  this  configuration  only  two  pressure 
traces  can  be  obtained:  one  at  the  crack  entrance,  and  the 
dther  at  the  tip. 

laalnHBMlatioa  . 

Data  Acquisition  System 

A  block  diagram  of  the  data  acquisition  system  for  the 
convective  flame  propagation  studies  is  shown  in  Fig.  4,  The 
data  acquisition  system  consists  of  three  major  parts:  I) 

ressure  measuring  system,  J)  flame  propagation  measuring 
sySttiti,  and  3)  transient  waveform  recording  system.  Details 
of  measurement  techniques  are  given  in  the  following  suc¬ 
tions.  The  data  acquisition  system  includes:  pressure  trans¬ 
ducer,  charge  amplifier:,  a  high-speed  movit  camera,  a  light- 
emitting  diode  (LED)  driving  unit,  a  motion  analyzer,  a  9- 
chanr.el  transient  waveform  recorder  (Physical  Data  Model 
#515.234),  a  ^-channel  PM  tape  recorder,  an  oscilloscope  and 
an  x-y  plotter.  The  maximum  sampling  rate  of  the  Physical 
Data  System  is  2x  10‘  sample:/  ),  and  ‘he  maximum  am¬ 
plitude  resolution  is  p.lV#  with  a.  4096  Word  memory  per 
channel.  The  recording  speed  of  the  HP  type  recorder  is  60 
i(r./s  at  a  carrier  frequency  qf  :,KA  kHz.  The  iape  can  be 
played  back  "gt  a  much  lower  speed  (1 H  In  /s)  to  expand  the 
duration  of  the  transient  signals.  Data  recorded  on  the 
Physical  Data  Systim  can  be  displayed  on  an  oscll’oscoye,  or 
plotted  on  an  x-.v  plotter  for  a  hard  copy.  Data  recorded  on 
tape  is  digitized  before  displey  or  plotting .  • 


Fig.  3  PrwgeHant  sample  ior  detailed  obamatioo  of  Igetlhm  front. 


Fig.  4  Block  diagram  of  data  acqulaltioa  a?*tem. 


Pressure  Measurements 

Four  ports  are  provided  along  the  length  of  the  crack  for 
transient  pressure  measurements.  The  first  port  is  located  in 
the  main  chamber,  very  close  to  the  crack  entrance.  The  other 
three  are  located  at  48,  138,  and  188  mm  from  the  first 
transducer  port.  Piezoelectric  quartz  transducers  were  used  to 
measure  the  pressure.  The  transducers  have  a  rise  time  of  1.5 
jis  and  a  natural  frequency  of  300  kHz,  and  can  accurately 
record  pressures  up  to  1000  atm.  The  transducers  were 
mounted  in  a  water-cooled  adapter  that  prevented  drifting  or 
damage  due  to  excessive  heat.  In  addition,  a  thin  layer  of 
silicon  rubber  insulation  was  placed  on  the  transducer  surface 
to  further  protect  the  transducer  diaphragm  from  high- 
temperature  gases.  The  transducers  were  mounted  about  2 
mm  below  the  bottom  surface  of  the  chamber. 

Charge  signals  from  the  pressure  transducers  were  carried 
through  an  insulated,  high-impedance,  coaxial  cable  to  a 
charge,  amplifier,  The  output  (voltage  signal)  of  the  charge 
amplifiers  was  recorded  simultaneously  on  the  transient 
waveform  recorder  and  on  the  magnetic  tape.  In  order  to 
accurately  determine  the  pressure  level,  calibration  signals 
equivalent  to  69  atm  (1000  psig)  were  also  recorded  for  each 
channel.  Calibration  signals  for  each  transducer  were  checked 
periodically  by  comparing  them  with  signals  obtained  from 
measuring  a  known  pressure  of  a  high-pressure  nitrogen  or  air 
tank. 
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Ftumt-PropataUon  Mtasurtmtius 

A  high-speed,  16  mm  motion  picture  camera  was  used  to 
observe  the  complete  igffitibh  event  in  (he  chamber.  With  a 
quarter-frame  optical  head,  site  camera  was  capable  of 
filming  at  a  maximum  of  about  *4,000  pictures/s.  The  Hycam 
camera  is  equipped  with  a  dual  light-emitting  diode  (LED) 
system;  one  LED  records  timing  marks  and  the  other  records 
a  common-time  signal  on  the  film  The  common-time  signal  is 
generated  at  the  start  of  the  event;  and  is  used  to  correlate 
events  recorded  on  the  transient  waveform  recorder  or  tape 
with  those  on  the  film.  The  LED  driver  unit,  a  common-time 
pulse  and  high-frequency  (Up  to  10  kHz)  timing  generator, 
was  used  to  operate  the  LED  time  system.  The  frequency  of 
the  timing  pulse  in  the  present  investigation  was  set  at  l  or  10 
kHz. 

The  Vanguard  motion  analyzer  was  used  to  analyze  the 
film,  frame  by  frame.  The  motion  analyzer  screen  is  equipped 
with  two  cross  hairs  to  (raver::  in  two  perpendicular  direc¬ 
tions  in  a  plane.  Two  micrometer  dials,  accurate  to  a 
thousandth  of  an  inch,  are  used  to  record  the  movement  of  the 
cross  hairs.  The  analyzer’s  frame  counter  is  designed  to  keep 
track  of  the  number  of  frames  analyzed.  In  order  to  obtain 
the  flame-propagation  speed  for  a  fixed  y  location  (centerline 
of  the  crack),  x  readings  (along  the  crack)  were  taken  for  each 
frame  until  the  flame  reached  the  crack  tip.  Readings  of  the 
micrometer  dials  were  scaled  by  using  the  known  value  of  the 
initial  length  of  the  crack  and  the  total  distance  traversed  by 
the  crosswire.  The  time  interval  between  each  picture  was 
obtained  from  the  time  marks  on  the  film. 

Remotely  Controlled  ignition  Circuit 

The  transient  ignition  processes  in  a  solid  propellant  crack 
involve  high  pressures  and  very  high  pressurization  rates 
( - 10*  atm/s).  Because  pressurization  rates  of  such  a 
magnitude  are  potentially  hazardous,  experiments  were 
carried  out  in  a  test  cell  with  l'/j-ft  thick  concrete  walls,  and 
the  ignition  event  was  remotely  controlled. 

A  block  diagram  of  the  remotely  controlled  ignition  and 
high-speed  photography  system  is  shown  in  Fig.  S.  The 
solenoid  that  triggers  the  percussion  primer  is  activated  by  an 
event  switch  built  into  the  high-speed  camera.  After  the 
camera  is  switched  on,  about  one  second  is  required  to  reach 
the  desired  framing  rate.  This  delay  time  can  be  related  to  the 
footage  of  the  film  that  passes  through  the  camera  before  the 


Fig.  S  Block  diagram  of  mnoleiy  controlled  igaiUoo  Md 
photograpO)'  system. 


attainment  of  the  prescribed  framing  rate.  The  event  switch  is 
controlled  by  a  footage  indicator  Daring  operation,  the 
footage  indicator  is  set  to  a  preselected  valii?  corresponding  to 
the  desired  framing  rate.  For  example,  the  indicator  is  set  at 
100  ft  for  a  framing  rate  of  40,000  pictures/s.  After  a 
specified  length  of  film  (preset  on  the  footage  indicator) 
passes  through  the  camera,  the  event  switch  closes,  causing  12 
V  ac  to  be  applied  across  the  relay,  and  the  relay  is  activated. 
The  relay,  in  turn,  activates  two  other  switches;  one  switch 
triggers  the  solenoid  (by  causing  a  current  to  flow  through  it) 
and  initiates  the  ignition  event,  and  the  other  closes  the 
common-time  switch  on  a  light-emitting  diode  (LED)  driver 
unit.  At  the  instant  at  which  the  common-time  switch  closes, 
the  LED  driver  unit  generates  a  2  ms  pulse  that  is  recorded  on 
the'  transient  waveform  recorders  and  is  simultaneously 
marked  on  the  film.  These  common-time  marks  are  used  for 
time  correlation  of  the  data  recorded  on  the  tape  and  film. 

To*  Aoctfiffi 

Test  samples  were  prepared  from  cast  propellant  grains 
obtained  from  propellant  processing  laboratories.  Propellant 
slabs  (approximately  215  x  26  x  17  mm)  were  cut  from  the  cast 
propellant  with  a  remotely  controlled  milling  machine.  A 
slightly  oversized  propellant  slab  was  glued  into  a  brass 
retainer.  The  sample  was  left  to  dry  for  24  h  in  order  to 
achieve  a  good  adhesive  bond  between  the  propellant  and  the 
brass  retainer.  The  sample  was  machined  again  to  the  desired 
dimensions.  A  crack  of  the  desired  dimension  was  then 
precisely  cut  through  the  center  of  the  crack  sample  with  a 
jeweler’s  slitting  saw.  Propellant  samples  for  the  alternative 
configuration  (for  detailed  flame-front  observation)  were 
obtained  in  a  similar  manner.  Oversized  propellant  slabs  were 
glued  onto  a  steel  plate  and  machined  to  the  desired  dimen¬ 
sions.  The  entrance  region  of  the  crack  was  rounded  slightly 
to  ensure  smooth  flow  development  inside  the  crack.  After 
machining,  the  samples  were  carefully  cleaned  with  a  high- 
pressure  air  jet  to  remove  any  loose  particles  in  the  crack 
cavity. 

During  assembly  of  the  test  chamber,  a  part  of  the 
propellant  side  surface  in  contact  with  the  sacrificial 
plexiglass  window  and  the  combustion  chamber  was  coated 
with  a  thin  layer  of  (lame  retardant  to  help  prevent  the  flame 
from  penetrating  between  the  propellant  and  the  contact 
surfaces.  A  layer  of  rubber-based  adhesive  sealant  wa^  ap¬ 
plied  to  all  contact  surfaces  of  the  various  components  of  the 
lest  chamber  to  prevent  leaks  during  experiments.  After  the 
crack  sample  was  mounted  into  the  crack  combustion 
chamber,  the  crack  geometry  was  carefully  measured  and 
recorded.  High-intensity  camera  lights  were  used  during. the 
test  to  facilitate  motion-picture  recording.  In  order  to  observe 
the  ignition  event  remotely  and  to  detect  possible  leaks  during 
the  test,  a  video  monitor  with  a  TV  camera  was  used. 

R*sult>  and  Discussion 

Four  types  of  AP-based  composite  solid  prc-ellan;;  were 
used  in  this  investigation.  Table  I  lists  the  propellants  and 
some  of  their  properties.  Propellants  A  and  P  have  almost 
identical  pressure  exponents  in  the  Saint  Robert's  burning  rate 
relationship;  however,  propellant  A  has  a  higher  burning  rate 
since  it  has  a  larger  pre-exponential  factor.  Propellant  B 
(AP/PBAA-EPOM)  wsis  used  in  most  of  the  tests.  Propellants 
C  and  D  have  identical  binder  (HTPB)  and  oxidizer  weight 
fractions  (71"  );  however,  the  oxidizer  particle  size  lit 
propellant  C  is  20  ^m,  whereas  the  oxidizer  penicle  size  <n 
propellant  D  is  200  nm.  The  effect  of  or.iditcr  pat  tide  size  On 
the  ignition  and  Jlame  ptopagarion  processes  can  be  obtained 
from  a  comparison  of  the  results  of  these  twu  propellents. 

The  physical  process*",  that  take  place  during  ihe  test  can  he 
described  generally  as  follows.  Discharge  of  urodun  gate 
from  the  igniter  system  pressurizes  the  test  chamber,  causing 
hot  ignition  gases  to  penetrate  the  crack  cavity.  TV,  prpvure- 
wave  fronts  move  along  the  crack,  and  the  hot  arses  transfer 
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Table  I  Proprilaai  proper* k* 


Propellant  type  A  B 

Composition  AP-bascd  AP/PBAA-EPON 

Weight  percent  of  oxidizer  •  •  *  t>: 

Average  particle  sire  d ,  jun  •  •  • 4  7S 

Pre-exponential  factor  tn  Saint  Roiieri's  1.62  0.5591 

burning  rate  expression  a,  mm/s/faimi* 

Pressure  exponent  in  Saint  Robert's  0.4106  0.41 

burning  rale  expression,  n 

Flame  temperature  7).  K  3000  1920 

Propellant  density,  kg/m*  1710  1000 

*AP  weight  percent  and  particle  site  for  propellant  A  not  availaMe  to  the  authors. 
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AP/HTPB 

72 

20 

0.8441 


AP/HTPB 

72 

200, 

0  5849 


heat  to  (he  crack  surface,  resulting  in  an  increase  in  propellant 
surface  temperature.  As  the  pressure  in  the  chamber  continues 
to  increase,  more  hot  gases  are  driven  into  the  crack,  and 
there  is  also  an  increase  in  the  convective  heat-transfer  rate. 
Eventually,  the  ignition  condition  is  realized,  and  the 
propellant  begins  to  burn.  Following  the  flame-spreading 
period,  the  gasification  of  the  crack  surface  may  generate  a 
pressure  higher  than  that  in  the  chamber.  Finally,  the  gases 
flow  out  of  the  crack  into  the  main  chamber. 

As  mentioned  in  the  section  on  instrumentation,  pressure 
measurements  were  made  at  four  axial  locations  along  the 
crack.  The  ignition-front  propagation  speed  was  deduced 
from  films  of  the  flame-front  propagation  made  by  a  high¬ 
speed  camera.  Because  of  the  complexity  involved  with  ac¬ 
curate  quantitative  measurement  of  temperature  for  highly 
transient  flow  situations,  no  attempt  was  made  to  measure  the 
temperature  in  the  crack.  (The  time  required  for  the  flame  to 
propagate  from  the  crack  entrance  to  the  crack  tip  was  less 
than  1  ms  in  most  cases.)  High  gas  temperatures  also  make 
this  kind  of  measurement  difficult.  In  all  test  firings  the  initial 
pressure  in  the  crack  chamber  was  1  atm,  and  the  initial 
temperature  of  the  propellant  was  about  295  K.  Results  of 
pressure  and  flame-propagation  measurements  are  discussed 
next,  followed  by  results  of  the  parametric  study. 

A  typical  set  of  time-correlated  pressure  traces  is  shown  in 
Fig.  6.  In  this  figure  a  dual  time  base  is  used  to  obtain 
maximum  information  in  one  plot.  The  initial  portions  of  (he 
P-t  traces  have  been  expanded  to  i'lushate  the  pressurization 
processes  at  various  locations  along  the  ctack.  Gage  G2  is 
located  at  the  crack  entrance,  and  G5  near  the  crack  tip.  The 
qualitative  nature  of  the  curves  is  similar,  with,  the  following 
important  differences:  1)  the  first  discernible  pressure  rise  for 
each  pressure  gage  occurred  consecutively  from  the  crack 
entrant  to  the  cruclc  tip,  i.e.,  from  02  t.V  GJ;  2)  the 
pressurization  rate  increased  consecutively  downstream  from 
the  crack  entrance;  and  3)  the  maximum  pressure  occurred  in 
the  interior  of  the  crack.  1  he  time  delay  between  the  first 
discernible  pressure  rise  at  downstream  locations  and  that  at 
the  crack  entrance  Is  caused  by  the  finite  time  required  for  the 
hoi  product  gasw  to  travel  from  ths  entrance  to  the  do  dating 
.  initial  pressurization,  the  increase  in  pressurization  rate  at 
do  .(stream  locations  is  caused  by  th?  coalescence  of  the 
(ravelins  pressure  waves.  This  can  be  further  explained  by 
noting  that  at  the  igniter  is  discharged,  the  pressure  of  the 
chamber  rises,  causing  vyzak  pressure  waves  ’o  travel 
downstream  of  the  crack.  A*  the  pressure  waves  travel  down 
the  crack,  they  are  followed  by  stronger  pressure  waves 
carted  by  rabidly  ^creasing  pressute  in  ti,e  chamber.  A*  the 
lame  tin:?,  s^rne  of  the  propellant  in  the  upstream  region 
begins  to  gisify,  thus  creating  higher  pressure  behind  the 
pressure  front.  The  result  of  the  cote.bined  effect  of  the  high 
pressure  regions  behind  the  /rout  it  a  steepening  of  the 
prwiurc  front  us  it  moves  downstream. 

1 ■  The  pressure  true  (it  any  location  can  be  subdivided  into 
three  regions:  I!  the  initial  rapid  pressurization  region,  which 
is  Whfolled  mainly  by  chamber  conditions;  ?.)  the  slowly 
rising  pressure  region;  anil  3)  the  depressurization  region.  The 
initial  uprising  portion  of  the  pressure-time  tiace  is  quite 


linear,  in  most  cases,  the  flame  reaches  the  crack  tip  during 
the  initial  uprising  portloit  of  pressure-time  trace  at  the  tip. 
Pressure  in  regions  2  and  3  is  controlled  by  combustion  and 
flow  processes  in  the  crack.  As  the  burning  of  the  propellant 
in  the  crack  continues,  the  pressure  in  the  cavity  rises  above 
that  at  the  crack  opening;  later,  this  pressure  difference  causes 
the  gases  to  flow  out  of  the  crack.  Pressure  in  the  crack 
continues  to  rise  until  the  mass  flow  rate  of  the  gases  from  the 
crack  is  greater  than  that  generated  by  the  burning  of  the 
propellant  in  the  cavity.  Since  the  pressure  in  this  region  is 
governed  by  combustion  processes  inside  the  crack,  and 
because  gage  G2  is  located  at  the  crack  opening  in  the  main 
chamber,  the  highest  peak  pressure  occurs  in  the  interior  of 
the  crack.  Depressurization  is  caused  by  the  increase  in  crack 
cavity  volume,  as  well  as  by  the  enlargement  of  the  crack 
opening,  but  the  burning  surface  area  remains  constant. 

The  common-time  signal,  which  is  used  to  correlate 
pressure  measurement  and  film,  is  not  shown  in  Fig.  6  because 
of  the  10-12  ms  time  lag  between  the  common-time  signal  and 
first  discernible  pressure  rise  at  the  crack  entrance  location.  In 
order  to  include  the  common-time  signal,  pressure-time  traces 
would  have  to  be  compressed  substantially,  which  is  not 
desirable.  The  time-correlated  pressure  and  ignition  data 
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\  Mon -d taxational  Axial  DUcanct,  x/t 

Fig.  I  Uut-»gaurn  polyaoaiil  fit  (o  KMint  IgaiUoa  (raw 
locatkM  vt  llm*  da**. 

indicate  thin  the  appearance  of  hot  igniter  gates  at  the  crack 
entrance  coincides  with  t^e  first  discernible  pressure  rise  at 
that  point.  Therefore,  the  initial  time  i  »0  in  this  study  can  he 
defined  as  either  the  first  discernible  pressure  rise  at  the  crack 
entrance  location,  or  the  first  appearance  of  hot  gases  on  the 
film. 

Figure  7  shows  typical  continuous  photographs  of  the 
ignition-front  propagation  in  a  crack.  The  instantaneous 
location  of  the  ignition  front  x  was  obtained  by  analyzing  the 
film  of  the  ignition  event  on  a  motion-picture  analyzer,  as 
discussed  in  the  last  section.  The  time  marks  on  the  film  were 
used  to  determine  the  time  t  corresponding  to  each  reading.  A 
least-  ;quares  polynomial  was  fitted  to  the  x-r  data,  and  an 
analytic  derivative  of  this  curve  was  used  to  obtain  the 
ignition-front  propagation  rate  «dx/d /) .  This  procedure 
made  it  possible  to  obtain  flame-propagation  rates  from  the 
measured  ignition-front  locations.  Polynomials  between  the 
third  and  fifth  order  were  used,  and  the  polynomial  with  the 
best  fit  was  chosen  to  represent  the  analytic  function  relating 
the  instantaneous  flame- front  location  to  time.  Figure  I 
shows  a  typical  fourth-order  polynomial  fit  to  the  measured 
ignition-front  location  vs  time  data.  As  can  be  seen  from  this 
figure,  the  curve  fits  the  data  quite  well.  The  sum  of  the 
square  of  errors  for  this  polynomial  fit  is  4. 143  x  10 ' J .  In  this 
paper,  she  remainder  of  the  figures  of  ignition-front  location 
vs  axial  distance  show  only  the  fitted  polynomial.  Errors  are 
of  the  same  order  at  that  given  in  Fig.  i. 
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The  effect  of  chamber-pressurization  rate  on  location  and 
velocity  of  the  ignition  front  for  propellant  B  is  shown  in  Fig. 
9.  The  crack  specimen  was  195  mm  long,  and  the  gap  width 
was  0.89  mm.  The  time  for  the  flame  to  propagate  to  a  given 
axial  location  and  the  velocity  of  the  ignition  front  at  that 
location  are  plotted  as  a  function  of  the  nondimcnsional  axial 
coordinate  fot  pressurization  rates  of  i.i  x  10’,  2.6x10’, 
and  4.2  x  10’  atm/s.  The  time  required  for  the  flame  to  reach 
the  crack  tip  is  836,  511*  and  391  **s,  respectively.  That  is,  as 
the  chamber  prcuurizudon  rata  increases,  the  ignition  front 
propagates  much  more  quickly  into  the  crack  because  the 
rapidly  increasing  pressure  in  the  chamber  acts  as  the  driving 
force  for  the  penetration  of  hot  gases  into  the  crack.  This 
higher  pressurization  causes  an  increass  in  the  velocity  of  the 
gas  flowing  into  the  crack,  and  results  in  a  higher  rate  of  heat 
transfer  to  the  crack  wails.  Faster  ignition-front  propagation 
is  a  net  effect  of  these  conditions.  It  should  be  noted  that 
because  of  the  presence  of  bright  gases  near  the  crack  en¬ 
trance,  it  is  extremely  difficult  to  identify  the  time  at  which 
the  ignition  front  at  the  crack  opening  first  appears.  This 
explains  why  the  plots  do  not  show  data  near  x-0. 

It  is  interesting  to  note  that  the  nature  of  the  curves  is  the 
same  in  all  three  cases.  By  extrapolating  the  time  vs  axial 
location  curve  to  x  •  0,  the  time  lag  betwee  *  the  arrival  of  hot 
gases  in  the  chamber  and  the  initiation  of  the  ignition  front 
near  the  crack  entrance  is  evident.  The  time  tag  decreases  as 
(he  pressurization  rate  is  increased.  The  ignition-front  velocity 
vs  the  nondimensional  axial  coordinate  curve  shows  that  after 
the  ignition  front  is  created  near  the  crack  entrance,  it  ac¬ 
celerates,  reaches  a  maximum,  and  then  decelerates.  This 
deceleration  near  the  closed  end  is  more  pronounced  with  a 
higher  pressurization  rate.  At  a  low  pressurization  rate,  the 
velocity  remains  relatively  uniform  over  the  bulk  of  the  crack 
length,  initial  acceleration  is  caused  by  the  preheating  by  the 
hot  gases  that  precede  the  itpiilon  from.  Deceleration  is 
believed  to  be  caused  by  the  end  effects.  Deceleration  of  the 
ignition  front  is  in  agreement  with  the  observation  made  by 
Taylor"  concerning  ignition  of  granular  propellants  with  a 
closed-end  boundary.  Maximum  ignition-front  propagation 
velocities  are  410,  673,  and  886  m/s  for  the  pressurization 
rates  of  i.i  x  10*,  2.6x10’,  and  4.2x10’  atm/s,  res  pec- 
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tively.  it  is  important  to  note  that  the  effect  of  d/Vdf,  as 
depicted  in  Fig.  9,  is  valid  only  for  the  range  of  values  in¬ 
dicated.  As  observed  earlier  by  the  authors,4  under  very  high 
d/Vdf  conditions  the  crack  tip  may  ignite  before  the  ignition 
front  propagates  to  the  tip. 

The  effect  of  gap  width  on  instantaneous  location  and 
velocity  of  the  ignition  front  for  propellant  B  is  shown  in  Fig. 
10.  The  length  of  the  crack  is  195  mm,  and  the  crack-gap 
widths  are  0.308, 0.889,  and  1 .27  mm.  As  the  crack-gap  width 
is  increased,  the  ignition  front  penetrates  the  crack  more 
quickly  and  reaches  the  tip  sooner  because  the  larger  gap 
width  offers  less  resistance  to  the  penetration  of  the  flow  into 
the  cavity.  However,  maximum  velocity  increases  as  the  gap 
width  decreases.  This  may  be  explained  as  follows.  Since  the 
time  delay  required  for  the  ignition  front  to  establish  itself 
near  the  crack  entrance  is  longer  for  cracks  with  smaller  gap 
widths,  the  preheating  effect  caused  by  convective  heat 
transfer  is  more  pronounced.  Once  gasification  begins, 
pressure  behind  the  ignition  front  (due  to  gasification)  will 
also  be  higher  in  narrow  cracks.  It  is  believed  that  the  com¬ 
bined  effec'  of  preheating  and  higher  pressure  causes  the 
ignition  front  to  accelerate  relatively  rapidly  in  narrow  cracks. 

Figure  11  shows  the  effect  of  gap  width  on  maximum 
pressure  within  the  crack  cavity.  It  is  dear  that  as  the  gap 
width  decreases,  the  maximum  pressure  in  'he  cavity  in¬ 
creases.  As  discussed  earlier,  maximum  pressure  in  the  crack 
cavity  is  determined  uy  ignition  and  combustion  processes 
inside  the  crack.  Since  crack  depths  are  identical  for  all 
samples  (23.4  mm),  burning  surface  areas  will  also  be 
identical  as  long  at  crack  lengths  for  all  samples  are  the  same. 
As  gap  width  is  decreased,  the  burning  surface  area  remains 
constant,  but  the  opening  area  of  the  crack  and  volume  of  the 
crack  gap  decrease;  therefore,  maximum  pressure  is  higher  for 
smaller  gap  widths. 

The  sffect  of  the  burning  rate  on  ignition-front  propagation 
was  studied  by  comparing  results  of  tests  conducted  on 
propellants  C  and  D  as  shown  in  Fig.  12.  The  initial  crack 
length  and  gap  width  were  193  and  0.89  mm,  respectively;  the 
Initial  chamber  prwsurixation  rate  was  3.2  x  10 5  atm/s.  This 
comparison  shows  that  the  ignition-front  propagation  rate  is 
faster  for  propellant  C  with  its  higher  burning  rate  (AP 
particle  tire  -20  Mm),  than  it  is  for  propellant  D  with  its 
lower  burning  rate  (AP  particle  sue  -200  Mm).  In  these  tests 
the  effect  of  surface  roughness  was  somewhat  suppressed 
during  the  ignition-front  propagation  period  because  of  the 
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fact  that  surfaces  were  machined.  The  rate  of  llame 
propagation  for  a  propellant  with  a  higher  burning  rate  is 
faster  because  the  high  burning  rate  produces  more  gases 
behind  the  ignition  front  and  causes  the  local  pressure 
gradient  near  the  front  to  increase. 

The  effect  of  propellant  type  was  investigated  further  by 
comparing  the  results  of  high-energy  propellant  A  with  those 
of  propellant  B.  Even  though  both  are  AP-based  propellants, 
the  flame  temperature  and  burning  rate  of  propellant  A  are 
much  higher  than  those  of  propellant  B.  Figure  13  shows  a 
comparison  of  the  ignition-front  propagation  rates  of  the  two 
propellants.  The  crack  samples  were  193  mm  long,  with  a  gap 
width  of  0.89  m.  i.  The  Initial  chamber  pressuriiation  rate  war 
3.2  x10s  atm/s.  The  more  energetic  propellant  A  has  a 
shorter  flrme-spreading  period  and  a  higher  maximum 
ignition-front  velocity.  Besides  having  a  higher  burning  rate, 
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the  product  guts  generated  from  propellant  A  also  have  a 
higher  temperature  which  further  facilitate*  the  ignition-front 
propagation  procesa.  Table  2  shows  the  effect  of  propellant 
burning  rate  on  maximum  pressure  in  the  cavity.  As  the 
burning  rate  is  increased,  the  maximum  pressure  in  the  crack 
also  increases  because  higher  burning  rate  propellants 
generate  more  gases,  and  for  cracks  of  identical  initial  length 
and  width  (i.c.,  same  initial  volume  of  the  crack  cavity),  this 
results  in  higher  pressure. 

Figure  14  shows  the  effect  of  crack  length  on  ignition-front 
velocity.  Propellant  B  was  used  in  these  tests,  and  the  gap 
width  was  kept  constant  at  0.89  mm.  The  crack  lengths  were 
4),  130.  and  193  mm.  As  the  crack  length  is  increased,  the 
time  required  for  the  ignition  front  to  reach  the  crack  lip 
increases,  and  the  average  ignition-front  propagation  velocity 
also  increases.  The  times  required  for  reaching  the  tip  are  232, 
476,  and  311  ms,  and  the  maximum  ignition-front  velocities 
are  433,  330,  and  673  m/s  for  crack  lengths  of  63,  130,  and 
193  mm,  respectively.  Maximum  pressure  in  the  crack  also 
increases  as  the  crack  length  is  increased.  Higher  flame- 
propagation  velocities  for  longer  cracks  are  the  result  of  the 
reduced  end  effects  which  allows  the  flame  front  to  accelerate. 
The  ignition  delay  near  the  entrance  portion  of  the  crack  is 
greater,  however,  for  longer  cracks.  This  may  be  the  result  of 
crack  entrance  deformation  caused  by  pressure  exerted  on  the 
front  end  of  the  crack  sample,  in  general,  these  observations 
are  in  agreement  with  the  limited  data  reported  by  Bobolev  et 
al.'» 

In  order  to  observe  hot  gas  penetration  and  ignition-front 
propagation  processes  in  detail,  experiments  were  conducted 
using  the  alternate  crack  configuration  (see  Fig.  3).  Figure  IS 


Sw-eiMulMit  Aalal  Dtauiet,  xi l 
Fig.  14  Effect  if  ciack  teagth  oa  oiumt  IgaMoa 


float 


Moa-etOMiUiul  AaUl  DUcoaca,  x/L 
Fig.  IS  Time  history  plot  of  bet-got  front  locoUoo  on*  Ignition  front 
toceltoa. 


shows  a  typical  time  history  plot  of  the  hot  gas  front  location 
and  the  ignition  front.  The  pressurization  rate  for  this  case 
was  4000  atm/s,  and  the  crack  gap  was  1.3  mm.  A  low 
pressurization  rate  was  used  in  this  test  in  order  to  detect  the 
hot  gas  front  clearly.  It  was  observed  that  I)  the  hot  gases 
precede  the  ignition  front  along  the  crack  and  reach  the  tip 
much  earlier  than  the  convective  ignition  front;  2)  the  hot 
gases  near  the  crack  tip  are  luminous  for  a  short  period  of 
time,  with  luminosity  disappearing  later  (this  can  be  at¬ 
tributed  to  the  quenching  of  the  hot  gases);  and  3)  the  flame 
from  is  nOmmiform  near  the  erhek  entrar.ee,  but  becomes 
quite  uniform  as  it  propagates  along  the  crack.  The 
nonuniformity  of  the  name  front  near  the  crack  entrance 
region  can  be  attributed  to  the  nonuniform  flow  at  the  en¬ 
trance  of  a  rectangular  cross-sectional  crack  cavity. 
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Tht  development  of  convective  burning  in  an  isolated  crack 
has  bt«n  studied  experimentally  under  a  wide  range  of 
operating  conditions  and  propellant  geometries.  Four  types  of 
AP -based  composite  solid  propellant  were  studied.  A  solid 
propellant  igniter  system  was  developed  to  closely  simulate 
actual  rocket  conditions.  A  high-speed  motion  picture  camera 
was  used  to  measure  the  flame  propagation  rate.  Several 
important  observations  and  conclusions  from  this  study  are 
summarised  as  follows. 

1)  The  initial  pressure  distribution  in  the  crack  is  controlled 
by  the  chamber  pressurization  rate;  however,  the  maximum 
pressure  in  the  crack  cavity  is  controlled  by  the  initial  crack 
geometry  and  the  ignition  and  combustion  processes  in  the 
crack. 

2)  The  hot  gases  precede  and  reach  the  crack  tip  sooner  than 
the  ignition  front.  The  ignition-front  propagation  speed 
increases  near  the  crack  entrance,  reaches  a  maximum,  and 
then  decelerates  near  the  closed  end  of  the  crack. 

3)  As  the  chamber  pressurization  rate  or  burning  rate  of  the 
propellant  is  increased,  the  time  required  for  the  ignition  front 
to  reach  the  crack-tip  decreases  and  the  maximum  flame-front 
propagation  velocity  increases. 

4)  As  the  gap  width  of  the  crack  is  decreased,  the  time 
required  for  the  ignition  front  to  reach  the  tip  increases,  but 
the  maximum  ignition-front  propagation  velocity  also  in¬ 
creases  slightly.  As  the  crack  length  is  decreased,  both 
maximum  ignition-front  propagation  velocity  and  the  time 
required  for  ignition  front  to  reach  the  tip  decrease. 

3)  The  maximum  pressure  observed  in  the  crack  cavity 
increases  when  the  gap  width  is  reduced,  when  the  crack 
length  is  increased,  or  when  the  propellant  burning  rate  is 
increased. 

Acknowledgments 

A  portion  of  the  results  reported  in  this  paper  were  obtained 
under  the  research  program  sponsored  by  the  Aerother- 
mochcmistry  Division  of  the  Naval  Weapons  Center,  China 
Lake,  California.  The  support  and  advice  of  Dr.  R.  L.  Derr 
and  Mr.  C.  F.  Price  arc  appreciated.  Some  of  the  results 
reported  here  were  obtained  under  the  project  sponsored  by 
the  Power  Program,  Office  of  Naval  Research,  Arlington, 
Virginia,  Contract  N00014-79-C-0762.  The  support  and 
advice  of  Dr.  R.  S.  Miller  is  acknowledged.  The  assistance  of 
J.  E.  Willi  and  D.  Davies  of  The  Pennsylvania  State 
University  is  also  acknowledged. 

References 

'Bradley,  H.  H.  Jr.  and  Hoags.  T.  L.,  “Convcctivt  Burning  in 
Propellant  Defects:  A  Literature  Review,"  Naval  Weapons  Center, 
China  Lake,  Calif.,  NWC  TP  <007,  Feb.  1978. 


’Belyaev,  A.  F.,  Bobo  lev,  v.  K.,  Korotkov,  A.  I.,  Sulimov,  A.  A., 
and  Chuiko,  S.  V.,  Transition  from  Oaflagration  to  Detonation  in 
Condensed  Phases,  translated  fr«m  Russian  by  Israel  Program  of 
Scientific  Translations,  Jerusalem,  197). 

5  Kuo,  K.  K.,  Chen.  A.  T,,  and  Davis,  T.  R.,  "Convective  Burning 
in  Solid-Propellant  Cracks,"  At  A  A  Journal,  Vol.  16,  June  I97g,  pp. 
600-607. 

4 Kumar,  M.  and  Kuo,  K.  K.,  "Ignition  of  Solid  Propellant  Crack 
Tip  Under  Rapid  Pressurization,”  AM  A  Journal,  Vol.  IB,  July  I9M, 
pp.  123-13). 

’Kuo.  K.  K.,  Kumar.  M.,  Kovacic,  S.  M.,  Wills.  J.  E„  and  Chang, 
T.  Y.,  "Combustion  Processes  in  Solid  Propellant  Cracks,"  Annual 
Report  (from  The  Pennsylvania  State  University)  to  the  Naval 
Weapons  Center,  China  Lake,  Calif.,  1979. 

‘Kuo,  K.  K.,  Kovalctn,  R.  L.,  and  Ackman.  S.  J.,  "Convective 
Burning  in  Isolated  Solid  Propellant  Cracks,”  Naval  Weapons 
Center,  China  Lake,  C*!if„  NWC  TP  <049.  Feb.  1979. 

’Kuo,  K.  K.,  McCItre,  D.  R.,  Chen,  A.  T.,  and  Lucas,  F.  G., 
"Transient  Combustion  h  Solid  Propellant  Cracks,"  Naval  Weapons 
Center.  China  Lake.  Calir.,  NWC  TP  59*3.  Oct.  1977. 

'Prentice,  J.  L..  "Ftanhdown  in  Solid  Propellants,"  U.  S.  Naval 
Ordnance  Test  Station,  China  take.  Calif.,  NAVWEPS  Kept.  7964, 
NOTSTP  3009,  Dec.  1961. 

'Prentice,  J.  L.,  "Combustion  in  Solid  Propellant  Grain  Defects: 
A  Study  of  Burning  in  Single-  and  Multi-Pore  Charges,"  Naval 
Weapons  Center,  China  Lake,  Calif.,  NWC  TM  3IS2.  June  1977. 

"Margolin,  A.  D.  and  Margulis,  V.  M.,  “Penetration  of  Com¬ 
bustion  into  an  Isolated  Pore  in  an  Explosive,"  Combustion,  Ex¬ 
plosion,  and  Shock  H’aver.  Vol.  5,  No.  1, 1969.  pp.  15-16. 

"Godai,  T„  "Flame  Propagation  into  the  Crack  of  Solid 
Propellant  Grain."  AIAA  Journal,  Vol.S,  July  1970.  pp.  1322-1326. 

"Bobolev,  V.  K.,  Margolin,  A.  D.,  and  Chuiko,  S.  V.,  "The 
Mechanism  by  which  Combusion  Products  Penetrate  into  t  ho  Pores  of 
a  Charge  of  Explosive  Material,"  DokLodv  Akodtmii  Sauk  USSR, 
Vol.  162.  May  I96?,pp.  368-391. 

11  Payne,  C.  E..  "Flame  Propagation  in  Propellant  Cracks,"  Air 
Force  Rocket  Propulsion  l  aboratory,  AFRPL-TR -69-66,  1969. 

"Krasnov,  Y.  K.,  Margulis,  V.  M„  Margolin,  A.  D.,  and  Pokhil, 
P.  F..  “Rate  of  Penetration  of  Combustion  into  the  Pores  of  an 
Explosive  Charge."  Combustion,  Explosion,  and  Shock  riaew,  Vol. 
6,  No.  3.  1970,  pp.  262-265. 

’’Belyaev,  A.  f„  Korotkov,  A.  I..  Sulimov,  A.  A„  Sukoyan,  M. 
K.,  and  Obmenin.  A.  V.,  “Development  of  Combustion  in  an 
Isolated  Pore,"  Combustion.  Explosion,  and  Shock  H’at'es.  Vol.  5, 
Jan. -March  1969,  pp.  4-9. 

'* Belyaev,  A.  F„  Bobolev,  V.  K.,  Korotkov,  A.  I.,  Sulimov,  A.  A„ 
and  Chuiko,  S.  V.,  "Development  of  Burning  in  a  Single  Pore," 
Transition  of  Combustion  of  Condensed  Systems  to  Detonation, 
Chap.  5,  Pt.  A.  Sec.  22,  Science  Publisher,  1973,  pp.  1 15-134. 

”  Kumar,  M.,  "A  Study  of  Flame  Spreading  and  Combustion 
Processes  in  Solid  Propellant  Cracks,"  Ph.D.  Thesis,  The  Penn¬ 
sylvania  State  University.  Aug.  1980. 

"Taylor,  J.  W,,  "The  Burning  of  Secondary  Explosive  Powers  by 
a  Convective  Mechanism,"  Transactions  of  the  Faraday  Society,  Vol. 
58, 1962,  pp.  561-568. 


APPENDIX  E 
<Ref.  No.  7) 


A  COHPRBgglVB  IGNITION  MODEL  FOR  COMPOSITE  SOLID  PROP 


A.  K.  Kulkarni,  H.  Kumar, ^  and  K.  K.  Kuo’ 
Department  of  Mechanical  Engineering 
The  Pennsylvania  State  University 
University  Park,  PA  16802 


May  25,  1981 


(Abstfadt  SuBiitted  for  presentation  at  the  Eighteenth  JANNAF 
Combustion  Meeting,  Pasadena,  Ca, ,  Oct.  19-23,  1981.) 


*This  research  represents  a  part  of  results  obtained  under  Contract  No. 
N00014-79-C-0762 ,  sponsored  by  the  Power  Program,  Office  of  Naval  Research, 
Arlington,  Va.,  under  the  management  of  Dr.  R.  S.  Miller. 

Assistant  Professor 
s 

Associate  Professor 


149 


Abstract 

Ignition  of  a  solid  props 11 ant  is  a  coaplex  physico-chemical  phenoaenon 
which  involves  interaction  between  various  processes  such  as  heat  transfer, 
fluid  mechanics,  phase  change,  mass  diffusion  of  cheat! cal  species,  chemical 
kinetics,  etc.  The  theoretical  modeling  of  the  ignition  process  is  further 
complicated  by  such  factors  as  the  heterogeneous  nature  of  the  propellant  (and 
the  overall  ignition  process) ,  intricate  chemical  kinetics,  determination  of  the 
controlling  mechanism  of  the  ignition  process,  and  selection  of  appropriate 
ignition  criterion.  A  review  of  the  solid  propellant  ignition  was  recently 
conducted  by  the  authors.1  It  was  noted  that  the  ignition  process  under  rapidly 
increasing  pressure  conditions  which  is  typical  of  ignition  transients  in  rocket 
motors  end  the  development  of  deflagration-to-detonation  transition  (DDT)  process 
in  fractured  propellant  grains,  has  not  been  adequately  studied  in  the  past. 
Furthermore,  most  ignition  models  employ  a  number  of  simplifying  assumptions 
which  may  not  be  valid  under  actual  operating  conditions.  Nor,  has  the  influence 
of  chemical  kinetics  on  ignition  processes  been  treated  in  detail  to  date.  All 
of  the  previous  models,  except  that  of  Kumar  and  Hermahce,2  have  considered  the 
composite  propellant  to  be  one -dimensional .  Most  models  make  an  a  priori  assump¬ 
tion  about  the  site  of  ignition,  i.e. ,  in  the  solid-phase,  gas-phase,  or  at  the 
interface.  It  is  apparent  that  in  order  to  fully  understand  the  ignition  process 
and  to  accurately  predict  the  ignition  delay  time,  a  comprehensive  theoretical 
model  is  needed.  The  objective  of  this  study  is  to  develop  a  generalized  ignition 
model.  Specific  features  of  the  proposed  model  are: 

1.  Two-dimensional  (axisymnetric)  geometry,  allowing  a  more  complete 
description  of  the  heterogeneous  propellant  ignition. 

2.  Inclusion  of  a  pressurization  (dP/dt)  term  in  the  governing 
equation  for  the  gases  surrounding  the  propellant  in  order  to 
simulate  actual  rocket  motor  ignition  conditions. 

3.  No  a  priori  assumption  of  a  solid-phase,  heterogeneous,  or 
gas-phase  reaction  mechanism  to  specify  the  ignition  site. 

(The  current  model  allows  chemical  reactions  in  all  regions, 
including  the  gas  phase,  interfaces,  and  subsurface.  ) 

4.  Consideration  of  detailed  chemical  kinetics  information  in  the 
formulation  of  the  model . 

The  physical  model  considers  an  oxidizer  particle  embedded  in  a  fuel  binder 
matrix;  the  sizes  of  the  particle  and  the  surrounding  binder  are  determined 
statistically  from  the  oxidizer  particle-size  distribution  and  the  fuel/oxidant 
ratio  of  the  propellant .  The  shape  of  the  oxidizer  particle  is  approximated  by 
a  cylindrical  pellet.  The  mathematical  model  consists  of  governing  equations 
for  the  solid-phase  and  gas-phase  regions.  In  the  solid  phase,  two  energy 
equations  are  considered;  one  for  the  oxidizer  and  the  other  for  the  fuel. 

Source  terms  in  these  equations  include  contributions  due  to  indepth  radiation 
absorption,  photochemical  heat  release,  and  heat  absorption  due  to  pyrolysis. 

The  equations  are  coupled  to  the  gas-phase  conservation  equations  through  the 
heat  flux  balance  at  the  solid-gas  interface.  Gas-phase  behavior  is  described 
by  mass,  energy,  and  species  conservation  equations.  The  measured  pressure-time 
traces  near  the  propellant  sample  surface  is  considered  to  be  an  input  to  the  model. 
For  AP-based  propellants,  five  different  species  (NHj,  HCIO^,  gaseous  fuel,  gaseous 
oxidizer,  and  products)  are  considered  to  be  present  in  the  gas  phase.  Five 


ISO 


chemical  reactions  tr«  considered  for  this  feaily  of  propdlTantb,  viz., 
exothermic  degradation  of  solid  AP  particle  into  oxidizer  gases,  dissociative 
s^liaifilik.of  tf  into  NH.  end  HC10.,  reaction  between  preaixed  NH*  end  HCIO,, 
fuel  pyrolysis,  and  diffusion  flea*  resulting  from  chemical  reaction  between 
oxidibif  and  fuel  gases  to  form  products.  The  required  initial  and  boundary 
conditions  tor  the  governing  equations  are  specified  in  detail  in  the  paper. 

The  proposed  model  could  serve  as  a  framework  for  the  formulation  of 
ignition  models  for  a  wide  range  of  conditions  and  propellants,  since  portions 
of  this  model  can  be  easily  simplified  and/or  replaced  for  specific  applications. 
It  can  also  aid  in  the  identification  of  the  necessary  empirical  input  required 
for  a  given  propellant.  Numerical  solution  of  this  model  is  in  progress;  results 
wili  be  presented  at  a  later  date. 
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Investigation  of  solid  propellant  ignition  under  rapid  pres suri ration  is 
essential  to  understand  fundamentals  of  ignition  mechanism  and  deflagration-to- 
detonation  transition  (DDT)  processes.  In  soae  of  our  previous  research  conducted 
at  The  Pennsylvania  State  University,  it  was  observed  that  the  tip  region  of  a 
solid  propellant  crack  can  be  ignited  by  compression  waves  generated  by  rapid 
chamber  pressurization.  It  was  also  observed  that  a  propellant  sample  placed  at 
the  tip  of  an  inert  crack  can  be  ignited  in  the  same  Manner.  The  measured  ignition 
delay  under  these  conditions  is  in  the  order  of  1  ms  or  leas.  This  short  ignition 
delay  tisw  implies  a  high  rate  of  energy  deposition  onto  the  propellant.  The 
mechanism  involved  could  be  used  either  for  interpreting  the  evolution  of  DDT  in 
fractured  propellant  grains  or  fbr  developing  future  efficient  igniter  systems. 

The  present  investigation  deals  with  both  experimental  and  theoretical  studies 
of  solid  propellant  ignition  under  rapid  pressure  loading.  Specific  objectives 
of  this  study  are: 

1.  To  observe  the  detailed  ignition  phenomenon  and  to  measure  the  ignition 
delay  time  and  instantaneous  heat  flux  to  the  propellant  sample  surface; 

2.  To  study  the  affect  of  pressurization  rate  (dP/dt),  igniter  gas  tempera - 
ture  (Tf),  and  propellant  type  on  ignition  delay;  and 

3.  To  predict  the  ignition  delay,  using  measured  values  of  heat  flux  to 
the  propellant  sample  surface  and  to  compare  calculated  and  measured 
values . 

Test  firings  have  been  conducted  to  experimentally  determine  the  ignition  delay 
of  solid  propellants.  The  combustion  chamber  consists  of  a  solid  propellant 
igniter  system  to  generate  high-temperature ,  high-velocity  product  gases.  Product 
gases  flow  into  the  main  cavity  of  the  test  rig  with  an  exit  nozzle  at  the  opposite 
end.  A  small  sample  of  propellant  is  situated  at  the  tip  of  an  inert  CTack  channel, 
which  is  perpendicular  to  the  main- flow  direction  of  igniter  gases.  As  a  result 
of  chamber  pressurization,  hot  gases  penetrate  the  crack  and  transfer  heat  to  the 
propellant  sample.  The  ignition  event  is  observed  through  a  transparent  plexiglass 
window . 

The  data  acquisition  system  consists  of  a  high-speed  (up  to  44,000  pictures 
per  second)  movie  camera  and  photodiode  system  to  detect  ignition,  a  piezoelectric 
pressure  transducer  system  to  measure  pressure-time  history,  and  a  thin-film 
thermocouple  located  at  the  crack  tip  adjacent  to  the  propellant  sample  to  deduce 
heat  flux  to  the  propellant  surface.  Transient  signals  are  recorded  on  a  multi¬ 
channel  waveform  recorder  with  a  maximum  sampling  rate  of  2  MHz.  The  digitized 
output  of  the  waveform  recorder  is  transferred  to  a  minicomputer  foT  data  storage 
and  processing. 

In  the  theoretical  analysis,  the  composite  propellent  is  considered  to  consist 
of  A P  oxidizer  particles  surrounded  by  the  fuel  binder.  The  equivalent  particle 
sizes  are  determined  stochastically.  The  energy  transfer  to  the  propellant  is 
modeled  by  two  coupled  transient  2D  (axisymaetric)  heat  conduction  equations  for 
a  heterogeneous  control  volume  composed  of  an  oxidizer  particle  and  fuel  binder. 

The  energy  equation  includes  both  surface  and  subsurface  solid-phase  reactions,  as 


well  «a  indepth  radiation  absorption.  Tbs  host  flux  deduced  fro*  temperature 
measurements  Is  ussd  as  a  boundary  condition  for  tha  partial  differential  equa¬ 
tions.  The  numerical  solution  is  obtained  by  using  an  implicit  iterative  scheme. 


influence  of  other  parameters  on  the  ignition  process. 


154 


DISTRIBUTION  LIST 


Dr.  L.V.  Schmidt 


No.  Copies 
1 


Assistant  Secretary  of  tht  Navy 
(R.E,  and  S)  Room  5E  731 
Pentagon 

Washington*  O.C.  20350 

Or.  A.L.  Slafkosky 
Sclantlflc  Advisor 
Commandant  of  the  Marine  Corps 
Code  RD-1 

Washington*  O.C.  20380 

Dr.  Richard  S.  Miller 
Office  of  Naval  Research 
Code  413 

Arlington,  VA  22217 

Mr.  David  Siegel 
Office  of  Naval  Research 
Code  260 

Arlington,  VA  22217 

Or.  R.J.  Marcus 
Office  of  Naval  Research 
Western  Office 
1030  East  Green  Street 
Pasadena,  CA  91106 

Or.  Larry  Peebles 
Office  of  Naval  Research 
East  Central  Regional  Office 
666  Sumner  Street,  Bldg.  114-D 
Boston,  MA  022IQ 

Or.  Phillip  A,  Miller 
Office  of  Naval  Research 
San  Francisco  Area  Office 
One  Hal  1  Idle  Plaza,  Suite  501 
San  Francisco,  CA  94102 

Mr.  Otto  K.  Helney 
AFATL  -  OLOL 
Elgin  AFB,  FL  32542 

Mr.  R,  Geisler 
ATTN:  MKP/MS24 
AFRPL 

Edwards  AFB,  CA  93523 


Dr.  F.  Roberto 
Code  AFRPL  MKPA 
Edwards  AFB,  CA  93523 


No.  Copies 
1 


Or.  L.H.  Caveny  1 

Air  Force  Office  of  Scientific 
Research 

Directorate  of  Aerospace  Sciences 
Bolling  Air  Force  Base 
Washington,  O.C.  20332 

Mr.  Donald  L.  Ball  1 

Air  Force  Office  of  Scientific 
Research 

Directorate  of  Chemical  Sciences 
Bolling  Air  Force  Base 
Washington,  O.C.  20332 

Or.  John  S.  Wilkes,  Jr.  1 

FJSRL/NC 

USAF  Academy,  CO  3o8Uo 

Or.  R.L.  Lou  1 

Aerojet  Strategic  Propulsion  Co. 

P.0.  Box  15699C 
Sacramento,  CA  95813 

Or.  V.j.  Keenan  1 

Anal-Syn  Lab  Inc. 

P.0.  Box  547 
Paoli,  PA  19301 

Or.  Philip  Howe  T 

Army  Ballistic  Research  Labs 
ARRADCOM 
Code  DROAR-BLT 

Aberdeen  Proving  Ground,  MO  21C05 

Mr.  L.A.  Watermeier  1 

Amy  Ballistic  Research  Labs 
ARRADCOM 
Code  OROAR-SLI 

Aberdeen  Proving  Ground,  MO  21005 

Or.  W.W.  Wharton  1 

Attn:  ORSMI-RKL 

Commander 

U.S.  Army  Missile  Command 
Redstone  Arsenal ,  AL  35898 


155 


OYti 


6/31 


DISTRIBUTION  LIST 


Or.  R.G.  Rhoades,  , 

Commander 

Artnv  Missile1  Command 
ORSMI-R 

Redstone  Arsenal ,  AL  35898 

Or.  W.O.  Stephens 
Atlantic  Research  Coro. 

Pine  Ridge  Plant 
7511  Wellington  Rd. 

Gainesville  VA  22065 

Dr.  A.W.  Barrows 
Ballistic  Research  Laboratory 
USA  ARRAOCOM 
OROAR-BLP 

Aberdeen  Proving  Ground,  MD  21005 

Or.  C.M.  Frey 
Chemical  Systems  Division 
P.O.  8ox  358 
Sunnyvale,  CA  94086 

Professor  F.  Rodriguez 
Cornell  University 
School  of  Chemical  Engineering 
01  in  Hall,  Ithaca,  N.Y.  14853 

Defense  Technical  Information 

Center 

DTIC-DOA-2 

Cameron  Station 

Alexandria,  VA  22314 

Dr.  Rocco  C.  Musso 
Hercules  Aerospace  Division 
Hercules  Incorporated 
Alleghany  Ballistic  Lab 
P.O.  Box  210 
Washington,  D.C.  21502 

Or.  Ronald  L.  Simmons 
Hercules  Inc.  Eglin 
AFATL/OLDL 

EglinAFB,  FL  32542 


No.  Copies 
1 


12 


No,  Copies 

Dr.  E.H.  Debutts  1 

Hercules  Inc. 

Baccus  Works 
P.O.  Box  98 
Magna ,  UT  84044 

Dr.  James  H.  Thacher  1 

Hercules  Inc.  Magna 

Baccus  Works 

P.O.  Box  98 

Magna,  UT  84044 

Mr.  Theordore  M.  Gilliland  1 

Johns  Hopkins  University  API 
Chemical  Propulsion  Info.  Agency 
Johns  Hopkins  Road 
Laurel,  MO  20810 

Dr.  R.  McGuire  1 

Lawrence  Livermore  Laboratory 
University  of  California 
Code  L-324 

Livermore,  CA  94550 

Dr.  Jack  Linsk  1 

Lockheed  Missiles  &  Space  Co. 

P.O.  Box  504 

154 

Dr.  B.G.  Craig  1 

Los  Alamos  National  Lab 

P.O.  Box  1663 

NSP/DOD ,  MS-245 

Los  Alamos,  Nfl  87545 


Dr.  R.L.  Rabie 

WX-2,  MS-952 

Los  Alamos  National  Lab. 

P.O.  Box  1663 

Los  Alamos  NM  37545 

E!S$  SlartSi^ier^Tf ic  Lab. 

P.O.  Box  1663 

Los  Alamos,  NM  37545 


1 


'■  77 . “•  -r'ZZ- 


I 


.  OYN 


DISTRIBUTION  L'.tST’ 


'I-  I,  ■: 


No.  Copies 


Nr.  R.  Brown 

Naval  Air  Systems  Command 
Code  330 

Washington,  D.C.  20361 

Or.  H.  Rosenwasser 
Naval  Air  Systems  Command 
AIR-31 OC 

Washington,  D.C.  20360 

Mr.  8.  Sobers 

Naval  Air  Systems  Command 

Cods  03P25 

Washington,  O.C.  20360 

Or.  L.R.  Rothsteln 
Assistant  Director  : 

Naval  Explosives  Cev, 
Engineering  Dept.  ' 

Naval  Weapons  Station. 
Yorktown,  VA  23691 

Or.  Lionel  Dickinson 
Naval  Explosive  Ordnance 
Oisposal  Tech.  Center 
Code  D 

Indian  Head,  MO  20640 

Mr.  C.L.  Adams 
Naval  Ordnance  Station 
Code  PM4 

Indian  Head,  MD  20640 

Mr.  S.  Mitchell 
Naval  Ordnance  Station 
Cede  5253 

Indian  Head,  MD  20640 

Dr.  William  Tolies 
Oean  of  Research 
Naval  Postgraduate  School 
Monterey,  CA  93940 

■Naval  Research  Lab.. 

Code. 61  CO 

Washi ng  ton ,  O.C.  20375 


No,  Copies 

Dr.  J.  Schnur  .1  >' 

Naval  Research  Lab. 

Code  6910  ' 

Washington,  D.C.  20375 

Mr,  R.  Beauregard  ,  1 

Naval  Sea  Systems  Command 
SEA  64E 

Washington,  D.C.  20362 

Mr.  G.  Edwards  1 

Naval  Sea  Systems  Command 
Code  62R3 

Washington,  D.C.  20362 

Mr.  John  Boyle  ,  1 

Materials  Branch 
Naval  Ship  Engineering  Center 
Philadelphia,  PA  19112 

Dr.  H.G.  Adolph  1 

Naval  Surface  Weapons  Center 
Cude  Rll 
White  Oak 

Silver  Spring,  MD  20910 

Dr.  T.Dk  Austin  1 

Naval  Surface  Weapons  Center 
Code  R16 

Indian  Head,  MD  20640 

Or.  T.  Hall  1 

Code  R-11 

Naval  Surface  Weapons  Center 
White  Oak  Laboratory 
Silver  Spring,  MD  20910 

Mr.  G.L.  Mackenzie  1 

Naval  Surface  Weapons  Center 
Code  R101 

Indian  Head,  MD  20640 

Dr.  K.F .  Mueller  1 

Naval  Surface  'Weapons  Center 
Code  R1 1 
White  Oak 

Silver  Spring,  MD  20910 


1 


f  • ,  • ,  157 

6/31  . 


DISTRIBUTION  LIST 


Mr.  J.  Murrin  1 

Naval  Sea  Systems  Command 
Code  62R2  “ 

Washington,  D.C.  20362 


Dr.  A.  Nielsen-  1 

Naval  Weapons  Center 
Code  385 

China  Lake,  CA  93555 


Qr.  D.J.  Pastine.  1 

vNaval  Surface  Weapons  Cneter 
Code  R04 
White  Oak 

Silver  Spring,  MO  20910 

Mr.  L.  Roslurtd  1 

Naval  Surface  Weapons  Center 
Code  R122 

White  Oak,  Silver  Spring 
MD  20910 

Mr.  M.  Stosz  1 

Naval  Surface  Weapons  Center 
Code  R121 
White  Oak 

Silver  Spring,  MO  20910 

Dr.  E.  Zimmet  1 

Naval  Surface  Weapons  Center 
Code  R13 
White  Oak 

Silver  Spring,  MD  20910 

Or.  0.  R.  Oerr  1 

Naval  Weapons  Center 
Code  388 

China  Lake,  CA  93555 

Mr.  Lee  N.  Gilbert  1 

Naval  Weapons  Center 
Code  3205 

China  Lake,  CA  93555 

Dr.  E.  Martin  1 

Naval  Weapons  Center 
Code  3858 

China  Lake,  CA  93555 

Mr.  R.  McCarten  1 

Naval  Weapons  Center 
Code  3272 

China  Lake,  CA  93555 


Or.  R.  Reed,  Jr.  1 

Naval  Weapons  Center 
Code  388  ;  ; 

China  Lake,  CA  93555 

Dr.  L.  Smith  1 

Naval  Weapons  Center 
Code  3205 

China  Lake,  CA  93555 

Or.  B;  Douda  " 

Naval.  Weapons  Support  Center 
Code  5042 

Crajne,  Indiana  47C22 

Or.  A.  Faulstich  1 

Chief  of  Naval  Technology 
MAT  Code  0716  /' 

Washington,  O.C.  20360 

LCOR  J.  Walker  1 

Chief  of  Naval  Material 
Office  of  Naval  Technoloay 
MAT,  Code  0712 
1  Washingf.cn  ,  O.C .  2C3bO 

Mr.  ooe  McCartney  1 

Naval  Ocean  Systems  Center 
San  01  ego,  CA  92152 

Or.  Yamamoto  1 

Marine  Sciences  (Division 
Nava!  Ocean  Systems  Center 
San  Diego,  CA  91.232 

Or.  G.  3osmajiai-!  1 

Applied  Chemistry.  Division  .’ 
Naval  Ship  Research  l  Development 
Center 

Annual  is,  MD  21401 

Or.  H.  Shuey  1 

Rohn  and  Ha*s  Company 
••Huntsv.fi-3s,  Alabama  3660V 


V:  <  < 

t  i  •' 

Dr.  J.F. 

Kincaid 

§ 

•  Str-stagic 
Office  . 

-  Systems 

Proje 

1!  ■  < 

tv 

AO*  Off  tmijfi 

t  of  the 

Mcvy, 

'  "  n.ooo  Bor' 

if.  .: 

•i'fishir.gto 

r'm'O.C. 

20:75 

tonic  Systems  Project  Office  1 
Pro.^i  5 ion  yM't  ■„  ; 

Ciis  $02701  i  .  *'!  ■■  •'  '  , 
r-c psrvr  of  'the  'la'/y 
We chino tori.,  O.Ck'  $0376 

Tlv^ocl^ortCrt  l 

Strategic  Systems;  Preset  Office 
Department  of  the  Ma’vy 
Roe;"  1  ;"•*  3  • 

,  , ..Wasihyi'jton'v  O.C.  2037,6 

0-*.  p . A ,  ftaniton  '  1 

1  h  :.i 

nun :>/ 1 .lie  Division., 
huotcviiVh ,  41 .\baoa  35807  ■.• 

Mr.  6. F.  Man, 1 
1 h ’»  :.ol  Oorco, ratten 
'  lAntSville  D;.i '*i Sion 
Ho'KC'Vil.'i?,  Alabama  25307 

.  ”r.  E.S.  S'.itt&ft , ’  1 

Thicko"!  eviration 
.  ElUor  Division  , 

P.0.  Box  24  i 
.  'Slkfon,.  K0  21 9?1 


Dr.  3,  ,T 

ruickoi 

Wasatch 

Division 

!■■:  34')  P 

.0.  to*  .524 

Di  ighan 

City, 'Ur.  -3d  302 

Or.  TiF. 

Davidson .  " 

Jf'  ;-oicj1  Di  •••so tor 


Sfititrl 


Gcwormsliv  ;■■:■'•  “2:;'?  Grc.»;V 
»  A  !?.-■• 

i  ,'i .  -w.j/.,  ■ 

.  Odocr  -r-  Vt.ii  04403 


Or.  C.Vi.  Vrisseri 
Thiokol  Elk  ton  Division 
P.0.  Be/.  241 
El kton ,  ;M0  21921 

Dr.  J .,C.  •  Hinsh.*>i.  .1 

Thicko  1  l.'ssotch  Divis'or. 

P.0.  Sc; ,524 

Brighan  City,  Utah,.-  33402 

U.S.  Arrr/  Research  Office  1 

C  hen  i  c  1*  U  Biological  Sc  i  er.  c  o  s 

Division 

P.0.  3-j"  12211 

Research  Triangle  Park  •. 

MC  27703' 

Dr.  R.P.  Walker  '  1 

USA  ARAADCOM  . 

DRDAR-LOC  , 

Dover,  :,’J  07CG1 

Or.  T,  Siroicn  1 

Muritti cr.s  Di ..rectors te 
Pror-el  1  ants  anil .  C/.rl  soiVes 
Defence  “.guipient  .Staff 
British  Embassy 
3100’  MaSsacvotatls  Ave. 

Washing  ten,  O.C.  20003. 

LTC  3.  Loving  i 

AFROl/UC 

Edwards  AFg,  CA  93523 

Professor  Alan  N.  Gent  1 

Institute  of  Polymer  Science 
University  of  Akron 
Akron,  OH  44325 

Mr.  J.  M,  Frankie  1 

Army  Ballistic  Research  Labs 
ARRAOCOM 
Code  DP.0AR-8LI 

Aberdeen  Proving  Ground,  MO  21005 


■rfjwwag'j*  vaisay;asaihia^^aaag^a  MBMiiif  imunar-ai£«i'-ra«Wa^ 


6/81 


DISTRIBUTION  LIST 


No.  Copies 


Dr.  Ingo  W.  May 

Army  Ballistic  Research  Labs 

ARRADCOM 

Code  DROAR-BLI 

Aberdeen  Proving  Ground,  MD  21005 


Professor  N.W.  Tschoegl 
California  Institute  of  Tech 
Dept,  of  Chemical  Engineering 
Pasadena,  CA  91125 


Professor  M,D.  Nicol 
University  of  California 
Oept.  of  Chemistry 
405  Hllgard  Avenue 
Los  Angeles,  CA  90024 


Professor  A.  G.  Evans 
University  of  California 
Berkeley,  CA  94720 


Professor  T.  Litovitz 
Catholic  Univ.  of  America 
Physics  Department 
520  Michigan  Ave.,  N.E. 
Washington,  D.C.  20017 


Professor  W*  G.  Knauss 
Graduate  Aeronautical  Lab 
California  Institute  of  Tech. 
Pasadena,  CA  31125 


Professor  Edward  Price 
Georgia  Institute  of  Tech. 
School  of  Aerospace  Engin. 
Atlanta,  Georgia  30332 


Dr.  Kenneth  0.  Hartman 
Hercules  Aerospace  Division 
Hercules  Incoroorated 
P.O.  8ox  210 
Cumberland,  MD  21502 


Or.  Thor  L.  Smith 
IBM  Research  Lab 
042.282 

San  Jose,  CA  95193 


No.  Copies 


Or.  J,  P.  Marshall 
Oept.  52-35,  Bldg.  204/2 
Lockheed  Missile  &  Space  Co. 
3251  Hanover  Street 
Palo  Alto,  CA  94304 


Ms.  Joan  L.  Janney 
Los  Alamos  National  Lab 
Mail  Stop  920 
Los  Alamos,  NM  37545 


Dr.  J.  M.  Walsh  1 

Los  Alamos  Scientific  Lab 
Los  Alamos,  NM  87545 


Professor  R,  W.  Armstrong 
Univ.  of  Maryland 
Department  of  Mechanical  Eng. 
College  Park,  MD  20742 


Prof.  Richard  A.  Reinhardt  1 
Naval  Postgraduate  School 
Physics  &  Chemistry  Dept. 
Monterey,  CA  93940 


Dr.  R.  Bernecker  1 

Naval  Surface  Weapons  Center 
Code  R13 

White  Oak,  Silver  Spring,  MO  20910 


Dr.  M,  J;  Kamlet  1 

Naval  Surface  Weapons  Center 
Code  Rll 

White  Oak,  Silver  Soring,  MD  20910 


Professor  J.  D.  Achenbach 
Northwestern  University 
Dept,  of  Civil  Engineering 
Evanston,  IL  60201 


Dr.  N,  L.  8asdekas 
Office  of  Naval  Research 
Mechanics  Program,  Code  432 
Arlington,  V A  22217 


Professor  Kenneth  Kuo  1 

Pennsylvania  State  Univ. 

Dept,  of  Mechanical  Engineering 
University  Park,  PA  1 6302 


160 


OYN 


6/81 


DISTRIBUTION  LIST 


No.  Copies 


No.  Copies 


Or.  S.  Sheffield  1 

Sandia  Laboratories 

Division  2513 

P.O.  Box  5800 

Albuquerque,  NM  87185 

Dr.  M.  Farber  1 

Space  Sciences,  Inc. 

135  Maple  Avenue 
Monrovia,  CA  91016 

Dr.  Y.  M.  Gupta  1 

SRI  International 
333  Ravenswood  AVenue 
Menlo  Park,  CA  94025 

Mr.  M.  Hill  1 

SRI  International 
333  Ravenswood  Avenue 
Menlo  Park,  CA  94025 

Professor  Richard  A.  Schapery  1 

Texas  A&M  Univ. 

Dept  of  Civil  Engineering 
College  Station,  TX  77843 

Dr.  Stephen  Swanson  1 

Univ.  of  Utah 

Dept,  of  Mech.  &  Industrial 
Engineering 
[  MEB  3008 

j  Salt  Lake  City,  UT  34112 

I  Mr.  J,  D.  Byrd  1 

•  Thiokol  Corp.  Huntsville 

I  Huntsville  Div. 

|  Huntsville,  AL  35807 

I  Professor  G.  D.  Duvall  1 

|  Wasnington  State  University 

[  Dept,  of  Physics 

[  Pullman,  WA  99153 

Prof.  T.  Dickinson 
Washington  State  University 
Dept,  of  Physics 
Pullman,  WA  99163 


1 


